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CONTROLLERS, OBSERVERS, AND 
APPLICATIONS THEREOF 

CROSS REFERENCE TO RELATED 

APPLICATIONS 5 

This application is a continuation-in-part (CIP) application 
and claims the benefit of U.S. Utility application Ser. No. 
10/351,664 filed Jan. 27, 2003 published as 2003/0199997. 
This application also claims the benefit of U.S. Provisional 10 
Application 60/373,404 filed Apr. 8, 2002; U.S. Provisional 
Application 60/718,393 filed Sep. 19, 2005; U.S. Provisional 
Application 60/718,581 filed Sep. 19, 2005; U.S. Provisional 
Application 60/718,899 filed Sep. 20, 2005; U.S. Provisional 
Application 60/728,928 filed Oct. 20, 2005; and U.S. Provi- ^ 
sional Application 60/728,929 filed Oct. 20, 2005, all of 
which are incorporated herein by reference. 

This work was supported at least in part by NASA under 
government contract number NGT3-52387. Accordingly, the 
U.S. government may have certain rights herein. 20 

TECHNICAL FIELD 

The systems, methods, application programming inter- 
faces (API), graphical user interfaces (GUI), computer read- 25 
able media, and so on described herein relate generally to 
controllers and more particularly to scaling and parameteriz- 
ing controllers, and the use of observers and tracking which 
facilitates improving controller design, tuning, and optimiz- 
ing. 30 

BACKGROUND 

A feedback (closed-loop) control system 10, as shown in 
Prior Art FIG. 1, is widely used to modify the behavior of a 35 
physical process, denoted as the plant 110, so it behaves in a 
specific desirable way over time. For example, it may be 
desirable to maintain the speed of a car on a highway as close 
as possible to 60 miles per hour in spite of possible hills or 
adverse wind; or it may be desirable to have an aircraft follow 40 
a desired altitude, heading and velocity profile independently 
of wind gusts; or it may be desirable to have the temperature 
and pressure in a reactor vessel in a chemical process plant 
maintained at desired levels. All these are being accom- 
plished today by using feedback control, and the above are 45 
examples of what automatic control systems are designed to 
do, without human intervention. 

The key component in a feedback control system is the 
controller 120, which determines the difference between the 
output of the plant 110, (e.g. the temperature) and its desired 50 
value and produces a corresponding control signal u (e.g., 
turning a heater on or off). The goal of controller design is 
usually to make this difference as small as possible as soon as 
possible. Today, controllers are employed in a large number 
of industrial control applications and in areas like robotics, 55 
aeronautics, astronautics, motors, motion control, thermal 
control, and so on. 

Classic Controllers 

Classic Control Theory provides a number of techniques an 
engineer can use in controller design. Existing controllers for 60 
linear, time invariant, and single-input single output plants 
can be categorized into three forms: the proportional/integral/ 
derivative (PID) controllers, transfer function based (TFB) 
controllers, and state feedback (SF) controllers. The PID 
controller is defined by the equation 65 

( 1 ) 


2 

where u is the control signal and e is the error between the set 
point and the process output being controlled. This type of 
controller has been employed in engineering and other appli- 
cations since the early 1920s. It is an error based controller 
that does not require an explicit mathematical model of the 
plant. The TFB controller is given in the form of 


U(s) = G c (s)E(s), G c (s) = -j-^r 


where U(s) and E(s) are Laplace Transforms of u and e 
defined above, and n(s) and d(s) are polynomials in s. The 
TFB controller can be designed using methods in control 
theory based on the transfer function model of the plant, 
Gp(s). A PID controller can be considered a special case of a 
TFB controller because it has an equivalent transfer function 
of 


ki (3) 

GAs) — kn 4 " — 4 - krfS 
s 

The State Feedback (SF) Controller 

The SF controller can be defined by 

u=r+Kx (4) 

and is based on the state space model of the plant: 

x(t)=Ax(t)+Bu(t) ,y(t)=Cx(t) +Du(t) (5) 

When the state x is not accessible, a state observer (SO): 

x =Ax+Bu+L iy-y) (6) 

is often used to find its estimate, x. Here r is the set point for 
the output to follow. 

Controller Tuning 

Over the years, the advances in control theory provided a 
number of useful analysis and design tools. As a result, con- 
troller design moved from empirical methods (e.g., ad hoc 
tuning via Ziegler and Nichols tuning tables for PID) to 
analytical methods (e.g., pole placement). The frequency 
response method (Bode and Nyquist plots) also facilitated 
analytical control design. 

Conventionally, controllers are individually designed 
according to design criteria and then individually tuned until 
they exhibit an acceptable performance. Practicing engineers 
may design controllers, (e.g., PID) using look-up tables and 
then tune the controllers using trial and error techniques. But 
each controller is typically individually designed, tuned, and 
tested. 

Tuning controllers has perplexed engineers. Controllers 
that are developed based on a mathematical model of the plant 
usually need their parameters to be adjusted, or “tuned” as 
they are implemented in hardware and tested. This is because 
the mathematical model often does not accurately reflect the 
dynamics of the plant. Determining appropriate control 
parameters under such circumstances is often problematic, 
leading to control solutions that are functional but ill-tuned, 
yielding lost performance and wasted control energy. 

Additionally, and/or alternatively, engineers design using 
analytical (e.g., pole placement) techniques, but once again 
tune with trial and error techniques. Since many industrial 
machines and engineering applications are built to be inher- 
ently stable, acceptable controllers can be designed and tuned 
using these conventional techniques, however, acceptable 
performance may not approach optimal performance. 


u =K P e+K I $e+K D e 



US 8,060,340 B2 


3 

One example conventional technique for designing a PID 
controller included obtaining an open-loop response and 
determining what, if anything, needed to be improved. By 
way of illustration, the designer would build a candidate 
system with a feedback loop, guess the initial values of the 
three gains (e.g., kp, kd, ki) in PID and observe the perfor- 
mance in terms of rise time, steady state error and so on. Then, 
the designer might modify the proportional gain to improve 
rise time. Similarly, the designer might add or modify a 
derivative controller to improve overshoot and an integral 
controller to eliminate steady state error. Each component 
would have its own gain that would be individually tuned. 
Thus, conventional designers often faced choosing three 
components in a PID controller and individually tuning each 
component. Furthermore, there could be many more param- 
eters that the design engineer must tune if a TFB or a state 
feedback state observer (SFSOB) controller is employed. 

Another observation of control design is that it is not por- 
table. That is, each control problem is solved individually and 
its solution cannot be easily modified for another control 
problem. This means that the tedious design and tuning pro- 
cess must be repeated for each control problem. The use of 
state observers is useful in not only system monitoring and 
regulation but also detecting as well as identifying failures in 
dynamical systems. Since almost all observer designs are 
based on the mathematical model of the plant, the presence of 
disturbances, dynamic uncertainties, and nonlinearities pose 
great challenges in practical applications. Toward this end, 
the high-performance robust observer design problem has 
been topic of considerable interest recently, and several 
advanced observer designs have been proposed. Although 
satisfactory in certain respects, a need remains for an 
improved strategy for an observer and incorporation and use 
of such in a control system. 

State Observers 

Observers extract real-time information of a plant’s inter- 
nal state from its input-output data. The observer usually 
presumes precise model information of the plant, since per- 
formance is largely based on its mathematical accuracy. 
Closed loop controllers require both types of information. 
This relationship is depicted in 3200 of FIG. 32. Such pre- 
sumptions, however, often make the method impractical in 
engineering applications, since the challenge for industry 
remains in constructing these models as part of the design 
process. Another level of complexity is added when gain 
scheduling and adaptive techniques are used to deal with 
nonlinearity and time variance, respectively. 

Disturbance Estimation Observes and Disturbance Rejection 

Recently, disturbance rejection techniques have been used 
to account for uncertainties in the real world and successfully 
control complex nonlinear systems. The premise is to solve 
the problem of model accuracy in reverse by modeling a 
system with an equivalent input disturbance d that represents 
any difference between the actual plant P and a derived/ 
selected model P„ of the plant, including external distur- 
bances w. An observer is then designed to estimate the dis- 
turbance in real time and provide feedback to cancel it. As a 
result, the augmented system acts like the model P„ at low 
frequencies, making the system accurate to P„ and allowing a 
controller to be designed for P„. This concept is illustrated in 
3900 of FIG. 39. 

The most common of these techniques is the disturbance 
observer (DOB) structure (Endo, S., H. Kobayashi, C. J. 
Kempf, S. Kobayashi, M. Tomizuka and Y. Hori (1996). 
“Robust Digital Tracking Controller Design for High-Speed 
Positioning Systems.” Control Eng. Practice, 4:4, 527-536; 
Kim, B. K., H.-T. Choi, W. K. Chung and I. H. Suh (2002). 


4 

“Analysis and Design of Robust Motion Controllers in the 
Unified Framework.” J of Dynamic Systems, Measurement, 
and Control, 124, 313-321; Lee, H. S. and M. Tomizuka 
(1996). “Robust Motion Controller Design for High-Accu- 
5 racy Positioning Systems.” IEEE Trans. Ind. Electron., 43 : 1 , 
48-55; Tesfaye, A., H. S. Lee and M. Tomizuka (2000). “A 
Sensitivity Optimization Approach to Design of a Distur- 
bance Observer in Digital Motion Control.” IEEE/ASME 
Trans, on Mechatronics, 5:1, 32-38; Umeno, T. and Y. Hori 
to (1991). “Robust Speed Control of DC Servomotors Using 
Modern Two Degrees of Freedom Controller Design”. IEEE 
Trans. Ind. Electron., 38:5, 363-368). It uses simple binomial 
Q-filters, allowing the observer to be parameterized, i.e. tuned 
by a single bandwidth parameter. A model deliberately dif- 
1 5 ferent from P is also suggested in E. Schrij ver and J. van Dijk, 
“Disturbance Observers for Rigid Mechanical Systems: 
Equivalence, Stability, and Design,” Journal of Dynamic Sys- 
tems, Measurement, and Control, vol. 124, no. 4, pp. 539-548, 
2002 to facilitate design, but no guidelines are given other 
20 than it should be as simple as possible, cautioning stability 
and performance may be in danger. Another obstacle is that a 
separate observer must be designed to provide state feedback 
to the controller. In existing research, derivative approximates 
are used in this way but their effect on performance and 
25 stability has yet to be analyzed. Furthermore, the controller 
design is dependent on the DOB design, meaning that deriva- 
tive approximates can not be arbitrarily selected. 

Multiple DOBs were used to control a multivariable robot 
by treating it as a set of decoupled single-input single-output 
30 (SI SO) systems, each with disturbances that included the 
coupled dynamics (Bickel, R. and M. Tomizuka (1999). “Pas- 
sivity-Based Versus Disturbance Observer Based Robot Con- 
trol: Equivalence and Stability.”/, of Dynamic Systems, Mea- 
surement, and Control, 121, 41-47; Hori, Y., K. Shimura and 
35 M. Tomizuka (1992). “Position/Force Control of Multi-Axis 
Robot Manipulator Based on the TDOF Robust Servo Con- 
troller For Each Joint.” Proc. of ACC, 753-757; Kwon, S. J. 
and W. K. Chung (2002). “Robust Performance of the Mul- 
tiloop Perturbation Compensator.” IEEE/ASME Trans. 
40 Mechatronics, 7:2, 190-200; Schrij ver, E. and J. Van Dijk 
(2002) Disturbance Observers for Rigid Mechanical Sys- 
tems: Equivalence, Stability, and Design.” J. of Dynamic Sys- 
tems, Measurement, and Control, 124, 539-548. 

Another technique, referred to as the unknown input 
45 observer (UIO), estimates the states of both the plant and the 
disturbance by augmenting a linear plant model with a linear 
disturbance model (Burl, J. B. (1999). Linear Optimal Con- 
trol, pp. 308-314. Addison Wesley Longman, Inc., Calif.; 
Franklin, G. F., J. D. Powell andM. Workman (1998). Digital 
50 Control of Dynamic Systems , Third Edition, Addison Wesley 
Longman, Calif.; Johnson, C. D. (1971). “Accommodation of 
External Disturbances in Linear Regulator and Servomecha- 
nism Problems.” IEEE Trans. Automatic Control , AC-1 6:6, 
635-644; Liu, C.-S., and H. Peng (2002). “Inverse-Dynamics 
55 Based State and Disturbance Observer for Linear Time-In- 
variant Systems.” J. of Dynamic Systems, Measurement, and 
Control, 124, 375-381; Profeta, J. A. Ill, W. G. Vogt andM. H. 
Mickle (1990). “Disturbance Estimation and Compensation 
in Linear Systems.” IEEE Trans. Aerospace and Electronic 
60 Systems, 26:2, 225-231; Schrij ver, E. and J. van Dijk (2002) 
“Disturbance Observers for Rigid Mechanical Systems: 
Equivalence, Stability, and Design.” J. of Dynamic Systems, 
Measurement, and Control, 124, 539-548). Unlike the DOB 
structure, the controller and observer can be designed inde- 
65 pendently, like a Luenberger observer. However, it still relies 
on a good mathematical model and a design procedure to 
determine observer gains. An external disturbance w is gen- 
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erally modeled using cascaded integrators (1/s*). When they 
are assumed to be piece-wise constant, the observer is simply 
extended by one state and still demonstrates a high degree of 
performance. 

Extended State Observer (ESO) 5 

In this regard, the extended state observer (ESO) is quite 
different. Originally proposed by Han, J. (1999). “Nonlinear 
Design Methods for Control Systems.” Proc. 14 th IFAC 
World Congress , in the form of a nonlinear UIO and later 
simplified to a linear version with one tuning parameter by to 
Gao, Z. (2003). “Scaling and Parameterization Based Con- 
troller Tuning.” Proc. of ACC, 4989-4996, the ESO combines 
the state and disturbance estimation power of a UIO with the 
tuning simplicity of a DOB. One finds a decisive shift in the 
underlying design concept as well. The traditional observer is 15 
based on a linear time-invariant model that often describes a 
nonlinear time- varying process. Although the DOB and UIO 
reject input disturbances for such nominal plants, they leave 
the question of dynamic uncertainty mostly unanswered in 
direct form. The ESO, on the other hand, addresses both 20 
issues in one simple framework by formulating the simplest 
possible design model P rf =l/s w for a large class of uncertain 
systems. P^ is selected to simplify controller and observer 
design, forcing P to behave like it at low frequencies rather 
than P w . As a result, the effects of most plant dynamics and 25 
external disturbances are concentrated into a single unknown 
quantity. The ESO estimates this quantity along with deriva- 
tives of the output, giving way to the straightforward design of 
a high performance controller. 

Active Disturbance Rejection Control (ADRC) 30 

Originally proposed by Han, J. (1999). “Nonlinear Design 
Methods for Control Systems.” Proc. 14 th IFAC World Con- 
gress , a nonlinear, non-parameterized active disturbance 
rejection control (ADRC) is a method that uses an ESO. A 
linear version of the ADRC controller and ESO were param- 35 
eterized for transparent tuning by Gao, Z. (2003). “Scaling 
and Parameterization Based Controller Tuning.” Proc. of 
ACC, 4989-4996. Its practical usefulness is seen in a number 
of benchmark applications already implemented throughout 
industry with promising results (Gao, Z., S. Hu and F. Jiang 40 
(2001). “A Novel Motion Control Design Approach Based on 
Active Disturbance Rejection.” Proc. of 40 th IEEE Confer- 
ence on Decision and Control ; Goforth, F. (2004). “On 
Motion Control Design and Tuning Techniques.” Proc. of 
ACC; Hu, S. (2001). “On High Performance Servo Control 45 
Solutions for Hard Disk Drive.” Doctoral Dissertation , 
Department of Electrical and Computer Engineering, Cleve- 
land State University; Hou, Y., Z. Gao, F. Jiang and B. T. 
Boulter (2001). “Active Disturbance Rejection Control for 
Web Tension Regulation.” Proc. of 40 th IEEE Conf. on Deci- 50 
sion and Control ; Huang, Y., K. Xu and J. Han (2001). “Flight 
Control Design Using Extended State Observer and Nons- 
mooth Feedback.” Proc. of 40th IEEE Conf. on Decision and 
Control; Sun, B and Z. Gao (2004). “A DSP-Based Active 
Disturbance Rejection Control Design for a 1 KW H-Bridge 55 
DC-DC Power Converter.” To appear in: IEEE Trans, on Ind. 
Electronics ; Xia, Y., L. Wu, K. Xu, and J. Han (2004). “Active 
Disturbance Rejection Control for Uncertain Multivariable 
Systems With Time-Delay, 2004 Chinese Control Confer- 
ence). It was also applied to a fairly complex multivariable 60 
aircraft control problem (Huang, Y., K. Xu and J. Han (2001). 
“Flight Control Design Using Extended State Observer and 
Nonsmooth Feedback.” Proc. of 40th IEEE Conf. on Decision 
and Control). 

What is needed is a control framework for application to 65 
systems throughout industry that are complex and largely 
unknown to the personnel often responsible for controlling 


6 

them. In the absence of required expertise, less tuning param- 
eters are needed than current approaches, such as multi-loop 
PID, while maintaining or even improving performance and 
robustness. 

Linear Active Disturbance Rejection Controller (LADRC) 

In addition to the above controllers, a more practical con- 
troller is the recently developed from Active Disturbance 
Rejection Controller (ADRC). Its linear form (LADRC) for a 
second order plant is introduced below as an illustration. The 
unique distinction of ADRC is that it is largely independent of 
the mathematical model of the plant and is therefore better 
than most controllers in performance and robustness in prac- 
tical applications. 

Consider an example of controlling a second order plant 

y=-ay-by +w+bu (7) 

where y and u are output and input, respectively, and w is an 
input disturbance. Here both parameters, a and b, are 
unknown, although there is some knowledge ofb, (e.g., b 0 «b, 
derived from the initial acceleration of y in step response). 
Rewrite (7) as 

y=-ay- by+ w+(b-b 0 )u+b 0 u=f+b 0 u (8) 

where f=-ay-by+w+(b-b 0 )u. Here f is referred to as the 
generalized disturbance, or disturbance, because it represents 
both the unknown internal dynamics, -ay-by+(b-b 0 )u and 
the external disturbance w(t). 


If an estimate of f, j can be obtained, then the control law 
reduces the plant to y=(;f-j)+u 0 which is a unit-gain double 
integrator control problem with a disturbance (f-jp). 

Thus, rewrite the plant in (8) in state space form as 

' *i =x 2 (9) 

*2 = *3 + b 0 u 

*3 = h 

,y = x i 

with x 3 =;f added as an augmented state, and h=f is seen as an 
unknown disturbance. Now j can be estimated using a state 
observer based on the state space model 

x = Ax + Bu + Eh (10) 


y = Cz 

where 



0 

1 

0 


0 


0 

A = 

0 

0 

1 

, B = 

bo 

,C=[1 0 0], E = 

0 


.0 

0 

0 . 


. 0 . 


.1. 


Now the state space observer, denoted as the linear extended 
state observer (LESO), of (10) can be constructed as 

z =Az+Bu+L(y-y) 
y=Cz 

and if f is known or partially known, it can be used in the 
observer by taking h=;f to improve estimation accuracy. 
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z=Az+Bu+L(y-y)+Eh specific. However, all of these and other previous methods 

apply to systems where the model is known. 

y=Cz (lia) Model inaccuracy can also create tracking problems. The 


The observer can be reconstructed in software, for example, 
and L is the observer gain vector, which can be obtained using 
various methods known in the art like pole placement, 

( 12 ) 

where [ ] T denotes transpose. With the given state observer, 
the control law can be given as: 

~Z3 + «o (13) 

w = — - — 

bo 

Ignoring the inaccuracy of the observer, 

?=(f-z 3 )+ W 0~ W 0 (14) 

which is an unit gain double integrator that can be imple- 
mented with a PD controller 

u 0 =k p {r-z l )-k < f 2 (15) 

Tracking Control 

Command following refers to the output of a controlled 
system meeting design requirements when a specified refer- 
ence trajectory is applied. Oftentimes, it refers to how closely 
the output y compares to the reference input r at any given 
point in time. This measurement is known as the error e=r-y. 

Control problems can be categorized in two major groups; 
point-to-point control and tracking control. Point-to-point 
applications call for a smooth step response with minimal 
overshoot and zero steady state error, such as when control- 
ling linear motion from one position to the next and then 
stopping. Since the importance is placed on destination accu- 
racy and not on the trajectory between points, conventional 
design methods produce a controller with inherent phase lag 
in order to produce a smooth output. Tracking applications 
require precise tracking of a reference input by keeping the 
error as small as possible, such as when controlling a process 
that does not stop. Since the importance is placed on accu- 
rately following a changing reference trajectory between 
points, the problem here is that any phase lag produces unac- 
ceptably large errors in the transient response, which lasts for 
the duration of the process. Although it does not produce a 
response without overshoot, it does produce a much smaller 
error signal than the point-to-point controller. The signifi- 
cance is in its ability to reduce the error by orders of magni- 
tude. A step input may be used in point-to-point applications, 
but a motion profile should be used in tracking applications. 

Various methods have been used to remove phase lag from 
conventional control systems. All of them essentially modify 
the control law to create a desired closed loop transfer func- 
tion equal to one. As a result, the output tracks the reference 
input without any phase lag and the effective bandwidth of the 
overall system is improved. The most common method is 
model inversion where the inverse of the desired closed loop 
transfer function is added as a prefilter. Another method pro- 
posed a zero Phase Error Tracking Controller (ZPETC) that 
cancels poles and stable zeros of the closed loop system and 
compensates for phase error introduced by un-cancelable 
zeros. Although it is referred to as a tracking controller, it is 
really a prefilter that reduces to the inverse of the desired 
closed loop transfer function when unstable zeros are not 
present. Other methods consist of a single tracking control 
law with feed forward terms in place of the conventional 
feedback controller and prefilter, but they are application 


performance of model-based controllers is largely dependent 
5 on the accuracy of the model. When linear time-invariant 
(LTI) models are used to characterize nonlinear time- varying 
(NTV) systems, the information becomes inaccurate over 
time. As a result, gain scheduling and adaptive techniques are 
developed to deal with nonlinearity and time variance, 
to respectively. However, the complexity added to the design 
process leads to an impractical solution for industry because 
of the time and level of expertise involved in constructing 
accurate mathematical models and designing, tuning, and 
maintaining each control system. 

15 There have been a number of high performance tracking 
algorithms that consist of three primary components: distur- 
bance rejection, feedback control, and phase error compen- 
sation implemented as a prefilter. First, disturbance rejection 
techniques are applied to eliminate model inaccuracy with an 
20 inner feedback loop. Next, a stabilizing controller is con- 
structed based on a nominal model and implemented in an 
outer feedback loop. Finally, the inverse of the desired closed 
loop transfer function is added as a prefilter to eliminate phase 
lag. Many studies have concentrated on unifying the distur- 
25 bance rejection and control part, but not on combining the 
control and phase error compensation part, such as the RIC 
framework. Internal model control (IMC) cancels an equiva- 
lent output disturbance. B. Francis and W. Wonham, “The 
Internal Model Principal of Control Theory ” Automatica, vol 
30 1 2, 1976, pp. 457-465. E. Schrijver and J. van Dijk, “Distur- 
bance Observers for Rigid Mechanical Systems Equivalence, 
Stability, and Design,” Journal of Dynamic Systems, Mea- 
surement, and Control , vol. 124, December 2002, pp. 539- 
548 uses a basic tracking controller with a DOB to control a 
35 multivariable robot. The ZPETC has been widely used in 
combination with the DOB framework and model based con- 
trollers. 

Thus, having reviewed controllers and observers, the appli- 
cation now describes example systems and methods related to 
40 controllers and observers. 

Web Processing Applications 

Web tension regulation is a challenging industrial control 
problem. Many types of material, such as paper, plastic film, 
cloth fabrics, and even strip steel are manufactured or pro- 
45 cessed in a web form. The quality of the end product is often 
greatly affected by the web tension, making it a crucial vari- 
able for feedback control design, together with the velocities 
at the various stages in the manufacturing process. The ever- 
increasing demands on the quality and efficiency in industry 
50 motivate researchers and engineers alike to explore better 
methods for tension and velocity control. However, the highly 
nonlinear nature of the web handling process and changes in 
operating conditions (temperature, humidity, machine wear, 
and variations in raw materials) make the control problem 
55 challenging. 

Accumulators in web processing lines are important ele- 
ments in web handling machines as they are primarily respon- 
sible for continuous operation of web processing lines. For 
this reason, the study and control of accumulator dynamics is 
60 an important concern that involves a particular class of prob- 
lems. The characteristics of an accumulator and its operation 
as well as the dynamic behavior and control of the accumu- 
lator carriage, web spans, and tension are known in the art. 

Both open-loop and closed-loop methods are commonly 
65 used in web processing industries for tension control pur- 
poses. In the open-loop control case, the tension in a web span 
is controlled indirectly by regulating the velocities of the 
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rollers at either end of the web span. An inherent drawback of 
this method is its dependency on an accurate mathematical 
model between the velocities and tension, which is highly 
nonlinear and highly sensitive to velocity variations. Never- 
theless, simplicity of the controller outweighs this drawback 
in many applications. Closing the tension loop with tension 
feedback is an obvious solution to improve accuracy and to 
reduce sensitivity to modeling errors. It requires tension mea- 
surement, for example, through a load cell, but is typically 
justified by the resulting improvements in tension regulation. 

Most control systems will unavoidably encounter distur- 
bances, both internal and external, and such disturbances 
have been the obstacles to the development of high perfor- 
mance controller. This is particularly true for tension control 
applications and, therefore, a good tension regulation scheme 
must be able to deal with unknown disturbances. In particular, 
tension dynamics are highly nonlinear and sensitive to veloc- 
ity variations. Further, process control variables are highly 
dependent on the operating conditions and web material char- 
acteristics. Thus, what are needed are systems and methods 
for control that are not only overly dependent on the accuracy 
of the plant model, but also suitable for the rejection of sig- 
nificant internal and external disturbances. 

Jet Engine Control Applications 

A great deal of research has been conducted towards the 
application of modern multivariable control techniques on 
aircraft engines. The majority of this research has been to 
control the engine at a single operating point. Among these 
methods are a multivariable integrator windup protection 
scheme (Watts, S. R. and S. Garg (1996). “An Optimized 
Integrator Windup Protection Technique Applied to a Turbo - 
fan Engine Control ” AIAA Guidance Navigation and Control 
Conj. ’), a tracking filter and a control mode selection for 
model based control (Adibhatla S. and Z. Gastineau (1994). 
“Tracking Filter Selection .And Control Mode Selection For 
Model Based Control T AIAA 30th Joint Propulsion Confer- 
ence and Exhibit ), an H m method and linear quadratic Gaus- 
sian with loop transfer recovery method (Watts, S. R. and S. 
Garg (1995). “A Comparison Of Multivariable Control 
Design Techniques For A Turbofan Engine Control.” Inter- 
national Gas Turbine and Aeroengine Congress and Expo .), 
and a performance seeking control method (Adibhatla, S. and 
K. L. Johnson (1993). “Evaluation of a Nonlinear Psc Algo- 
rithm on a Variable Cycle Engine.” AIAA/SAE/ASME/ASEE 
29 th Joint Propulsion Conference and Exhibit ). Various 
schemes have been developed to reduce gain scheduling 
(Garg, S. (1997). “A Simplified Scheme for Scheduling Mul- 
tivariable Controllers .” IEEE Control Systems) and have even 
been combined with integrator windup protection and H m 
(Frederick, D. K., S. Garg and S. Adibhatla (2000). “Turbofan 
Engine Control Design Using Robust Multivariable Control 
Technologies. IEEE Trans, on Control Systems Technology). 

Conventionally, there have been a limited number of con- 
trol techniques for full flight operation (Garg, S. (1997). “A 
Simplified Scheme for Scheduling Multivariable Control- 
lers.” IEEE Control Systems; and Polley, J. A., S. Adibhatla 
and P. J. Hoffman (1988). “Multivariable Turbofan Engine 
Control for Full Conference on Decision and Control Flight 
Operation.” Gas Turbine and Expo). However, there has been 
no development of tuning a controller for satisfactory perfor- 
mance when applied to an engine. Generally, at any given 
operating point, models can become inaccurate from one 
engine to another. This accuracy increases with model com- 
plexity, and subsequently design and timing complexity. As a 
result, very few of these or similar aircraft design studies have 
led to implementation on an operational vehicle. 
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The current method for controlling high performance jet 
engines remains multivariable proportional -integral (PI) con- 
trol (Edmunds, J. M. (1979). “Control System Design Using 
Closed-Loop Nyquist and Bode Arrays.” Int. J on Control, 
5 30:5, 773-802, and Polley, J. A, S. Adibhatla and P. J. Hoff- 
man (1 988). “Multivariable Turbofan Engine Control for Full 
Conference on Decision and Control. Flight Operation.” Gas 
Turbine and Expo). Although the controller is designed by 
implementing Bode and Nyquist techniques and is tunable, a 
10 problem remains due to the sheer number of tuning param- 
eters compounded by scheduling. 

Health Monitoring and Fault Detection 

The terms “health”, “fault”, “diagnosis”, and “tolerance” 
15 are used in broad terms. Some literature defines a fault as an 
unpermitted deviation of at least one characteristic property 
or variable by L. H. Chiang, E. Russell, and R. D. Braatz, 
Fault Detection and Diagnosis in Industrial Systems , 
Springer- Verlag, February 2001. Others define it more gen- 
20 erally as the indication that something is going wrong with the 
system by J. J. Gertler, “Survey of model -based failure detec- 
tion and isolation in complex plants,” IEEE Control Systems 
Magazine, December 1988. 

Industry is increasingly interested in actively diagnosing 
25 faults in complex systems. The importance of fault diagnosis 
can be seen by the amount of literature associated with it. 
There are a number of good survey papers by (J. J. Gertler, 
“Survey of model-based failure detection and isolation in 
complex plants,” IEEE Control Systems Magazine , Decem- 
30 ber 1988. V. Venkatasubramanian, R. Rengaswamy, K. Yin, 
and S. N. Kavuri, “A review of process fault detection and 
diagnosis part i: Quantitative model-based methods,” Com- 
puters and Chemical Engineering, vol. 27, pp. 293-3 1 1, April 
2003, (P. M. Frank, “Fault diagnosis in dynamic systems 
35 using analytical and knowledge-based redundancy: a survey 
and some new results ,” Automatica, vol. 26, no. 3, pp. 459- 
474, 1990, K. Madani, “A survey of artificial neural networks 
based fault detection and fault diagnosis techniques,” Inter- 
national Joint Conference on Neural Networks , vol. 5, pp. 
40 3442-3446, July 1999, P. M. Frank, “Analytical and qualita- 
tive model-based fault diagnosis-a survey and some new 
results,” European Journal of Control, 1996, P. M. Frank and 
X. Ding, “Survey of robust residual generation and evaluation 
methods in observer-based fault detection,” Journal of Pro- 
45 cess Control, 1997, J. Riedesel, “A survey of fault diagnosis 
technology [for space power systems],” in Proceedings of the 
24 th Intersociety IECEC- 89. Conversion Engineering Con- 
ference, 1989, pp. 183-188, A. Willsky, “A survey of design 
methods for failure detection in dynamic systems,” NASA 
50 STI/Recon Technical Reports, vol. 76, pp. 11 347-+, 1975, 
M. Kinnaert, “Fault diagnosis based on analytical models for 
linear and nonlinear systems — a tutorial,” Department of 
Control Engineering and System Analysis, Universite Libre 
de Bruxelles, Tech. Rep., 2004.) and books by (L. H. Chiang, 
55 E. Russell, and R. D. Braatz, Fault Detection and Diagnosis 
in Industrial Systems , Springer- Verlag, February 2001, M. 
Blanke, M. Kinnaert, J. Junze, M. Staroswiecki, J. Schroder, 
and J. Lunze, Diagnosis and Fault-Tolerant Control , 
Springer- Verlag, August 2003, R. Patton, P. M. Frank, and R. 
60 N. Clark, Issues of Fault Diagnosis for Dynamic Systems , 
Springer- Verlag Telos, 2000, S. Simani, C. Fantuzzi, and R. 
Patton, Model-based Fault Diagnosis in Dynamic Systems 
Using Identification Techniques. Springer- Verlag, January 
2003, E. Russell, L. H. Chiang, and R. D. Braatz, Data- 
65 Driven Methods for Fault Detection and Diagnosis in Chem i- 
cal Processes ( Advances in Industrial Control). Springer- 
Verlag, 2000, M. Basseville and I. V. Nikiforov, Detection of 



US 8,060,340 B2 


11 

Abrupt Changes: Theory and Application. Prentice-Hall, Inc, 
April 1993.) which collect many of the issues and solutions 
for faults. 

There are four main categories of fault diagnosis. Fault 
detection is the indication that something is going wrong with 
the system. Fault isolation determines the location of the 
failure. Failure identification is the determination of the size 
of the failure. Fault accommodation and remediation is the act 
or process of correcting a fault. Most fault solutions deal with 
the first three categories and do not make adjustments to 
closed loop systems. The common solutions can be catego- 
rized into a six major areas: 

1. Analytical redundancy by (J. J. Gertler, “Survey of 
model-based failure detection and isolation in complex 
plants,” IEEE Control Systems Magazine, December 
1988, A. Will sky, “A survey of design methods for fail- 
ure detection in dynamic systems,” NASA STI/Recon 
Technical Report N , vol. 76, pp. 1 1 347-+, 1975, E. Y. 
Chow and A. S. Willsky, “Analytical redundancy and the 
design of robust failure detection systems,” IEEE Trans- 
actions on Automatic Control, October 1982.) 

2. Statistical analysis by (L. H. Chiang, E. Russell, and R. 
D. Braatz, Fault Detection and Diagnosis in industrial 
Systems , Springer- Verlag, February 2001, E. Russell, L. 
H. Chiang, and R. D. Braatz, Data-Driven Methods for 
Fault Detection and Diagnosis in Chemical Processes 
{Advances in Industrial Control). Springer- Verlag, 
2000, M. Basseville and I. V. Nikiforov, Detection of 
Abrupt Changes: Theory and Application. Prentice- 
Hall, Inc, April 1993.) 

3. Knowledge/fuzzy logic systems 

4. Neural networks by (K. Madani, “A survey of artificial 
neural networks based fault detection and fault diagno- 
sis techniques,” International Joint Conference on Neu- 
ral Networks , vol. 5, pp. 3442-3446, July 1999, J. W. 
Hines, D. W. Miller, and B. K. Hajek, “Fault detection 
and isolation: A hybrid approach,” in American Nuclear 
Society Annual Meeting and Embedded Topical Meeting 
on Computer-Based Human Support Systems Technol- 
ogy, Methods and Future , Philadelphia, Pa., Oct. 
29-Nov. 2 1995.) 

5. Hybrid solutions by (J. W. Hines, D. W. Miller, and B. K. 
Hajek, “Fault detection and isolation: A hybrid 
approach,” in American Nuclear Society Annual Meet- 
ing and Embedded Topical Meeting on Computer-Based 
Human Support Systems: Technology, Methods and 
Future , Philadelphia, Pa., Oct. 29-Nov. 2 1995.) 

6. Fault tolerant control by (M. Kinnaert, “Fault diagnosis 
based on analytical models for linear and nonlinear sys- 
tems — a tutorial,” Department of Control Engineering 
and System Analysis, Universite Libre de Bruxelles, 
Tech. Rep., 2004, M. Blanke, M. Kinnaert, J. Junze, M. 
Staroswiecki, J. Schroder, and J. Lunze, Diagnosis and 
Fault-Tolerant Control , Springer- Verlag, August 2003) 

Some of these methods attempt to remove the need for 
accurate mathematical models yet require other implicit mod- 
els. Analytical redundancy by (E. Y. Chow and A. S. Willsky, 
“Analytical redundancy and the design of robust failure 
detection systems,” IEEE Transactions on Automatic Con- 
trol , October 1982.), the most popular method, relies heavily 
on mathematical models. 

Often detailed model information is not available although 
diagnostics of the dynamic control system are still important. 
A less developed but important problem is characterizing 
what can be determined from input output data with few 
assumptions about the plant. 
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Without adequate knowledge of the plant, disturbances, 
faults, and modeling errors, it is difficult to build an effective 
estimator. For the most part, each of these issues has been 
approached independently. 

5 

SUMMARY 

This section presents a simplified summary of methods, 
systems, and computer readable media and so on for scaling 
to and parameterizing controllers to facilitate providing a basic 
understanding of these items. This summary is not an exten- 
sive overview and is not intended to identify key or critical 
elements of the methods, systems, computer readable media, 
and so on or to delineate the scope of these items. This 
15 summary provides a conceptual introduction in a simplified 
form as a prelude to the more detailed description that is 
presented later. 

The application describes scaling and parameterizing con- 
trollers. With these two techniques, controller designing, tun- 
20 ing, and optimizing can be improved. In one example, sys- 
tems, methods, and so on described herein facilitate reusing a 
controller design by scaling a controller from one application 
to another. This scaling may be available, for example, for 
applications whose plant differences can be detailed through 
25 frequency scale and/or gain scale. While PID controllers are 
used as examples, it is to be appreciated that other controllers 
can benefit from scaling and parameterization as described 
herein. 

Those familiar with filter design understand that filters may 
30 be designed and then scaled for use in analogous applications . 
Filter designers are versed in the concept of the unit filter 
which facilitates scaling filters. In example controller scaling 
techniques, a plant transfer function is first reduced to a unit 
gain and unit bandwidth (UGUB) form. Then, a known con- 
35 troller for an appropriate UGUB plant is scaled for an analo- 
gous plant. Since certain plants share certain characteristics, 
classes of UGUB plants can be designed for which corre- 
sponding classes of scaleable, parameterizable controllers 
can be designed. 

40 Since certain classes of plants have similar properties, it is 

possible to frequency scale controllers within classes. For 
example, an anti-lock brake plant for a passenger car that 
weighs 2000 pounds may share a number of characteristics 
with an anti-lock brake plant for a passenger car that weighs 
45 2500 pounds. Thus, if a UGUB plant can be designed for this 
class of cars, then a frequency scaleable controller can also be 
designed for the class of plants. Then, once a controller has 
been selected and engineered for a member of the class (e.g., 
the 2000 pound car), it becomes a known controller from 
50 which other analogous controllers can be designed for other 
similar cars (e.g., the 2500 pound car) using frequency scal- 
ing. 

This scaling method makes a controller “portable”. That is 
a single controller can be used as the “seed” to generate 
55 controllers for a large number of different plants that are 
similar in nature. The remaining question concerns how to 
account for differences in design requirements. Controller 
parameterization addresses this issue. The example param- 
eterization techniques described herein make controller coef- 
60 ficients functions of a single design parameter, namely the 
crossover frequency (also known as the bandwidth). In doing 
so, the controller can be tuned for different design require- 
ments, which is primarily reflected in the bandwidth require- 
ment. 

65 The combination of scaling and parameterization methods 
means that an existing controller (including PID, TFB, and 
SFSOB) can be scaled for different plants and then, through 
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the adjustment of one parameter, changed to meet different 
performance requirements that are unique in different appli- 
cations. 

Certain illustrative example methods, systems, computer 
readable media and so on are described herein in connection 
with the following description and the annexed drawings. 
These examples are indicative, however, of but a few of the 
various ways in which the principles of the methods, systems, 
computer readable media and so on may be employed and 
thus are intended to be inclusive of equivalents. Other advan- 
tages and novel features may become apparent from the fol- 
lowing detailed description when considered in conjunction 
with the drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Prior Art FIG. 1 illustrates the configuration of an output 
feedback control system. 

FIG. 2 illustrates a feedback control configuration. 

FIG. 3 illustrates an example controller production system. 

FIG. 4 illustrates an example controller scaling method. 

FIG. 5 illustrates an example controller scaling method. 

FIG. 6 compares controller responses. 

FIG. 7 illustrates loop shaping. 

FIG. 8 illustrates a closed loop simulator setup. 

FIG. 9 compares step responses. 

FIG. 10 illustrates transient profile effects. 

FIG. 11 compares PD and LADRC controllers. 

FIG. 12 illustrates LESO performance. 

FIG. 13 is a flowchart of an example design method. 

FIG. 14 is a schematic block diagram of an example com- 
puting environment. 

FIG. 15 illustrates a data packet. 

FIG. 16 illustrates sub-fields within a data packet. 

FIG. 17 illustrates an API. 

FIG. 18 illustrates an example observer based system. 

FIG. 19 is block diagram of a web processing system, that 
includes a carriage and a plurality of web spans, in accor- 
dance with an exemplary embodiment; 

FIG. 20 is a linear active disturbance rejection control 
based velocity control system, in accordance with an exem- 
plary embodiment; 

FIG. 21 is an observer based tension control system, in 
accordance with an exemplary embodiment; 

FIG. 22 illustrates a desired exit speed in association with 
a carriage speed on a web processing line, in accordance with 
an exemplary embodiment; 

FIG. 23 shows a simulated disturbance introduced in a 
carriage of a web processing system, in accordance with an 
exemplary embodiment; 

FIG. 24 shows a simulated disturbance introduced in a 
process and exit-side of a web processing system, in accor- 
dance with an exemplary embodiment; 

FIG. 25 shows simulated velocity and tension tracking 
errors for carriage roller by utilizing a LADRC, in accordance 
with an exemplary embodiment; 

FIG. 26 shows simulated velocity tracking errors for a 
carriage roller by IC, LBC and LADRC 1 , in accordance with 
an exemplary embodiment; 

FIG. 27 shows a simulated control signal for carriage roller 
by IC, LBC and LADRC 1, in accordance with an exemplary 
embodiment; 

FIG. 28 shows a simulated tension tracking error by LBC, 
LADRC 1, and LADRC2, in accordance with an exemplary 
embodiment; 
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FIG. 29 illustrates a methodology for design and optimi- 
zation of a cohesive LADRC, in accordance with an exem- 
plary embodiment. 

FIG. 30 is a schematic of a turbo fan in the Modular 
5 Aero -Propulsion System Simulation (MAPSS) package, in 
accordance with an exemplary embodiment. 

FIG. 31 is a component-level model of a turbofan engine 
within the MAPSS package, in accordance with an exemplary 
embodiment. 

FIG. 32 illustrates a closed loop control system that 
employs an observer. 

FIG. 33 illustrates an ADRC for a first order system, in 
accordance with an exemplary embodiment. 

FIG. 34 illustrates an ADRC for a second order system, in 
accordance with an exemplary embodiment. 

15 FIG. 35 illustrates a single-input single-output unity gain 
closed loop system, in accordance with an exemplary 
embodiment. 

FIG. 36 illustrates a multiple single-input single-output 
loop system, in accordance with an exemplary embodiment. 
20 FIG. 37 is a graph that shows the ADRC controller’s 
responses comparing controlled variables at Operating Point 
#1 for various levels of engine degradation, in accordance 
with an exemplary embodiment. 

FIG. 38 is a graph that shows the Nominal controller’s 
25 responses comparing controlled variables at Operating Point 
#1 for various levels of engine degradation, in accordance 
with an exemplary embodiment. 

FIG. 39 illustrates a disturbance rejection model. 

FIG. 40 illustrates a current discrete estimator system, in 
accordance with an exemplary embodiment. 

FIG. 41 illustrates an open-loop tracking error plot, in 
accordance with an exemplary embodiment. 

FIG. 42 illustrates a model of a canonical form system with 
disturbance, in accordance with an exemplary embodiment. 

FIG. 43 illustrates a plot of a response of an industrial 
35 motion control test bed to a square torque disturbance, in 
accordance with an exemplary embodiment. 

FIG. 44 illustrates a plot of a response of an industrial 
motion control test bed to a triangular torque disturbance, in 
accordance with an exemplary embodiment. 

40 FIG. 45 illustrates a plot of a response of an industrial 
motion control test bed to a sinusoidal torque disturbance, in 
accordance with an exemplary embodiment. 

FIG. 46 illustrates a block a diagram of a second order 
active disturbance rejection control (ADRC) system with 
45 phase compensation, in accordance with an exemplary 
embodiment. 

FIG. 47 illustrates a block a diagram of a second order 
ADRC system with tracking, in accordance with an exem- 
plary embodiment. 

FIG. 48 is a plot of tracking of a transient profile, in accor- 
50 dance with an exemplary embodiment. 

FIG. 49 is a diagram of a dynamic estimation system for 
fault and health monitoring, in accordance with an exemplary 
embodiment. 

FIG. 50 illustrates the input-output characteristics for sys- 
55 tern diagnostics, in accordance with an exemplary embodi- 
ment. 

FIG. 51 shows a structure of disturbances, health degrada- 
tion and faults, in accordance with an exemplary embodi- 
ment. 

60 FIG. 52 is a control design broken into the estimation law, 
rejection law and nominal control law, in accordance with an 
exemplary embodiment. 

LEXICON 

65 

As used in this application, the term “computer compo- 
nent” refers to a computer-related entity, either hardware, 
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firmware, software, a combination thereof, or software in 
execution. For example, a computer component can be, but is 
not limited to being, a process running on a processor, a 
processor, an object, an executable, a thread of execution, a 
program and a computer. By way of illustration, both an 
application running on a server and the server can be com- 
puter components. One or more computer components can 
reside within a process and/or thread of execution and a 
computer component can be localized on one computer and/ 
or distributed between two or more computers. 

“Computer communications”, as used herein, refers to a 
communication between two or more computers and can be, 
for example, a network transfer, a file transfer, an applet 
transfer, an email, a hypertext transfer protocol (HTTP) mes- 
sage, a datagram, an object transfer, a binary large object 
(BLOB) transfer, and so on. A computer communication can 
occur across, for example, a wireless system (e.g., IEEE 
802.1 1), an Ethernet system (e.g., IEEE 802.3), a token ring 
system (e.g., IEEE 802.5), a local area network (LAN), a wide 
area network (WAN), a point-to-point system, a circuit 
switching system, a packet switching system, and so on. 

“Logic”, as used herein, includes but is not limited to 
hardware, firmware, software and/or combinations of each to 
perform a function(s) or an action(s). For example, based on 
a desired application or needs, logic may include a software 
controlled microprocessor, discrete logic such as an applica- 
tion specific integrated circuit (ASIC), or other programmed 
logic device. Logic may also be fully embodied as software. 

An “operable connection” is one in which signals and/or 
actual communication flow and/or logical communication 
flow may be sent and/or received. Usually, an operable con- 
nection includes a physical interface, an electrical interface, 
and/or a data interface, but it is to be noted that an operable 
connection may consist of differing combinations of these or 
other types of connections sufficient to allow operable con- 
trol. 

“Signal”, as used herein, includes but is not limited to one 
or more electrical or optical signals, analog or digital, one or 
more computer instructions, a bit or bit stream, or the like. 

“Software”, as used herein, includes but is not limited to, 
one or more computer readable and/or executable instruc- 
tions that cause a computer or other electronic device to 
perform functions, actions and/or behave in a desired manner. 
The instructions may be embodied in various forms like rou- 
tines, algorithms, modules, methods, threads, and/or pro- 
grams. Software may also be implemented in a variety of 
executable and/or loadable forms including, but not limited 
to, a stand-alone program, a function call (local and/or 
remote), a servelet, an applet, instructions stored in a memory, 
part of an operating system or browser, and the like. It is to be 
appreciated that the computer readable and/or executable 
instructions can be located in one computer component and / 
or distributed between two or more communicating, co-op- 
erating, and/or parallel processing computer components and 
thus can be loaded and/or executed in serial, parallel, mas- 
sively parallel and other manners. It will be appreciated by 
one of ordinary skill in the art that the form of software may 
be dependent on, for example, requirements of a desired 
application, the environment in which it runs, and/or the 
desires of a designer/programmer or the like. 

“Data store”, as used herein, refers to a physical and/or 
logical entity that can store data. A data store may be, for 
example, a database, a table, a file, a list, a queue, a heap, and 
so on. A data store may reside in one logical and/or physical 
entity and/or may be distributed between two or more logical 
and/or physical entities. 
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To the extent that the term “includes” is employed in the 
detailed description or the claims, it is intended to be inclusive 
in a manner similar to the term “comprising” as that term is 
interpreted when employed as a transitional word in a claim. 
5 To the extent that the term “or” is employed in the claims 
(e.g., A or B) it is intended to mean “A or B or both”. When the 
author intends to indicate “only A or B but not both”, then the 
author will employ the term “A or B but not both”. Thus, use 
of the term “or” in the claims is the inclusive, and not the 
to exclusive, use. See Bryan A. Gamer, A Dictionary of Modem 
Legal Usage 624 (2d Ed. 1995). 

DETAILED DESCRIPTION 

15 Example methods, systems, computer media, and so on are 
now described with reference to the drawings, where like 
reference numerals are used to refer to like elements through- 
out. In the following description for purposes of explanation, 
numerous specific details are set forth in order to facilitate 
20 thoroughly understanding the methods, systems, computer 
readable media, and so on. It may be evident, however, that 
the methods, systems and so on can be practiced without these 
specific details. In other instances, well-known structures and 
devices are shown in block diagram form in order to simplify 
25 description. 

Scaling 

Controllers typically are not scalable and thus are not por- 
table between applications. However, controllers can be made 
portable via scaling as described in the example systems and 
30 methods provided herein. In general, a plant mathematically 
represented by a transfer function G (s), (where s is the 
Laplace Transform variable) can be scaled according to: 

Tj p {s)=kG p {s/u> p ) (16) 

35 where co^ is the plant frequency scale and k is the gain scale, 
to represent a large number of plants that differ from the 
original plant by a frequency scale, oo^, and a gain scale, k. 

Then, a corresponding controller G c (s) for the plant G^(s) 
can be scaled according to: 

40 

G c (s)=mG c (s/«> p ). (17) 

Consider a unit feedback control system 200 with the plant 
G^(s) 210 and controller G c (s) 220, as shown in FIG. 2. 
Assume that G c (s) 220 was designed with desired command 
45 following, disturbance rejection, noise rejection, and stability 
robustness. Now, consider a similar class of plants kG^s/oo^). 
For given oo^, using example systems and methods described 
herein, a suitable controller can be produced through fre- 
quency scaling. Thus define co^ as the frequency scale and k as 
50 the gain scale of the plant G p (s/(o p ) with respect to G^(s). Then 

G c (s)=(l/k)G c (s/a p ). (18) 

Referring to FIG. 3, an example system 300 that employs 
frequency scaling is illustrated. The system 300 includes a 
55 controller identifier 310 that can identify a known controller 
associated with controlling a known plant. The controller 
may have one or more scaleable parameters (e.g., frequency, 
gains) that facilitate scaling the controller. The controller 
identifier 310 may access a controller information data store 
60 3 30 and/or a plant information data store 340 to facilitate 
characterizing one or more properties of the known control- 
ler. By way of illustration, the controller identifier 310 may 
identify the frequency scale of the controller (oo c ) and/or the 
frequency scale (oo^) and transfer function (s) of a plant con- 
65 trolled by the known controller. 

The controller information data store 330 may store, for 
example, controller class information and/or information 
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concerning scaleable controller parameters. Similarly, the 
plant data store 340 may store, for example, plant information 
like transfer function shape, frequency scale, and so on. 

The system 300 may also include a controller scaler 320 
that produces a scaled controller from the identified scaleable 5 
parameter. The scaler 320 may make scaling decisions based, 
for example, on information in the controller information data 
store 330 (e.g., controller class, scaleable parameters, fre- 
quency scale), information in the plant information data store 
3 4 0 (e .g . plant class , plant transfer function, frequency scale), 1 0 
and so on. 

While illustrated as two separate entities, it is to be appre- 
ciated that the identifier 310 and scaler 320 could be imple- 
mented in a single computer component and/or as two or 
more distributed, communicating, co-operating computer 15 
components. Thus, the entities illustrated in FIG. 3 may com- 
municate through computer communications using signals, 
carrier waves, data packets, and so on. Similarly, while illus- 
trated as two separate data stores, the controller information 
data store 330 and the plant information data store 340 may be 20 
implemented as a single data store and/or distributed between 
two or more communicating, co-operating data stores. 

Aspects of controller scaling can be related to filter design. 

In filter design, with the bandwidth, the pass band, and stop 
band requirements given, filter design is straight forward. An 25 
example filter design method includes finding a unit band- 
width filter, such as an nth order Chebeshev filter H(s), that 
meets the pass band and stop band specifications and then 
frequency scaling the filter as H(s/oo 0 ) to achieve a bandwidth 
of oo 0 . 30 

Revisiting the system 200 in FIG. 2, to facilitate under- 
standing frequency scaling and time scaling as related to 
controllers, denote oo p as the frequency scale of the plant 
Gp^/oo^) with respect to G (s) 210, and x p =ll(o p , the corre- 
sponding time scale. Then denote k as the gain scale of the 35 
plant kGp(s) with respect to Gp(s) 21 0 . With these definitions 
in hand, differences in example industrial control problems 
can be described in terms of the frequency and gain scales. 

For example, temperature processes with different time con- 
stants (in first order transfer functions), motion control prob- 40 
lems with different inertias, motor sizes, frictions, and the like 
can be described in terms of the defined frequency and gain 
scales. 

These scales facilitate paying less attention to differences 
between controllers and applications and more attention to a 45 
generic solution for a class of problems because using the 
scales facilitates reducing linear time invariant plants, proper 
and without a finite zero, to one of the following example 
forms: 


111 11 1 ( 19 ) 

5 + 1 ’ 5 ’ 5 2 + 2%S + 1 ’ 5(5 + 1) ’ 5 2 ’ 5 3 + %iS 2 + %2 S + 1 ’ 

through gain and frequency scaling. For example, the motion 
control plant of Gp(s)=23.2/s(s+1.41) is a variation of a 
generic motion control plant Gp(s)=l/s(s+l) with a gain fac- 
tor ofk=11.67 and 0^=1 .41. 


order transfer function response. Additionally, equation (19) 
can be appended by terms like: 


5 + 1 s 2 +2£ z s+1 

s 2 + 2£s + r s 3 +^s 2 +&s+r ■" 


( 21 ) 


to include systems with finite zeros. Thus, while a set of 
examples is provided in equations (19) and (21), it is to be 
appreciated that a greater and/or lesser number of forms can 
be employed in accordance with the systems and methods 
described herein. Furthermore, in some examples, scaling can 
be applied to reflect the unique characteristics of certain prob- 
lems. For example, a motion control system with significant 
resonant problems can be modeled and scaled as 



scaling 11 


1 

5(5 + 1 ) 



(22) 


where the resonant frequencies satisfy o:) rp =mo p , co rz =mco^. 
Problems with multiple frequency scales, co , noo^, and moo^, 
can be referred to as multi-scale problems. With these defi- 
nitions in hand, an example controller scaling technique is 
now described. 

Assume G c (s) is a stabilizing controller for plant G^(s), and 
the loop gain crossover frequency is oo c (s), then the controller 

Tj c (s)=G c (s/to p )lk (23) 

will stabilize the plant G /7 (s)=kG /7l (s/co p ). The new controller 
new loop gain 

L(s)=G p (s)G c (s) (24) 

will have a bandwidth of <jo c <jo^, and substantially the same 
stability margins of 

L(s)=G p (s)G c (s) 

since 


L(s)=L(s/(ti p ). 

Note that the new closed-loop system has substantially the 
same frequency response shape as the original system except 
that it is shifted by 00 ^. Thus, feedback control properties like 
bandwidth, disturbance and noise rejection are retained, as is 
the stability robustness, from the previous design, except that 
frequency ranges are shifted by 00 ^. 

Now that controller scaling has been described, PID scal- 
ing can be addressed. According to the frequency scale prin- 
ciple discussed above, and assuming the original controller 
for G^s) is a PID, e.g., 


23.2 11.67 


5(5+1.41) s ( s 
1.41 U.41 


( 20 ) 


Ic ■ 

G c (5) = k p + — + k^s 


(25) 


65 

Equation (19) describes many example industrial control 
problems that can be approximated by a first order or a second 


then the new controller for the plant kG p (s/co p ) is obtained 
from (25) as 
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That is, the new PID gains, k p , Iq, an are obtained from the 
original ones as 


k p - 



^iCOp k d 

k ’ d 


(27) 1Q 


To demonstrate the practical application and tangible 
results possible from the method described above, in the 
following example, consider a plant that has a transfer func- 
tion of 


G p (s) = 


1 

S 2 + 5 + 1 


20 


and the PID control gains of kp=3, lq=l, and k d =2. Now, 
assume the plant has changed to 

25 


G P (S)~ 


U0- 


+ Io +1 


30 


The new gains are calculated from equation (30) as kp=3, 
k.=10, k y =0.2. Thus, rather than having to build, design, and 
tune the controller for the plant 


35 


G p (s) = - 


lloJ + To + 1 


from scratch, the PID designer was able to select an existing 
PID appropriate for the PID class and scale the PID. Thus, 
frequency scaling facilitates new systems and methods for 
controller design that take advantage of previously designed 
controllers and the relationships between controllers in 45 
related applications. 

In one example, the controller is a PID controller. The PID 
controller may have a plant frequency scale 00 ^ as a scaleable 
parameter. In another example, the method includes produc- 
ing the scaled controller. F or example, a computer component 50 
may be programmed to perform the frequency scaled control- 
ling. Additionally, computer executable portions of the 
method may be stored on a computer readable medium and/or 
be transmitted between computer components by, for 
example, carrier waves encoding computer executable 55 
instructions. 

In view of the exemplary systems shown and described 
below, example methodologies that are implemented will be 
better appreciated with reference to the flow diagrams of 
FIGS. 4, 5 and 13 . While for purposes of simplicity of expla- 60 
nation, the illustrated methodologies are shown and described 
as a series of blocks, it is to be appreciated that the method- 
ologies are not limited by the order of the blocks, as some 
blocks can occur in different orders and/or concurrently with 
other blocks from that shown and described. Moreover, less 65 
than all the illustrated blocks may be required to implement 
an example methodology. Furthermore, additional and/or 
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alternative methodologies can employ additional, not illus- 
trated blocks. In one example, methodologies are imple- 
mented as computer executable instructions and/or opera- 
tions, stored on computer readable media including, but not 
limited to an application specific integrated circuit (ASIC), a 
compact disc (CD), a digital versatile disk (DVD), a random 
access memory (RAM), a read only memory (ROM), a pro- 
grammable read only memory (PROM), an electronically 
erasable programmable read only memory (EEPROM), a 
disk, a carrier wave, and a memory stick. 

In the flow diagrams, rectangular blocks denote “process- 
ing blocks” that may be implemented, for example, in soft- 
ware. Similarly, the diamond shaped blocks denote “decision 
blocks” or “flow control blocks” that may also be imple- 
mented, for example, in software. Alternatively, and/or addi- 
tionally, the processing and decision blocks can be imple- 
mented in functionally equivalent circuits like a digital signal 
processor (DSP), an ASIC, and the like. 

A flow diagram does not depict syntax for any particular 
programming language, methodology, or style (e.g., proce- 
dural, object-oriented). Rather, a flow diagram illustrates 
functional information one skilled in the art may employ to 
program software, design circuits, and so on. It is to be appre- 
ciated that in some examples, program elements like tempo- 
rary variables, initialization of loops and variables, routine 
loops, and so on are not shown. 

Turning to FIG. 5, a flowchart for an example method 500 
for producing a controller is illustrated. The method 500 
includes, at 510, identifying a controller G c (s) that stabilizes 
a plant G^(s) where the controller has a frequency 00 c and, at 
520, producing a controller G c (s) by scaling the controller 
G c (s) according to G c (s)=G c (s/ oo^j/k, where the controller 
(j c (s) will stabilize the plant G (s)HcG ^( s/cg ), where 00 ^ is 
the frequency scale of the plant G /3 (s/co p ), and where k is the 
gain scale of the plant kG^(s). In one example, the controller 
is a PID controller of the form 


G c (s) 


ki 

= JL + - 


+ k d s , 


where k^ is a proportional gain, lq is an integral gain, and k d is 
a derivative gain. In another example, 


G c (s) = 




In yet another example, the PID gains k p , lq, and k d are 
obtained from the k^, k z and k^ according to 


k p - 



kiojp _ k d 

~T' d ~ k^Tp 


It is to be appreciated that this example method can be 
employed with linear and/or non-linear PIDs. 

Applying a unit step function as the set point, the responses 
of an original controller and a scaled controller are shown in 
FIG. 6, demonstrating that the response of the scaled control- 
ler is substantially the same as the response of the original 
controller, but scaled by t=1/ol) 0 . The gain margins of both 
systems are substantially infinite and the phase margins are 
both approximately 82.372 degrees. The 0 dB crossover fre- 
quency for both systems are 2.3935 and 23.935 r/s, respec- 
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TABLE I-continued 


21 

tively. Thus, the PID scaled by the example method is demon- 
strably appropriate for the application. 

While the method described above concerned linear PIDs, 
it is to be appreciated that the method can also be applied to 
scaling nonlinear PIDs. For example, PID performance can 5 
be improved by using nonlinear gains in place of the linear 
ones. For example, 

u=k p g p {e)+k^g i {e)dt+k d g d {e) (28) 

where g p (e), g z (e), and g^e) are nonlinear functions. The to 
non-linear PIDs can be denoted NPID. Nonlinearities are 
selected so that the proportional control is more sensitive to 
small errors, the integral control is limited to the small error 
region — which leads to significant reduction in the associate 
phase lag — and the differential control is limited to a large 15 
error region, which reduces its sensitivity to the poor signal to 
noise ratio when the response reaches steady state and the 
error is small. 

The NPID retains the simplicity of PID and the intuitive 
tuning. The same gain scaling formula (30) will also apply to 20 
the NPID controller when the plant changes from G^(s) to 
kGpCs/tOp). 

Scaling facilitates concentrating on normalized control 
problems like those defined in (22). This facilitates selecting 
an appropriate controller for an individual problem by using 25 
the scaling formula in (26) and the related systems and meth- 
ods that produce tangible, results (e.g., scaled controller). 
This further facilitates focusing on the fundamentals of con- 
trol, like basic assumptions, requirements, and limitations. 
Thus, the example systems, methods, and so on described 30 
herein concerning scaling and parameterization can be 
employed to facilitate optimizing individual solutions given 
the physical constraints of a problem. 

Parameterization 

Working with controllers can be simplified if they can be 35 
described in terms of a smaller set of parameters than is 
conventionally possible. Typically, a controller (and possibly 
an observer) may have many (e.g. 15) parameters. The sys- 
tems and methods described herein concerning parameteriza- 
tion facilitate describing a controller in terms of a single 40 
parameter. In one example, controller parameterization con- 
cerns making controller parameters functions of a single vari- 
able, the controller bandwidth o) c . 

Considering the normalized plants in (19) and assuming 
desired closed-loop transfer functions are: 45 


EXAMPLES OF (0 ( 

.-PARAMETERIZED CONTROLLERS 

Gc(s, (O c ) 

6l) c (s + 1) 


, s 2 4- 2£s 4- 1 

£i>p(s + 1) <j? c s 


S 


or. 

s(s 4- 2co c ) 

s -l- 2xo c s + 2io c 


Loop shaping design can also be parameterized. Loop- 
shaping refers to manipulating the loop gain frequency 
response, L(jtt))=G i? (jcL))G £ ,(jcL)), as a control design tool. One 
example loop-shaping method includes converting design 
specifications to loop gain constraints, as shown in FIG. 7 and 
finding a controller G c (joo) to meet the specifications. 

As an example of loop shaping, considering the plants of 
the form G^s), in Table I, the desired loop gain can be 
characterized as 


( S 4- CO | \ m 

Us) = G„(j)G c (,s) = (— 


1 1 




(30) 


where oo c is the bandwidth, and 

a>i<(O c> , (0 2 <(0„ mi^O, and ni^O (3 1) 

are selected to meet constrains shown in FIG. 7. In the 
example, both m and n are integers. In one example, default 
values for otq and oo 2 are 

(O^wyiO and (02=10(0,, (32) 

which yield a phase margin greater than forty-five degrees. 

Once appropriate loop gain constraints are derived and the 
corresponding lowest order L(s) in (33) is selected, the con- 
troller can be determined from 



1 1 


s 

— + 1 



cr p l (s) 


(33) 


An additional constraint on n is that 


1 1 




G~ p l (s) 


(34) 


oj it ( 29 ) 

5 + (O c ’ (s + (O c ) 2 ’ (5 + (O c ) 3 ’ 

then for second order plants, the damping ratio can be set to 
unity, resulting in two repeated poles at -oo c . The same tech- 
nique can also be applied to higher order plants. 

Applying pole-placement design to the first and second 
order plants in (22), a set of example to c parameterized con- 
trollers are obtained and shown in Table I. Information con- 
cerning the plants and the related controllers can be stored, for 
example, in a data store. 

TABLE I 


EXAMPLES OF (0 C - PARAMETERIZED CONTROLLERS 

Gp(s) 11 1 11 

s + 1 s s 2 + 2£s +1 s(s + 1) s 2 


is proper 

This design is valid for plants with a minimum phase. For a 
non-minimum phase plant, a minimum phase approximation 
of G /? _1 (s) can be employed. 

A compromise between oq and the phase margin can be 
made by adjusting co 1 upwards, which will improve the low 
frequency properties at the cost of reducing phase margin. A 
similar compromise can be made between phase margin and 
co 2 . 

Turning to FIG. 4, an example method 400 for scaling a 
controller is illustrated. The method 400 includes, at 410, 
identifying a known controller in a controller class where the 
known controller controls a first plant. The method 400 also 
includes, at 420, identifying a scaleable parameter for the 
known controller. At 430, the method 400 includes identify- 
ing a desired controller in the controller class, where the 
desired controller controls a second, frequency related plant 
and at 440, establishing the frequency relation between the 
known controller and the desired controller. At 450, the 


50 


55 


60 


65 
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method 400 scales the known controller to the desired con- 
troller by scaling the scaleable parameter based, at least in 
part, on the relation between the known controller and the 
desired controller. 

Practical Optimization Based on a Hybrid Scaling and 5 
Parameterization Method 

Practical controller optimization concerns obtaining opti- 
mal performance out of existing hardware and software given 
physical constraints. Practical controller optimization is mea- 
sured by performance measurements including, but not lim- 
ited to, command following quickness (a.k.a. settling time), 
accuracy (transient and steady state errors), and disturbance 
rejection ability (e.g., attenuation magnitude and frequency 
range). Example physical constraints include, but are not 15 
limited to, sampling and loop update rate, sensor noise, plant 
dynamic uncertainties, saturation limit, and actuation signal 
smoothness requirements. 

Conventional tuning relies, for example, on minimizing a 
cost function like H 2 and H 28 . However, conventional cost 20 
functions may not comprehensively reflect the realities of 
control engineering, and may, therefore, lead to suboptimal 
tuning. For example, one common cost function is math- 
ematically attractive but can lead to suboptimal controller 
tuning. Thus, optimizing other criteria, like oo c are considered. 25 

A typical industrial control application involves a stable 
single-input single-output (SISO) plant, where the output 
represents a measurable process variable to be regulated and 
the input represents the control actuation that has a certain 
dynamic relationship to the output. This relationship is usu- 30 
ally nonlinear and unknown, although a linear approximation 
can be obtained at an operating point via the plant response to 
a particular input excitation, like a step change. 

Evaluating performance measurements in light of physical 
limitations yields the fact that they benefit from maximum 35 
controller bandwidth co c . If poles are placed in the same 
location, then 00 c can become the single item to tune. Thus, 
practical PID optimization can be achieved with single 
parameter tuning. For example, in manufacturing, a design 
objective for an assembly line may be to make it run as fast as 40 
possible while minimizing the down time for maintenance 
and trouble shooting. Similarly, in servo design for a com- 
puter hard disk drive, a design objective may be to make the 
read/write head position follow the setpoint as fast as possible 
while maintaining extremely high accuracy. In automobile 45 
anti-lock brake control design, a design objective may be to 
have the wheel speed follow a desired speed as closely as 
possible to achieve minimum braking distance. 

In the three examples, the design goal can be translated to 
maximizing controller bandwidth co c . There are other indus- 50 
trial control examples that lead to the same conclusion. Thus, 
o.) c maximization appears to be a useful criterion for practical 
optimality. Furthermore, unlike purely mathematical optimi- 
zation techniques, oo c optimization has real world applicabil- 
ity because it is limited by physical constraints. For example, 55 
sending 00 c to infinity may be impractical because it may 
cause a resulting signal to vary unacceptably. 

As an example of how physical limitations may affect co c 
optimization, consider digital control apparatus that have a 
maximum sampling rate and a maximum loop update rate. 60 
The maximum sampling rate is a hardware limit associated 
with the Analog to Digital Converter (ADC) and the maxi- 
mum loop update rate is software limit related to central 
processing unit (CPU) speed and the control algorithm com- 
plexity. Typically, computation speeds outpace sampling 65 
rates and therefore only the sampling rate limitation is con- 
sidered. 
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As another example, measurement noise may also be con- 
sidered when examining the physical limitations of oo c opti- 
mization. For example, the oo c is limited to the frequency 
range where the accurate measurement of the process variable 
can be obtained. Outside of this range, the noise can be 
filtered using either analog or digital filters. 

Plant dynamic uncertainty may also be considered when 
examining the physical limitations of oo c optimization. Con- 
ventional control design is based on a mathematical descrip- 
tion of the plant, which may only be reliable in a low fre- 
quency range. Some physical plants exhibit erratic phase 
distortions and nonlinear behaviors at a relative high fre- 
quency range. The controller bandwidth is therefore limited 
to the low frequency range where the plant is well behaved 
and predictable. To safeguard the system from instability, the 
loop gain is reduced where the plant is uncertain. Thus, maxi- 
mizing the bandwidth safely amounts to expanding the effec- 
tive (high gain) control to the edge of frequency range where 
the behavior of the plant is well known. 

Similarly, actuator saturation and smoothness may also 
affect design. Although using transient profile helps to 
decouple bandwidth design and the transient requirement, 
limitations in the actuator like saturation, nonlinearities like 
backlash and hysteresis, limits on rate of change, smoothness 
requirements based on wear and tear considerations, and so 
on may affect the design. For example, in a motion control 
application with a significant backlash problem in the gear- 
box, excessively high bandwidth will result in a chattering 
gearbox and, very likely, premature breakdown. Thus, co c 
optimization, because it considers physical limitations like 
sampling rate, loop update rate, plant uncertainty, actuator 
saturation, and so on, may produce improved performance. 

In one controller optimization example, assume that the 
plant is minimum phase, (e.g., its poles and zeros are in the 
left half plane), that the plant transfer function is given, that 
the co c parameterized controllers are known and available in 
form of Table I, that a transient profile is defined according to 
the transient response specifications, and that a simulator 800 
of closed-loop control system as shown in FIG. 8 is available. 
It is to be appreciated that the closed loop control system 
simulator 800 can be, for example, hardware, software or a 
combination of both. In one example, the simulator incorpo- 
rates limiting factors including, but not limited to, sensor and 
quantization noises, sampling disturbances, actuator limits, 
and the like. 

With these assumptions, one example design method then 
includes, determining frequency and gain scales, 00 ^ and k 
from the given plant transfer function. The method also 
includes, based on the design specification, determining the 
type of controller required from, for example, Table I. The 
method also includes selecting the G c (s, 00 c ) corresponding to 
the scaled plant in the form of Table I. The method also 
includes scaling the controller to 



digitalizing G c (s/co /7 , oo c )/k and implementing the controller 
in the simulator. The method may also include setting an 
initial value of oo c based on the bandwidth requirement from 
the transient response and increasing oo c while performing 
tests on the simulator, until either one of the following is 
observed: 

a. Control signal becomes too noisy and/or too uneven, or 

b. Indication of instability (oscillatory behavior) 
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Consider an example motion control test bed for which the 
mathematical model of the motion system is 

y=(- 1 .41)4-23 .2T d )+23 .2 w (35) 

where y is the output position, u is the control voltage sent to 5 
the power amplifier that drives the motor, and T d is the torque 
disturbance. An example design objective for the example 
system could be rotating the load one revolution in one sec- 
ond with no overshoot. Thus, the physical characteristics of 
the example control problem are: 10 

1) lul<3.5 volt, 

2) sampling rate=l kHz, 

3) sensor noise is 0.1% white noise, 

4) torque disturbance up to 1 0% of the maximum torque, 

5 ) smooth control signal. 15 

The plant transfer function is 


GJs) = , k = 11.67 and co p 

y s ( s \ v 

"pUp +1 ) 


1.41. 


20 


Now consider the corresponding UGUB plant 


25 


Gp(s) = 


1 

s(s + 1) ' 


A PD design of 

u=k p (r-y)+k d {-y) 


with 


k p = (a c 2 and k d =2o ) c - 1 

makes the closed-loop transfer function 


35 


G c /(s) = 


0 + ^ c ) 2 ' 


40 


Considering the plant gain scale of k and the frequency scale 
of (i) p , the PD gains are then scaled as 

45 


k p = ~y = .086 oj 2 c and k d = 


2(x) c - 1 

koj p 


= .061(2w c - 1). 


50 

To avoid noise corruptions of the control signal, an approxi- 
mate differentiator 


s 


s 

10tc> c 


+1 
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is used where the corner frequency 10oo c is selected so that the 
differentiator approximation does not introduce problematic 60 
phase delays at the crossover frequency. Using a conventional 
root locus method, the one second settling time would require 
a closed-loop bandwidth of 4 rad/sec. The example single 
parameter design and tuning methods described herein facili- 
tate determining that an oo c of 20 rad/sec yields optimal per- 65 
formance under the given conditions. A comparison of the 
two designs is shown in FIG. 9. Note that a step disturbance of 
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1 volt is added at t=3 seconds to test disturbance rejection. 
Finally, a trapezoidal transient profile is used in place of the 
step command. The results are shown in FIG. 10. 
Parameterization of State Feedback and State Observer Gains 

As described in the Background section, the State Feed- 
back (SF) controller 

u=r+Kx (36) 

is based on the state space model of the plant: 

x(t)=Ax(t)+Bu(t) ,y(t)=Cx(t)+Du(t ) (3 7) 

When the state x is not accessible, a state observer (SO): 

x =Ax+Bu+L (y-y) (3 8) 

is often used to find its estimate, x. Here r is the setpoint for the 
output to follow. The state feedback gain K and the observer 
gain L are determined from the equations: 

eig (A+BK) =1 k c (s) and eig(A+LC)=K 0 (s) 

where X c (s) and X G (s) are polynomials of s that are chosen by 
the designer. Usually the K and L have many parameters and 
are hard to tune. 

The parameterization of state feedback and state observer 
gains are achieved by making 

k c (s)=(s+a c ) n and ^(sHs-kdJ" 

where (jo c and oo 0 are bandwidth of the state feedback system 
and the state observer, respectively, and n is the order of the 
system. This simplifies tuning since parameters in K and L are 
now functions of 00 c and oo 0 , respectively. 

Parameterization of Linear Active Disturbance Rejection 
Controller (LADRC) for a Second Order Plant 

Some controllers are associated with observers. Conven- 
tionally, second order systems with controllers and observers 
may have a large number (e.g., 1 5) of tunable features in each 
of the controller and observer. Thus, while a design method 
like the Han method is conceptually viable, its practical 
implementation is difficult because of tuning issues. As a 
consequence of the scaling and parameterization described 
herein, observer based systems can be constructed and tuned 
using two parameters, observer bandwidth (oo 0 ) and control- 
ler bandwidth (oo c ). 

State observers provide information on the internal states 
of plants. State observers also function as noise filters. A state 
observer design principle concerns how fast the observer 
should track the states, (e.g., what should its bandwidth be). 
The closed-loop observer, or the correction term L(y-^) in 
particular, accommodates unknown initial states, uncertain- 
ties in parameters, and disturbances. Whether an observer can 
meet the control requirements is largely dependent on how 
fast the observer can track the states and, in case of ESO, the 
disturbance f(t,xl,x2,w). Generally speaking, faster observ- 
ers are preferred. Common limiting factors in observer design 
include, but are not limited to dependency on the state space 
model of the plant, sensor noise, and fixed sampling rate. 

Dependency on the state space model can limit an appli- 
cation to situations where a model is available. It also makes 
the observer sensitive to the inaccuracies of the model and the 
plant dynamic changes. The sensor noise level is hardware 
dependent, but it is reasonable to assume it is a white noise 
with the peak value 0.1% to 1% of the output. The observer 
bandwidth can be selected so that there is no significant 
oscillation in its states due to noises. A state observer is a 
closed-loop system by itself and the sampling rate has similar 
effects on the state observer performance as it does on feed- 
back control. Thus, an example model independent state 
observer system is described. 
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Observers are typically based on mathematical models. 
Example systems and methods described herein can employ 
a “model independent” observer as illustrated in FIG. 18. For 
example a plant 1820 may have a controller 1810 and an 
observer 1830. The controller 1810 may be implemented as a 5 
computer component and thus may be programmatically tun- 
able. Similarly, the observer 1830 may be implemented as a 
computer component and thus may have scaleable param- 
eters that can be scaled programmatically. Furthermore, using 
analogous scaling and parameterizing as described herein, 10 
the parameters of the observer 1830 can be reduced to oo 0 . 
Therefore, overall optimizing of the system 1800 reduces to 
tuning oo c , and oo 0 . 

Consider a simple example for controlling a second order 
plant 

y=- ay- by+ w+b (39) 

where y and u are output and input, respectively, and w is an 
input disturbance. Here both parameters, a and b, are 
unknown, although there is some knowledge of b, (e.g., b 0 «b 20 
derived from the initial acceleration of y in step response). 
Rewrite (39) as 
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and if f is known or partially known, it can be used in the 
observer by taking h=;f to improve estimation accuracy. 

z =Az+Bu+L(y-y)+Eh 

t=Cz (43 b) 

The observer can be reconstructed in software, for example, 
and L is the observer gain vector, which can be obtained using 
various methods known in the art like pole placement, 

L=[|hP 2 P 3 ] r (44) 

where [ ] T denotes transpose. With the given state observer, 
15 the control law can be given as: 

-Z3 + «o (45) 

u = — 

bo 

Ignoring the inaccuracy of the observer, 


y-ay-by +w+(b-b 0 )u+b 0 u=f+b 0 u (40) 

where j ? =-ay-by+w+(b-b 0 )u. Here f is referred to as the 
generalized disturbance, or disturbance, because it represents 
both the unknown internal dynamics, -ay-by+(b-b 0 )u and 
the external disturbance w(t). 

If an estimate of f, j can be obtained, then the control law 

-} + U 0 
u = u 

bo 

reduces the plant to y = (j J -f)+u 0 which is a unit-gain double 
integrator control problem with a disturbance Cf-J). 

Thus, rewrite the plant in (40) in state space form as 

40 

Xi =x 2 (41) 

x 2 =x 3 + b 0 u 
x 3 = h 
, y = x i 

with x 3 =f added as an augmented state, and h=f is seen as an 
unknown disturbance. Now f can be estimated using a state 
observer based on the state space model 

x=Ax + Bu + Eh (42) 

y = Cz 
where 



0 

1 

0 


0 


0 

A = 

0 

0 

1 

,B = 

bo 

,C=[1 0 0], E = 

0 


.0 

0 

0 . 


. 0 . 


.1. 


Now the state space observer, denoted as the linear extended 
state observer (LESO), of (42) can be constructed as 

z=Az+Bu+L(y-y ) 
y=Cz 


y=(f-z 3 )+U(f*Uo (46) 

25 

which is an unit gain double integrator that can be imple- 
mented with a PD controller 

«o=y^i)-^2 (47) 

where r is the setpoint. This results in a pure second order 
closed-loop transfer function of 


c s 2 + kdS + k p 

Thus, the gains can be selected as 

k rf =2l c and]^=0) c 2 (49) 

where oo c and £, are the desired closed loop natural frequency 
and damping ratio, t, can be chosen to avoid oscillations. Note 
that -k d z 2 , instead ofk^(r-z 2 ), is used to avoid differentiating 
45 the set point and to make the closed-loop transfer function a 
pure second order one without a zero. 

This example, illustrated in FIG. 11, shows that distur- 
bance observer based PD control achieves zero steady state 
50 error without using the integral part of a PID controller. This 
example illustrates that disturbance observer based PD con- 
trol achieves zero steady state error without using the integral 
part of a PID controller. The example also illustrates that the 
design is model independent in that the design relies on the 
55 approximate value of b in (39). The example also illustrates 
that the combined effects of the unknown disturbance and the 
internal dynamics are treated as a generalized disturbance. By 
augmenting the observer to include an extra state, it is actively 
estimated and canceled out, thereby achieving active distur- 
60 bance rejection. This LESO based control scheme is referred 
to as linear active disturbance rejection control (LADRC) 
because the disturbance, both internal and external, repre- 
sented by is actively estimated and eliminated. 

The stability of controllers can also be examined. Let 
65 e z =x z -z z -, i=l, 2, 3. Combine equation (43) and (44) and sub- 
tract the combination from (42). Thus, the error equation can 
(43) be written: 
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for {A,B,C}. The parameters in L are functions of oo 0 . One 
e = A e e + Eh (50) example procedure for coo optimization based design is now 

where described. 

Given tolerable noise thresholds in the observer states, 


-A i o 

A e - A - LC = -A 0 1 

.-A o o. 

and E is defined in (42). The LESO is bounded input, bounded 
output (BIBO) stable if the roots of the characteristic poly- 
nomial of A e 

( 51 ) 

are in the left half plane (LHP) and h is bounded. This sepa- 
ration principle also applies to LADRC. 

The LADRC design from (43) to (46) yields a BIBO stable 
closed-loop system if the observer in (43) and (44) and the 
feedback control law (46) for the double integrator are stable, 
respectively. This is shown by combing equations (45) and 
(47) into a state feedback form of u=(l/b 0 )[-k^- k^-l]z=Fz, 
where F=(l/b 0 )[-k p -k^-l ] . Thus, the closed-loop system can 
be represented by the state-space equation of: 

x r A BF li * | [B E] 1| r I ( 52 ) 

_z\ [LC A-LC + Bf\[z\ [[5 0]\[h\ 


5 increase oo G until at least one of the thresholds is about to be 
reached or the observer states become oscillatory due to sam- 
pling delay. In general, the faster the ESO, the faster the 
disturbance is observed and cancelled by the control law. 

A relationship between co 0 and co c can be examined. One 
to example relationship is 

(0 o ~3^5(0 c (55) 

Equation (55) applies to a state feedback control system 
where oo c is determined based on transient response require- 
15 ments like the settling time specification. Using a transient 
profile instead of a step command facilitates more aggressive 
control design. In this example there are two bandwidths to 
consider, the actual control loop bandwidth co c and the equiva- 
lent bandwidth of the transient profile, o c . Part of the design 
20 procedure concerns selecting which of the two to use in (55). 
Since the observer is evaluated on how closely it tracks the 
states and oo c more indicative than oo c , on how fast the plant 
states move, co c is the better choice although it is to be appre- 
ciated that either can be employed. Furthermore, taking other 
25 design issues like the sampling delay into consideration, a 
more appropriate minimum 0 l) o is found through simulation 
and experimentation as 

(D o «5^10(0„ (56) 


where B=B/b 0 , and which is BIBO stable if its eigenvalues are 
in the LHP. By applying row and column operations, the 
closed-loop eigenvalues 


r A 

BF 

A + BF 

BF ) 

eig 


= eig 


{[ LC 

A-LC + BF 

{[ 0 

A-LC \) 


= eig(A + BF){Jeig(A-LC) 

= {roots of s z + kjs + k p ) U 
{roots of s 3 + A -s 2 + A 5 + A1 

Since r is the bounded reference signal, a nontrivial condi- 
tion on the plant is that h =j is bounded. In other words, the 
disturbance f must be differentiable. 

ESO Bandwidth Parameterization 

oo 0 parameterization refers to parameterizing the ESO on 
observer bandwidth oo 0 . Consider a plant (42) that has three 
poles at the origin. The related observer will be less sensitive 
to noises if the observer gains in (44) are small for a given co 0 . 
But observer gains are proportional to the distance for the 
plant poles to those of the observer. Thus the three observer 
poles should be placed at -oo 0 , or equivalently, 

(53) 

That is 


An example for optimizing LADRC is now presented. One 
example LADRC design and optimization method includes 
designing a parameterized LESO and feedback control law 
where co 0 and oo c are the design parameters. The method also 
includes designing a transient profile with the equivalent 
bandwidth of co c and selecting an oj 0 from (56). The method 
then includes setting oo c =oo 0 and simulating and/or testing the 
LADRC in a simulator. The method also includes incremen- 
tally increasing ol) c and oo 0 by the same amount until the noise 
levels and/or oscillations in the control signal and output 
exceed the tolerance. The method also includes incrementally 
increasing or decreasing oo c and oo 0 individually, if necessary, 
to make trade-offs between different design considerations 
like the maximum error during the transient period, the dis- 
turbance attenuation, and the magnitude and smoothness of 
the controller. 

In one example, the simulation and/or testing may not yield 
satisfactory results if the transient design specification 
described by oo c is untenable due to noise and/or sampling 
limitations. In this case, control goals can be lowered by 
reducing eo c and therefore to c and co 0 . It will be appreciated by 
one skilled in the art that this approach can be extended to 
Luenberg state observer based state feedback design. 

By way of illustration, reconsider the control problem 
example associated with equations (32), but apply the 
LADRC in (43) to (48). Note that b=23 .2 for this problem, but 
to make the design realistic, assume the designer’s estimate of 
b is b o =40. Now rewrite the plant differential equation (38) as 
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40 


45 


50 


55 


^=3(d 0 , (3 2 =3co 0 2 , |3 3 =(0 o 3 (54) 

It is to be appreciated that equations (53) and (54) are 
extendable to nth order ESO. Similarly, the parameterization 


y=(-1.41y+23.2T d )+(23.2-40)u+40u=f+40u 

60 The LESO is 


method can be extended to the Luenberger Observer for arbi- 
trary A, B, and C matrices, by obtaining {A,B,C} as observ- 
able canonical form of {A,B,C}, determining the observer 65 
gain, L, so that the poles of the observer are at -oo 0 and using 


-3o) 0 

1 

°l 

r° 

3(x) 0 

-3 (*£ 

0 

1 40 

3u? 0 1 

.-<4 

0 


1° 



the inverse state transformation to obtain the observer gain, L, 
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(60a) 


and zi=Z2~Pi{zi~y{t)) 

Zi -» y, Z2 -> j, and £2 = £3 ~ AUl - jW) 

23 -> / = -1.41> + 23.2r d + (23.2-40)«, as r -> 00 5 : 

Z n = £«+i - A (Zi - y{t )) + b 0 u. 

The control law is defined as * n+l = ~P n+l (^ ~ JW) + K 


u = ~4Q~ ^ w ° = M r ~zi)~ kdZ2 
with 

k d = 2£a) c ,f = Uandk p = u) 2 c 15 

where 00 c is the sole design parameter to be tuned. A trapezoi- 
dal transient profile is used with a settling time of one second, 
or 00^=4. From (56), co 0 is selected to be 40 rad/sec. The 20 
LADRC facilitates design where a detailed mathematical 
model is not required, where zero steady state error is 
achieved without using the integrator term in PID, where 
there is better command following during the transient stage 
and where the controller is robust. This performance is 25 
achieved by using a extended state observer. Example perfor- 
mance is illustrated in FIG. 12. 

Parameterization of LADRC for nth Order Plant 

It will be appreciated by one skilled in the art that observer 30 
based design and tuning techniques can be scaled to plants of 
arbitrary orders. For a general nth order plant with unknown 
dynamics and external disturbances, 

y (n) =S(t,y,y, . . . ,y (n ^ V) ,u,u, . . . u (n ~^ ,w)+bu (57) 

the observer can be similarly derived, starting from the state 
space equation 

*1 = x 2 (58) 

x 2 = x 3 

x n = x n+i + b 0 u 

*n+ 1 = h 

y = * 1 

with x n+1 =f added as an augmented state, and h=j mostly 
unknown. The observer of (43) in its linear form with the 
observer gain 

L=[M 2 ---!L + i] r (59) 

has the form 

zi=z 2 - A (zi - y(t)) (60) 

z .2 = Z3 - Mzi - y(t)) 

Zn = Zn+l - Pn(Zl ~ y(t)) + b 0 U 
Zn+l = -Pn+iiZi -y{t)) 

and if f is known or partially known, it can be used in the 
observer by taking h=;f to improve estimation accuracy. 


With the gains properly selected, the observer will track the 
states and yield 

z x (t)^y(t\z 2 (t)^y{t), ■ ■ ■ ^(O-TJ^Cf) 

z n+i {t)-^j{t,y,y, , / n ~ l \u,u, . . . w) (61) 

The control law can also be similarly designed as in (45) and 
(47), with 

Zn+l + «0 (62) 

« = 

bo 

which reduces the plant to approximately a unit gain cascaded 
integrator plant 

y (n) =(f-z n+1 )+u 0 ~u 0 (63) 

and 

«o=yr-Zi)-^z 2 - . . . -k dn l z n (64) 

where the gains are selected so that the closed-loop charac- 
teristic polynomial has n poles at -co c , 

s n +k dn ] s n ~ i + . . . +k dl s+k p =(s+ (x) c ) n (65) 

35 oo c is the closed-loop bandwidth to be optimized in tuning. 
The co 0 optimization can similarly be applied using 

5 M =p 1 5 M " 1 + . . . +|L_ 1 s+|L=(s+C0 o )" (66) 

The following example method can be employed to iden- 
40 tify a plant order and b 0 . Given a “black box” plant with input 
u and output y, the order, n, and b 0 can be estimated by 
allowing the plant to discharge energy stored internally so 
that it has a zero initial condition, (e.g., y(0)=y(0)= . . . 
y(«-i)(0)=0) and then assuming f(0)=0. The method includes 
45 applying a set of input signals and determining the initial 
slope of the response: y(0 + ), y(0 + ), .... The method also 
includes determining the slope y w (0 + ) that is proportional to 
u(0) under various tests, (e.g., y^\0 + )=ku(0)). Then the 
method includes setting n=i+ 1 and b 0 =k. 

50 Auto-Tuning Based on the New Scaling, Parameterization 
and Optimization Techniques 

Auto-tuning concerns a “press button function” in digital 
control equipment that automatically selects control param- 
eters. Auto-tuning is conventionally realized using an algo- 
55 rithm to calculate the PID parameters based on the step 
response characteristics like overshoot and settling time. 
Auto-tuning has application in, for example, the start up 
procedure of closed-loop control (e.g., commissioning an 
assembly line in a factory). Auto-tuning can benefit from 
60 scaling and parameterization. 

In some applications, dynamic changes in the plant during 
operations are so severe that controller parameters are varied 
from one operating point to another. Conventionally, gain- 
scheduling is employed to handle these situations. In gain- 
65 scheduling, the controller gains are predetermined for differ- 
ent operating points and switched during operations. 
Additionally, and/or alternatively, self-tuning that actively 
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adjusts control parameters based on real time data identifying 
dynamic plant changes is employed. 

Common goals of these techniques are to make the con- 
troller parameter determination automatic, given the plant 
response to a certain input excitation, say a step function and 
to maintain a consistent controller performance over a wide 
range of operations, (e.g. making the controller robust). 

Example systems, methods and so on described herein 
concerning scaling and parameterization facilitate auto-scal- 
ing model based controllers. When a transfer function model 
of a plant is available, the controller can be designed using 
either pole placement or loop shaping techniques. Thus, 
example scaling techniques described herein facilitate auto- 
mating controller design and tuning for problems including, 
but not limited to, motion control, where plants are similar, 
differing only in dc gain and the bandwidth, and adjusting 
controller parameters to maintain high control performance 
as the bandwidth and the gain of the plant change during the 
operation. 

In the examples, the plant transfer functions can be repre- 
sented as G (s)=kG (s/co ), where G^(s) is given and known 
as the “mother” plant and k and are obtained from the plant 

response or transfer function. Assuming the design criteria 
are similar in nature, differing only in terms of the loop gain 
bandwidth, oo c , the controller for similar plants can be auto- 
matically obtained by scaling the given controller, G c (s, oo c ), 
for G^(s). This is achieved by combining the controller scal- 
ing, defined in equation (26), and co c -parameterization to 
obtain the controller for G /7 (s)=kG /7 (s/co p ) as 

G c (s, to e )=G e (s/<i> p ,<i> e )/k ( 67 ) 

There are three parameters in (67) that are subject to tun- 
ing. The first two parameters, k and co^, represent plant 
changes or variations that are determined. The third param- 
eter, co c , is tuned to maximize performance of the control 
system subject to practical constraints. 

An example method for auto-tuning is now described. The 
auto -tuning method includes examining a plant G^(s) and the 
nominal controller G c (s, co c ). Given the plant G p (s) and the 
nominal controller G c (s, oo c ), the method includes performing 
off-line tests to determine k and co^ for the plant. The method 
also includes using equation (67) to determine a new control- 
ler for the plant, G^(s)=kG p (s/ co^), obtained in the previous 
act. The method also includes optimizing co c for the new plant. 

An example method for adaptive self-tuning is now 
described. The adaptive self-tuning procedure includes 
examining a plant G p (s )=kG^ (s/ oo^) , where k and a) p are sub- 
ject to change during plant operation. Given the plant 
G p (s )=kG^ (s/ co p ) , the method includes performing real time 
parameter estimation to determine k and co^ as they change. 
The method also includes determining when the performance 
of the control system is degraded beyond a pre-determined, 
configurable threshold and updating the controller using (67). 
The method also includes selectively decreasing oo c if the 
plant dynamics deviate significantly from the model kG / ,(s/ 
tOp), which causes performance and stability problems. The 
method also includes selectively increasing oo c subject to oo c - 
optimization constraints if the plant model can be updated to 
reflect the changes of the plant beyond k and (X) p . 

The LADRC technique does not require the mathematical 
model of the plant. Instead, it employs a rough estimate of the 
single parameter b in the differential equation of the plant 
(57). This estimation is denoted as b 0 and is the sole plant 
parameter in LADRC. As the dynamics of the plant changes, 
so does b. Thus, b 0 can be estimated by rewriting (57) as 

/ n H(t)+bu ( 69 ) 
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and assuming the zero initial condition, (e.g., y (z) (0)=0, i=l, 
2, . . . n-1 and f(0)=0). Then b 0 «b can be estimated by using 

b o =/ n \0 + )/u( 0 ) ( 70 ) 

5 where u(0) is the initial value of the input. It is to be appreci- 
ated that this method can be applied to both open loop and 
closed-loop configurations. For the auto-tuning purposes, the 
test can be performed off-line and a step input, u(t)=constant 
can be applied. The LADRC does not require b 0 to be highly 
io accurate because the difference, b-b 0 , is treated as one of the 
sources of the disturbance estimated by LESO and cancelled 
by control law. 

The b 0 obtained from the off-line estimation of b described 
above can be adapted for auto-tuning LADRC. An auto- 
15 tuning method includes, performing off-line tests to deter- 
mine the order of the plant and b 0 , selecting the order and the 
b 0 parameter of the LADRC using the results of the off-line 
tests, and performing a computerized auto -optimization. 

Using the controller scaling, parameterization and optimi- 
20 zation techniques presented herein, an example computer 
implemented method 1300 as shown in FIG. 13 can be 
employed to facilitate automatically designing and optimiz- 
ing the automatic controls (ADOAC) for various applica- 
tions. The applications include, but are not limited to, motion 
25 control, thermal control, pH control, aeronautics, avionics, 
astronautics, servo control, and so on. 

The method 1300, at 1310, accepts inputs including, but 
not limited to, information concerning hardware and software 
limitations like the actuator saturation limit, noise tolerance, 
30 sampling rate limit, noise levels from sensors, quantization, 
finite word length, and the like. The method also accepts input 
design requirements like settling time, overshoot, accuracy, 
disturbance attenuation, and so on. Furthermore, the method 
also accepts as input the preferred control law form like, PID 
35 form, model based controller in a transfer function form, and 
model independent LADRC form. In one example, the 
method can indicate if the control law should be provided in 
a difference equation form. At 1320, a determination is made 
concerning whether a model is available. If a model is avail - 
40 able, then at 1330 the model is accepted either in transfer 
function, differential equations, or state space form. If a 
model is not available, then the method may accept step 
response data at 1340. Information on significant dynamics 
that is not modeled, such as the resonant modes, can also be 
45 accepted. 

Once the method has received information input, the 
method can check design feasibility by evaluating the speci- 
fication against the limitations. For example, in order to see 
whether transient specifications are achievable given the limi- 
50 tations on the actuator, various transient profiles can be used 
to determine maximum values of the derivatives of the output 
base on which the maximum control signal can be estimated. 
Thus, at 1350, a determination is made concerning whether 
the design is feasible. In one example, if the design is not 
55 feasible, processing can conclude. Otherwise, processing can 
proceed to 1360. 

If the input information passes the feasibility test, then at 
1360, the method 1300 can determine an oo c parameterized 
solution in one or more formats. In one example, the oo c 
60 solution can then be simulated at 13 70 to facilitate optimizing 
the solution. 

In one example, to assist an engineer or other user, the 
ADOAC method provides parameterized solutions of differ- 
ent kind, order, and/or forms, as references. The references 
65 can then be ranked separately according to simplicity, com- 
mand following quality, disturbance rejection, and so on to 
facilitate comparison. 
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Computer Processing of Control Algorithms 

FIG. 14 illustrates a computer 1400 that includes a proces- 
sor 1402, a memory 1404, a disk 1406, input/output ports 
1410, and a network interface 1412 operably connected by a 
bus 1408. Executable components of the systems described 
herein may be located on a computer like computer 1400. 
Similarly, computer executable methods described herein 
may be performed on a computer like computer 1400. It is to 
be appreciated that other computers may also be employed 
with the systems and methods described herein. The proces- 
sor 1402 can be a variety of various processors including dual 
microprocessor and other multi -processor architectures. The 
memory 1404 can include volatile memory and/or non-vola- 
tile memory. The non-volatile memory can include, but is not 
limited to, read only memory (ROM), programmable read 
only memory (PROM), electrically programmable read only 
memory (EPROM), electrically erasable programmable read 
only memory (EEPROM), and the like. Volatile memory can 
include, for example, random access memory (RAM), syn- 
chronous RAM (SRAM), dynamic RAM (DRAM), synchro- 
nous DRAM (SDRAM), double data rate SDRAM (DDR 
SDRAM), and direct RAM bus RAM (DRRAM). The disk 
1406 can include, but is not limited to, devices like a magnetic 
disk drive, a floppy disk drive, a tape drive, a Zip drive, a flash 
memory card, and/or a memory stick. Furthermore, the disk 
1406 can include optical drives like, compact disk ROM 
(CD-ROM), a CD recordable drive (CD-R drive), a CD 
rewriteable drive (CD-RW drive) and/or a digital versatile 
ROM drive (DVD ROM). The memory 1404 can store pro- 
cesses 1414 and/or data 1416, for example. The disk 1406 
and/or memory 1404 can store an operating system that con- 
trols and allocates resources of the computer 1400. 

The bus 1408 can be a single internal bus interconnect 
architecture and/or other bus architectures. The bus 1408 can 
be of a variety of types including, but not limited to, a memory 
bus or memory controller, a peripheral bus or external bus, 
and/or a local bus. The local bus can be of varieties including, 
but not limited to, an industrial standard architecture (ISA) 
bus, a microchannel architecture (MSA) bus, an extended 
ISA (EISA) bus, a peripheral component interconnect (PCI) 
bus, a universal serial (USB) bus, and a small computer sys- 
tems interface (SCSI) bus. 

The computer 1400 interacts with input/output devices 
1418 via input/output ports 1410. Input/output devices 1418 
can include, but are not limited to, a keyboard, a microphone, 
a pointing and selection device, cameras, video cards, dis- 
plays, and the like. The input/output ports 1410 can include 
but are not limited to, serial ports, parallel ports, and USB 
ports. 

The computer 1400 can operate in a network environment 
and thus is connected to a network 1420 by a network inter- 
face 1412. Through the network 1420, the computer 1400 
may be logically connected to a remote computer 1422. The 
network 1420 can include, but is not limited to, local area 
networks (LAN), wide area networks (WAN), and other net- 
works. The network interface 1412 can connect to local area 
network technologies including, but not limited to, fiber dis- 
tributed data interface (FDDI), copper distributed data inter- 
face (CDDI), Ethemet/IEEE 802.3, token ring/IEEE 802.5, 
and the like. Similarly, the network interface 1412 can con- 
nect to wide area network technologies including, but not 
limited to, point to point links, and circuit switching networks 
like integrated services digital networks (ISDN), packet 
switching networks, and digital subscriber lines (DSL). 

Referring now to FIG. 15, information can be transmitted 
between various computer components associated with con- 
troller scaling and parameterization described herein via a 
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data packet 1500. An exemplary data packet 1500 is shown. 
The data packet 1500 includes a header field 1510 that 
includes information such as the length and type of packet. A 
source identifier 1520 follows the header field 1510 and 
5 includes, for example, an address of the computer component 
from which the packet 1500 originated. Following the source 
identifier 1520, the packet 1500 includes a destination iden- 
tifier 153 0 that holds, for example, an address of the computer 
component to which the packet 1500 is ultimately destined. 
10 Source and destination identifiers can be, for example, glo- 
bally unique identifiers (guids), URLS (uniform resource 
locators), path names, and the like. The data field 1540 in the 
packet 1500 includes various information intended for the 
15 receiving computer component. The data packet 1500 ends 
with an error detecting and/or correcting field 1550 whereby 
a computer component can determine if it has properly 
received the packet 1500. While six fields are illustrated in the 
data packet 1500, it is to be appreciated that a greater and/or 
20 lesser number of fields can be present in data packets. 

FIG. 16 is a schematic illustration of sub-fields 1600 within 
the data field 1540 (FIG. 15). The sub-fields 1600 discussed 
are merely exemplary and it is to be appreciated that a greater 
and/or lesser number of sub-fields could be employed with 
25 various types of data germane to controller scaling and 
parameterization. The sub-fields 1600 include a field 1610 
that stores, for example, information concerning the fre- 
quency of a known controller and a second field 1620 that 
stores a desired frequency for a desired controller that will be 
30 scaled from the known controller. The sub -fields 1600 may 
also include a field 1630 that stores a frequency scaling data 
computed from the known frequency and the desired fre- 
quency. 

35 Referring now to FIG. 17, an application programming 
interface (API) 1700 is illustrated providing access to a sys- 
tem 1710 for controller scaling and/or parameterization. The 
API 1700 can be employed, for example, by programmers 
1720 and/or processes 1730 to gain access to processing 
40 performed by the system 1710. For example, a programmer 
1720 can write a program to access the system 1710 (e.g., to 
invoke its operation, to monitor its operation, to access its 
functionality) where writing a program is facilitated by the 
presence of the API 1700. Thus, rather than the programmer 
45 1720 having to understand the internals of the system 1710, 
the programmer’ s task is simplified by merely having to learn 
the interface to the system 1710. This facilitates encapsulat- 
ing the functionality of the system 1710 while exposing that 
functionality. Similarly, the API 1700 can be employed to 
50 provide data values to the system 1710 and/or retrieve data 
values from the system 1710. 

For example, a process 1730 that retrieves plant informa- 
tion from a data store can provide the plant information to the 
system 1710 and/or the programmers 1720 via the API 1700 
55 by, for example, using a call provided in the API 1700. Thus, 
in one example of the API 1700, a set of application program 
interfaces can be stored on a computer-readable medium. The 
interfaces can be executed by a computer component to gain 
access to a system for controller scaling and parameteriza- 
60 tion. Interfaces can include, but are not limited to, a first 
interface 1740 that facilitates communicating controller 
information associated with PID production, a second inter- 
face 1750 that facilitates communicating plant information 
associated with PID production, and a third interface 1760 
65 that facilitates communicating frequency scaling information 
generated from the plant information and the controller infor- 
mation. 
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LADRC Applied to Web Processing 

In another embodiment, a linear Active Disturbance Rejec- 
tion Control (LADRC) can be employed to provide control on 
web processing lines. LADRC requires very little informa- 
tion of the plant dynamics, has only two parameters to tune, 
and has very good disturbance rejection capability. LADRC 
controllers are inherently robust against plant variations and 
are effective in a large range of operations. 

The mathematical model of a web process line and the to 
existing control methods are illustrated. The accumulator 
dynamics are used as a test bed in association with an exem- 
plary embodiment. Generally, a web processing line layout 
includes an entry section, a process section and an exit sec- 
tion. Operations such as wash and quench on the web are 15 
performed in the process section. The entry and exit section 
are responsible for web unwinding and rewinding operations 
with the help of accumulators located in each sections. 

With reference to FIG. 19, an exemplary exit accumulator 2 o 
1900 is illustrated. Accumulators are primarily used to allow 
for rewind or unwind core changes while the process contin- 
ues at a constant velocity. Dynamics of the accumulator 
directly affect the behavior of web tension in the entire pro- 
cess line. Tension disturbance propagates along both the 25 
upstream and downstream of the accumulator due to the 
accumulator carriage. 

Since there is no difference between the entry accumulator 
and exit accumulator, except that one is for unwinding and the 
other is for rewinding operations, the embodiment discussed 
relates to exit accumulators. However, it is to be understood 
that the systems and methods described herein can relate to an 
accumulator in substantially any location within substantially 
any system (e.g., a web process line, etc.). As shown, the exit 35 
accumulator 1900 includes a carriage 1902 and web spans 
1904, 1906, 1908, 1910, 1912, 1914, and 1916. It is to be 
understood that the web spans 1904-1916 are for illustrative 
purposes only and that the number of web spans can be N, 
where N is an integer equal to or greater than one. 40 

The dynamics of the carriage tension and the entry/exit 
rollers are summarized below: 

ic(t) = — 7^(v c (f) + Lv e (t) - v p (t)) 
x c (t)y N 

X c {t) = v c (f) 

LW = ^-(-Nt c (t)-Ff(v c (t)) + u c (t))-g 
v e {t) = ji-BfeVeit) 4- R 2 (t c (t) - t r ) 4- RK e u e (t) + R 2 6 e (t )) 

Vpit) = j(-Bf Vp (t) 4- R 2 (t c (t) - t r ) + RK p u p {t) + R 2 6 p {t)) 

where v c (t), v e (t) and v p (t) are the carriage velocity, exit-side 
and process- side web velocity, respectively. x c (t) is the car- 
riage position, t r is the desired web tension in the process line 60 
and t c (t) is the average web tension. u c (t), u e (t) and u^(t) are 
the carriage, exit-side and process-side driven roller control 
inputs, respectively. The disturbance force, F^t), includes 
friction in the carriage guides, rod seals and other external 
force on the carriage. K e and K p are positive gains. 8 e (t) and 65 
8 p (t) are disturbances on the exit side and process line. The 
constant coefficients in (71) to (75) are described in Table II. 


(71) 

(72) 

(73) 

(74) 

(75) 
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TABLE II 


PLANT COEFFICIENTS 


Values 

Descriptions 

M c 

7310 kg 

Mass of the carriage 

A 

3.27 x 10“ 4 m 2 

Cross sectional area of web 

E 

6.90 x 10 10 N/m 2 

Modulus of elasticity 

R 

0.1524 m 

Radius of exit and pro cess - 



side roller 

N 

34 

Number of web spans 

J 

2.1542 kg-m 2 

Moment of inertia 

v / 

35.037 x 10 5 N-s/m 

Viscous friction coefficient 

B/ 

2.25 x 10“ 3 N-m-s 

Bearing friction coefficient 


Existing Web Tension Control Methods 

The control design objective is to determine a control law 
such that the process velocities, v c (t), v e (t) and v (t), as well 
as the tension, t c (t), all closely follow their desired trajectories 
or values. It is assumed that v c (t), v c (t) and v^(t), are measured 
and available as feedback variables. 

Typically, proportional -integral -derivative (PID) control is 
the predominant method in industry, and such control is con- 
ventionally employed with web applications. In one example, 
an industry controller can employ a feed-forward method for 
the position and velocity control of the accumulator carriage, 
and the feed-forward plus proportional -integral (PI) control 
method for the exit-side driven roller and process-side driven 
roller velocity control. The control law can be described as: 


Udit) = M c (v d c it) + g+ jj~v d c it) + ^-4) 


(76) 


M0 = ~^[~J *0 + ^ (0 “ kpeeveit) - k ie Je„(t)dT^ 

Mpli t) = ^ —j~ V p (?) 4 - V p it) — kpp^vpit) ~ ^ip ^J'^vpiO^ 


ill) 

(78) 


where U c7 (t),u e7 (t) and u^/t) are the carriage, exit-side and 
process-side driven roller control inputs. v c d , v e d and \ p d are 
the desired velocity of carriage exit-side and process-side 
rollers, respectively; and v c d , v e d and v p d their derivatives . k pe 
and k pp are proportional gains and k ie , k ip are integral gains. 

An alternative control method based on Lyapunov method 
can also be employed: 


N i AE - N \ 


U e it) -- 


RK e 


O + - t)' m - T 6 ' sgn(e 4 

u ^=lkl^ pU> + K(,> ~ 

( AE R 2 \ R 2 

] ~[Nx c it)~T) etc{t) ~T 


(BO) 


(81) 


y p^vpit) ~ 


r S p sgnie V p) 


where y 3 , y e , and y p are the controller gains to be selected. 
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The following tension observer can be used to estimate 
6(t): 


( 2 AE 

N \ 

( 2 AE R 2 j 

(82) 5 

AMO 

~w c r c{t)+ 

[ “ T J (£v ' (r) - e ' ,p W 


= t cO 




m-4. 

. «!c(0 = fc(0 - 

1 

II 

10 

= M0 - 




= v p (r) - 

- v£(r), e„(t) = 

v,(r)-v?(r). 





15 


Since the velocities are generally controlled in open-loop 
by a conventional PI feed forward and control method, the 
industrial controller needs to retune the controller when the 
operating conditions are changed and external disturbance 
appears. In addition, the industrial controller has a poor per- 
formance in the presence of disturbance. 

The Lyapunov based controller (LBC) improves the indus- 
trial controller by adding auxiliary error feedback terms to get 
better performance and disturbance rejection. However, the 2 5 
LBC has its own shortcomings since it is designed specifi- 
cally to deal with disturbances, which are introduced in the 
model. Thus, when uncertainties appear in a real-world appli- 
cation, the LBC may require re-design of the controller. 

In view of the conventional systems and methods, the 30 
exemplary embodiment was developed in the framework of 
an alternative control design paradigm, where the internal 
dynamics and external disturbances are estimated and com- 
pensated in real time. Therefore, it is inherently robust against 
plant variations and effective in disturbances and uncertain- 35 
ties in real application. In tension regulation, both open-loop 
and closed-loop options will be explored. In the open-loop 
case, the tension is not measured but indirectly controlled 
according to Equation (7 1 ) by manipulating the velocity vari- 
ables . In the closed-loop case, a tension observer is employed 40 
in the tension feedback control. 

A New Solution to Velocity and Tension Regulation 

In developing new solutions for this difficult industry prob- 
lem, performance and simplicity are stressed. That is, the new 
controller must have a much better performance than the 45 
existing ones, and it should also be simple to design, imple- 
ment, and tune. In order to provide a comprehensive control 
structure, velocity and tension are both addressed. The three 
velocity loops are very similar in nature and finding a better 
solution would be a good first step. The tension problem is 50 
crucial because of its importance and its nonlinear dynamics. 
Based on the cost and performance considerations, two solu- 
tions are discussed herein: 1) if the tension model in (1) is 
reliable, it can be well controlled with fast and accurate veloc- 
ity loops; 2) industry users are quite willing to install tension 55 
sensors for direct tension feedback control in return for better 
tension performance. FIG. 20 illustrates an exemplary veloc- 
ity control system and FIG. 21 illustrates a tension control 
system. 

FIG. 20 illustrates a LADRC-based velocity control sys- 60 
tern 2000 that employs a linear extended state observer 
(LESO) 2002. An extended state observer (ESO) is a unique 
method to solve the fault estimation, diagnosis and monitor- 
ing problem for undesired changes in dynamic systems. As an 
overview, ESO uses minimal plant information while esti- 65 
mating the rest of the unknown dynamics and unknown faults . 
This requires an observer that uses minimal plant information 
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while still being able to estimate the essential information. In 
one example, for fault problems, the important information is 
the faults and disturbances. With minimal infonnation, the 
ESO is designed to estimate these unknown dynamic varia- 
tions that compose the faults. As implemented with fault 
diagnosis, these estimated dynamics are analyzed for changes 
that represent the fault or deterioration in health. Accordingly, 
the more that is known about a relationship between the 
dynamics and a specific fault the better the fault can be iso- 
lated. The basic idea for fault remediation is that estimated 
fault information is employed to cancel the effect of the faults 
by adjusting the control to reject faults. 

The ESO system can be employed in various forms of 
dynamic systems. These include but are not limited to elec- 
trical, mechanical, and chemical dynamic systems often con- 
cerned with control problems. The most advantage would be 
achieved if this solution closes the loop of the system to 
accommodate the estimated faults. However, without 
dynamically controlling the system this method would still 
provide a benefit for health status and fault detection without 
automatically attempting to fix the fault or optimize the 
health. In one example, ESO is employed in web processing 
systems, as discussed in detail below. Other applications can 
include power management and distribution. 

The ESO offers a unique position between common meth- 
ods. There are generally two ways that the health and fault 
diagnosis problem is approached. On one side of the spectrum 
the approach is model dependent analytical redundancy. The 
other side of the spectrum is the model-less approaches from 
fuzzy logic, neural networks and statistical component analy- 
sis. The ADRC framework offers a unique position between 
these two extremes without entering into hybrid designs. The 
ESO requires minimal plant information while estimating the 
rest of the unknown dynamics and unknown faults. Further- 
more, built into the solution is a novel scheme for automatic 
closed loop fault accommodation. 

Although a single velocity loop is illustrated, it is to be 
appreciated that the control system 2000 can be applied sepa- 
rately for all three velocity loops v c (t), v e (t) and v^(t). Velocity 
regulation in a process line is one of the most common control 
problems in the manufacturing industry. Since most pro- 
cesses are well-behaved, a PID controller is generally suffi- 
cient. Other techniques, such as pole-placement and loop 
shaping, could potentially improve the performance over PID 
but require mathematical models of the process . They are also 
more difficult to tune once they are implemented. An alterna- 
tive method is described below: 

The velocity equations (73)-(75) can be rewritten as 


M0 = f c (t) + b c u c {t) 

(S3) 

M0 = M0 + b e U e {t) 

(84) 

WO = fp{t) + b p u p (t) 

(85) 

where 


m = - F f (t) - M c g ) 

(86) 

/.(f) = J(-B/V.(r) + R 2 (r, - r„(r)) + R 2 6,(t)) 

(87) 

fpW = -j(-B f V p (t) + tffcf r) - t r ) + R 2 6p(t)) 

(88) 

The plants in (84)-(86) are all of the form 


v(t)=f(t)+bu(t) 

(89) 
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where v(t) is the measure to be controlled, u(t) is the control 
signal, and the value of b is known, approximately. ;f(t) rep- 
resents the combined effects of internal dynamics and exter- 
nal disturbance. 

The key to the control design is to compensate for f(t), and 5 
such compensation is simplified if its value can be determined 
at any given time. To make such a determination, an extended 
state observer can be applied. 

Writing the plant in (89) in a state space form 

( xi = x 2 + bu (90) 

< *2 - h 
{ y = xi 

Let x 1 =v, with x 2 =f added as an augmented state, and h=f 15 
as unknown disturbance. 

The state space model is 


| x = Az + Bu + Eh 

t y = Cx 

where 



,C=[10] 


(91) 


20 


25 


Now f can be estimated using a state observer based on the 
state space model 

Based on (91), a state observer, can be constructed as 


Z = Az + Bu + L{y - 5>) 
y = Cz 


(92) 


35 


where z— >x. If f is known or partially known, it can be used in 
the observer by taking h=j to improve estimation accuracy. 

( Z = Az + Bu + L(y -y) + Eh (92a) 

t m = Cz 


The observer reduced to the following sets of state equations 45 
is the LESO. 


( ti = Z 2 + Li (y -z\) + bu (93) 

\z2=L2(y-Zl) 50 

If partial information is used, the observer is then represented 
by 

55 

j Zi =Z 2 + L(y-zi) + bu (93a) 

1*2 = L(y-zi) + h 

By setting X(s)=lsI-(A-LC)ls 2 +L 1 s+L 2 equal to the desired 
error dynamics, (s+oo) 2 , the observer gains are solved as func- 60 
tions of a single tuning parameter, co 0 . 

As known, L 1 =2 ol> 0 , L 2 =oo 0 2 can be parameterized and 
assign eigenvalues of the observer to oo 0 . With a functioning 
LESO, which results in z x ->v and z 2 ->;f , the control law will 
be designed as 65 


(94) 
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This reduces the plant to an approximate integral plant 

v{i)=m-z 2 (t))+u 0 (t)~u 0 (t) (95) 

which can be easily controlled by 

u 0 {t)=k p {r{t)-zm (96) 

For the given set point r, an approximate closed-loop trans- 
fer function is created without the addition of zeros from the 
controller. 

y(s) _ k p (97) 

r(s) s + k p 

By setting equal to the desired transfer function, co t y(s+oj j, 
the controller gains are solved as functions of one tuning 
parameter, oo c . 

Set k p =o) c , where co^ is the desired closed-loop bandwidth. 

In this example, to show how z converges to f, it is calcu- 
lated from (89) that ;f=v-bu. After solving (92), (93) and (95) 
for Z 2 by superposition, the result is a filter version of f. 

cjI (98) 

Z 2 = CsvCs) - bu{s))~ -J 

(s + CJoY 

The LESO can be further simplified by substituting (93) 
into (92) to remove an algebraic loop and decouple Z 2 , allow- 
ing ADRC to be presented in PID form 

u=k p (r-z l )-Lj(y-z l )/b (99) 

where v(t) is the measure to be controlled, u(t) is the control 
signal, and the value of b is known, approximately. f(t) rep- 
resents the combined effects of internal dynamics and exter- 
nal disturbance. 

The disturbance observer-based PD controller can achieve 
zero steady state error without using an integrator. 

The unknown external disturbance and the internal uncer- 
tain dynamics are combined and treated as a generalized 
disturbance. 

By augmenting the observer an extra state, which can be 
actively estimated and canceled out the disturbance, thereby 
achieving active disturbance rejection. 

The PD controller can be replaced with other advanced 
controller if necessary. The tuning parameters are oo 0 and oj c . 

The only parameter needed is the approximate value of b in 
(89). 

Both open-loop and closed-loop solutions to tension regu- 
lation are discussed below. The open-loop system is simple 
and economic; whereas the closed-loop system is more pre- 
cise but requires an additional sensing device. 

Open-Loop Tension Regulation 

High quality velocity regulation can allow tension in a web 
based control system to be controlled via open-loop, if the 
model of the tension dynamics (71 ) is accurate. From (7 1 ), the 
tension can be computed as 


C* AE ( 1 ) (100) 

t c (t) = t c (0) + J (^v c (r) + — (v e (t) - Vp (r))J d t 

where t c (0) is the initial value of tension. For the open-loop 
control, let the desired velocities; v c d , v/^andv^/, be carefully 
chosen so that (96) yields 

t c {t)=t c d Mh (101) 

For a given initial condition t c (0) and a given time con- 
straint, t x . Then, if all three velocity loops are well-behaved, 
the actual tension should be close to the desired value. An 


u=(-z 2 +u 0 )/ b 
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example of this method is given in simulation below. Note 
that, for this purpose, the desired velocities must satisfy the 
following condition 


vg(r) = - 


v?(r)-v£(r) 

IT 2 -’'** 


( 102 ) 
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2002, ;f c (t) J e (t) and j^(t) can be obtained. Since the other 
components f(t) are all known in this problem, tension esti- 
mation from three velocity loops can be calculated based on 
(103)-(105). 

Finally, the tension observer output value is obtained from 
the average of three tension estimations. 


The above approach is a low cost, open-loop solution. As 
the operating condition changes, the tension dynamics (1) 
could vary, causing variations in tension. If the tension is not 10 
measured, such variations may go unnoticed until visible 
effects on the product quality appear. To maintain accurate 
tension control, industry users usually are willing to install a 
tension sensor, which regulates the tension in a feedback 
loop, as discussed below. 15 

Observer Based Closed Loop Tension Regulation 

FIG. 21 illustrates an observer-based closed-loop tension 
control system 2100, wherein the system employs block dia- 
grams for the velocity and tension control loops. In this man- 
ner, the tension and velocity can be controlled at relatively the 2 q 
same time to provide real time control of a web processing 
line. A velocity controller 2102 acts as a PID controller and 
receives information from all three velocity loops, v c (t), v e (t) 
and v (t), which represent the carriage, exit-side and process- 
side driven roller velocities respectively. The velocity con- 
troller 2102 receives proportional velocity data from a veloc- 25 
ity profile bank 2104, derivative velocity data from a tension 
controller 2106, and integral velocity data from a plant 2108. 

All three inputs allow the velocity controller 2102 to maintain 
desired target values for the control signal inputs (u c (t), u e (t) 
and u^(t)) into the plant 2108 for each of the carriage, exit-side 30 
and process-side driven rollers. 

In one example, a tension meter, such as a load cell can be 
used for closed-loop tension control. Conventionally one or 
more physical instruments are required to sense the tension, 
which require additional machine space, and need adjust- 35 
ment. Therefore, implementing tension control without a ten- 
sion sensor can provide an economic benefit. Accordingly, a 
tension observer 2110 is employed to act as a surrogate for a 
hardware tension sensor to provide closed-loop tension con- 
trol. In one embodiment, the tension observer 2110 receives 40 
roller control input values (u c (t),u e (t) and u (t)) from the 
velocity controller 2102 and roller velocity values (v c (t),v e (t) 
and v^(t)) from the plant 2108. The output from the tension 
observer 2110, t c (t), is coupled with the derivative value of the 
average web tension, t c rf (t), wherein both values are input into 
the tension controller 2106. The computation of the output 45 
value of the tension observer 2110 is given below. 

Recall in (73)-(75), tension is coupled in velocity loops 
(v c (t), v e (t) and v^(t)), and an Active Disturbance Rejection 
Control (ADRC) controller can be used to decouple the ten- 
sion from the velocity loops. Actually, tension is part of the 50 
f(t) component, which is estimated and canceled out in 
LESO, as illustrated in FIG. 20. 

Considering f(t) in three velocity loops, and if the other 
parts of f(t) are known, tension can be estimated through 
equations (8 6)- (8 9) and presented as: 


M c ( 1 \ 

fecit) = j/ c (r) + -(-GW - Meg)) 

%At) = - B f v e (t) 4 - R 2 t r 4 - R 2 6 e (t)) 

tepit) = ^r( jf p (t ) 4 - BfVpit) 4 - R 2 t r + R 2 6 p (t)) 

With a proper parameter setting, the LESO 2002 can guar- 
antee that z x — >v and z 2 —>f. That is to say, from the LESO 
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(103) 

(104) 60 
(105) 

65 


l r - - \ (106) 

tc(t) = ^(tcc(t) + t„(t) + tepit)) 

Web Processing Simulation and Comparison 

In this section, four types of control systems are compared 
via simulations, including: 1) the commonly used industrial 
controller (IC) shown in equations (76) to (78); 2) the LBC in 
equations (79) to (82); 3) the three ADRC controllers, 
described in (91)-(94), for the velocity loops with tension 
regulated in open-loop (LADRC1); and 4) the same LADRC 
velocity controllers with an additional LADRC controller for 
the tension feedback loop (LADRC2). 

Note that in IC and LADRC1, the tension is controlled 
open-loop, while LADRC2 closes the tension loop with a 
tension feedback. LBC relies on the tension estimator for its 
closed-loop tension control. 

The comparison of these controllers is carried out in the 
presence of disturbances. In addition, to demonstrate the fea- 
sibility of the proposed methods, they are implemented in 
discrete-time form with a sampling period of 10 ms. 

Three control schemes are investigated by conducting 
simulations on an industrial continuous web process line. The 
desired tension in the web span is 5 1 80N. The desired process 
speed is 650 feet per minute (fpm). A typical scenario of the 
exit speed and the carriage speed during a rewind roll change 
is depicted in FIG. 22. The objective of control design is to 
make the carriage, exit velocity, and process velocities 
closely track their desired trajectories, while maintaining the 
desired average web tension level. 

To make the simulation results realistic, three sinusoidal 
disturbances are injected. F^t) in (73) is a sinusoidal distur- 
bance with the frequency of 0.5 Hz and amplitude of 44N, and 
is applied only in three short specific time intervals: 20:30 
seconds, 106:126 seconds, and 318:328 seconds as shown in 
FIG. 23. 5 e (t) and ^(t), in equation (4) and (5), are also 
sinusoidal functions with the frequency of 0.2 Hz and the 
amplitude of 44N, and is applied throughout the simulation, 
as shown in FIG. 24. 

Following the parameterization and design procedure 
described above, co c and oo 0 are the two parameters need to be 
tuned. As known in the art, the relationship between 00 c and oo 0 
is oo 0 «3|Z|5go c . So we only have one parameter to tune, which 
is w c . 

The other important parameter needed is the approximate 
value of b in (89). For this problem, the best estimate of b in 
(83), (84) and (85) is as follows: 


b c = ^~ = 1.368 x 10 -4 , b e = jK e = 0.7057, 
b p = jK p = 0.7057, b t = A *E/5 = 3.76x 10 6 

FIG. 29 illustrates a methodology 2900 for design and 
optimization of a cohesive LADRC. At 2902, a parameterized 
LESO controller is designed where co, and Co, are design 
parameters. At 2904, an approximate value of b in different 
plant is chosen. For example, b c , b e , b p , and b t , which repre- 
sents disparate known values in disparate locations within a 
web processing system. At 2906, oo 0 is set to equal 5oo c . The 
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LADRC is simulated and/or tested. In one example, a simu- 
lator or a hardware set-up is employed. At 2908, the value of 
oo c is incrementally increased until the noise levels and/or 
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the proposed LADRC controllers have a distinct advantage in 
the presence of sinusoidal disturbances and a much better 
performance in tension control. 


TABLE IV 


SIMULATION COMPARISON 

Method 


Maximum Error 


Root Mean Square Error 

v c 

(m/s) 

v e 

(m/s) 

(m/s) 

% 

(N) 

v c 

(m/s) 

(m/s) 

(m/s) 

t c 

(N) 

IC 

5.0E-4 

8.5E-3 

8.5E-3 

8.8E+4 

1.0E-4 

1.0E-3 

1.0E-3 

71.0 

LBC 1 

1.2E-4 

2.7E-3 

1.4E-3 

12.8 

3.0E-5 

5.0E-4 

6.0E-4 

11.1 

LADRC 1 

8.0E-5 

1.5E-3 

2.0E-4 

4.1 

1.0E-5 

1.0E-5 

2.0E-4 

2.8 

LADRC2 

7.0E-5 

1.3E-3 

2.0E-4 

1.5E-2 

1.0E-5 

1.0E-5 

2.0E-4 

l.E-3 


oscillations in the control signal and output exceed a desired 
tolerance. At 2910, the ratio of co^ and co 0 is modified until a 
desired behavior is observed. 

The parameters of the four controllers are shown in Table 
III. 


TABLE III 


VALUES OF THE GAINS USED IN THE SIMULATION 

Method 

Velocity Loops 

Tension Loop 

IC 

k> = 0.1k^ = 0.1 
1^ = 100,1^=100 


LBC 

y 3 = 100, y e = 100, = 100 


LADRC 1 

15,0)^ = 40, GJ^-40, 


LADRC2 

Same as ADRC1 

($ ct = 12 


Here k pe , k pp , k ie and k ip are the gains in (76)-(78) for the 
IC. y 3 , y e , and y p are the gains in (79)-(81) for the LBC. b c , b e , 
and b p are specific values of b in (92) for the carriage, exit, and 
process velocity loops, respectively. Similarly, co oc , oo oe and 
co are the observer gains in equation (91); and oo cc , w ce and 
aj cp are the controller gains (k^) in equation (94). b,, (o ct , and 
coot are the corresponding ADRC parameters for the tension 
plant in (109). 

The velocity errors (v c , v e and v ) and tension tracking 
errors t c resulting from ADRC1 are shown in FIG. 25. Obvi- 
ously, the velocity and tension tracking errors are quite small, 
despite the fact that the controller design is not based on the 
complete mathematical model of the plant and there are sig- 
nificant disturbances in the process. 

The comparisons of IC, LBC and LADRC 1 are shown in 
FIGS. 26 and 27, in terms of the tracking errors and control 
signals for the carriage velocity loop. The carriage velocity 
errors indicate that LADRC 1 is much better than the other 
two methods and the control signal indicates that the LADRC 
controller actively responds to the disturbances. It is to be 
appreciated that utilizing the systems and methods disclosed 
herein, similar characteristics can be found in the exit and 
process velocity loops. 

Due to the poor results of the IC controller, only LBC, 
LADRC 1, LADRC2 are compared in the tension control 
results in FIG. 28. With a direct tension measurement, 
LADRC2 results in negligible tension errors. Furthermore, 
even in an open-loop control, LADRC 1 has a smaller error 
than LBC. This can be attributed to the high quality velocity 
controllers in LADRC 1. 

The velocity and tension errors of all four control systems 
are summarized in Table IV. Overall, these results reveal that 


A new control strategy is proposed for web processing 
20 applications, based on the active disturbance rejection con- 
cept. It is applied to both velocity and tension regulation 
problems. Although only one section of the process, includ- 
ing the carriage, the exit, and the process stages, is included in 
this study, the proposed method applies to both the upstream 
25 and downstream sections to include the entire web line. Simu- 
lation results, based on a full nonlinear model of the plant, 
have demonstrated that the proposed control algorithm 
results in not only better velocity control but also significantly 
30 less web tension variation. The proposed method can provide 
several benefits over conventional systems and methods. For 
example, 1) no detailed mathematical model is required; 2) 
zero steady state error is achieved without using the integrator 
term in the controller; 3) improved command following is 
35 achieved during the transient stage; 4) the controller is able to 
cope with a large range of the plant’s dynamic change; and 5) 
excellent disturbance rejection is achieved. 

Additional Forms of the Extended State Observer 

Although various observers are known, such as high gain 
40 observers, sliding mode observers, and extended state observ- 
ers (ESO), it is generally regarded that the extended state 
observer is superior in dealing with dynamic uncertainties, 
disturbances, and sensor noise. Controllers that use it depend 
on quick and accurate estimation in real time of the output and 
45 equivalent disturbance as well as their derivatives. 

Observers are used to estimate variables that are internal to 
the system under control, i.e. the variables are not readily 
available outputs. Observers use a model of the system with 
correction terms and are run in continuous time. In order for 
50 continuous functions of time to run in hardware, however, 
they are often discretized and run at fixed sample rates. Dis- 
crete observers are often referred to as estimators. 

The fundamental limiting factor of the controller and esti- 
mator is the sampling rate. Improving the ESO will improve 
55 the overall performance of the system. Up to this point, Euler 
approximations have been used to implement the ESO in 
hardware, which adversely affects its performance at slower 
sampling rates. As described in greater detail herein, several 
discrete variants of extended state observers are further iden- 
60 tified and analyzed. 

There are three main contributions, discrete implementa- 
tion of the ESO or (DESO), the Generalization of the ESO and 
DESO or (GESO), and discrete parameterization of the 
DESO and GESO. 

65 Performance enhancements are made to the ESO, both in 
formulation and in implementation. Although this is referred 
to as the DESO, a number of methods are disclosed. Here, the 
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system model is first discretized using any number of meth- 
ods; Euler, zero order hold (ZOH), and first order hold (FOH). 
Then, a predictive discrete estimator (PDE) from G. F. Fran- 
klin, J. D. Powell, and M. Workman, Digital Control of 
Dynamic Systems, 3rd ed., Menlo Park, Calif: Addison Wes- 
ley Longman, Inc., 1 998, pp. 328-337 is constructed from the 
discrete model and correction terms are determined in dis- 
crete time symbolically as a function of one tuning parameter. 
It is also formulated as a Current Discrete Estimator (CDE) 
from G. F. Franklin, J. D. Powell, and M. Workman, Digital 
Control of Dynamic Systems, 3rd ed., Menlo Park, Calif.: 
Addison Wesley Longman, Inc., 1998, pp. 328-337 to main- 
tain stable operation at lower sampling rates, a major limiting 
factor in controls. Typical discretization methods, such as a 
PDE, generate at least one sample of delay, whereas a CDE 
removes this delay by adding a current time step update to the 
estimated state. Next, Euler, zero order hold (ZOH), and first 
order hold (FOH) versions of all discrete matrices are deter- 
mined symbolically to retain the simplicity of single param- 
eter tuning. In the past, only an approximation using Euler 
integration was used. The problem is that the correction terms 
were determined in continuous time and become inaccurate 
when they are increased and at low sample rates. A second 
order example is used. Simple tests show that the CDE with 
ZOH performs the best. 

The DESO is then generalized to estimate systems of arbi- 
trary order, as well as to estimate multiple extended states. 
This is referred to as the generalized ESO or (GESO). This 
reformulation incorporates a disturbance model of arbitrary 
order, thus allowing the amount of disturbance rejection to be 
specified for different types of systems. Multiple extended 
states allow the estimation of higher order derivatives of the 
disturbance, which improves the estimation of the distur- 
bance, allowing it to be more accurately cancelled. In the past, 
disturbances were restricted to first order and estimated using 
one extended state. The standard ESO does not make use of 
this information. A number of advantages exist for the current 
discrete version of the GESO. First, it offers better estimation 
and accordingly higher stability. Another implementation 
benefit is the minimal code space and processing power 
changes in addition to the standard ESO. The GESO also 
improves the performance and increases the range of opera- 
tion while maintaining a similar level of complexity. 

The immediate application of the DESO and GESO can be 
applied to ADRC controllers. Due to the current and eventual 
wide-spread application of both ADRC, the powerful GESO 
also has great immediate and future potential. Many plants or 
other control applications have physical upper limits for 
sample time and they will benefit from a stable and accurate 
estimation at lower sampling rates. They may also have a need 
to estimate higher order disturbances and they will benefit 
from higher performance control. 

The preferred embodiments described herein allow these 
advanced control methods to be a practical solution for indus- 
try to transparently implement a high performance controller 
into their systems. The problem it solves is that the usability 
of the controller will no longer suffer dramatically as a result 
adding complexity to achieve higher performance. This 
means a significant reduction in time to design, implement, 
tune, and maintain each drive in every plant and/or every 
application. 

The preferred embodiment observers have been tested in 
simulation and hardware. Results on simple test applications 
and popular motion control problems have shown stable con- 
trol at lower sampling rates than what are possible with the 
standard ESO. It was applied to ADRC with a tracking con- 
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troller. The controller was tested in a realistic simulation and 
in hardware in a motion control servo -drive 
Discrete Implementation of the Extended State Observer 
(DESO) 

5 For the sake of simplicity, consider the continuous-time 
differential equation of a second order plant where u and y are 
the input and output, respectively, and b is a constant. 

y=g(y,y,t)+w+bu (107) 

10 Combining the internal dynamics g(y,y,t) with 311 external 
disturbance w to form a generalized disturbance f(y,y,t), the 
system is rewritten as 

y=f(y,y,w,t)+bu. (108) 

An augmented state space model is constructed 

x = Ax + Bu + Ef ( 109 ) 


y = Cx + Du 



'0 

1 

O' 


"O' 


'O' 

A = 

0 

0 

1 

,B = 

b 

,E = 

0 


_0 

0 

0 . 


.0. 


1 


C = [1 0 0], D = [0] 


where x=[y,y,fl r includes the disturbance to be estimated. 
Next, an observer is created from the state space model. 

3Q x=Ax+Bu+L(y-y) 

$=C£+Du (110) 

Note that J is ignored in (110) since it is unknown and is 
estimated by the correction term. The observer is rewritten to 
35 output the state 

x=[A-LC]x+ [B-LD,L ]u c 

y c =x (111) 

40 where u c =[u, y] T is the combined input and y c is the output. It 
is then decomposed into individual state equations for the 
purpose of implementation. For the sake of simplicity, the 
observer gain vector L is determined by placing the poles of 
the characteristic equation in one location. 

45 X(5)=l5/-(^-Zqi=(5+coJ 3 

L=[3 (0 o ,3 (0 o 2 ,(0 o 3 ] T (112) 

The state space model in (109) is first discretized (formu- 
lated in discrete-time) by applying Euler, ZOH, or FOH. 

50 

x(k+ l)=<&x(k)+Tu(k) 

y{k)=Hx{k)+Ju(k ) (113) 

A discrete observer is created from this model. 

55 x{k+l)=m{k)+Tu{k)+L/y(k)-y{k)) 

y{k)=Hx{k)+Ju(k) (114) 

This is known as a predictive discrete estimator (G. F. Fran- 
6Q klin, J. D. Powell, and M. Workman, Digital Control of 
Dynamic Systems, 3rd ed., Menlo Park, Calif.: Addison Wes- 
ley Longman, Inc., 1998, pp. 328-337) because the current 
estimation error y(k)-y(k) is used to predict the next state 
estimate x(k+l). 

65 However, by defining the predictive estimator gain vector 
as 


( 115 ) 
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the estimated state reduces to 
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x(k+ l)=<S>x(k)+Tu(k) (116) 

where the new state includes a current time step update, 
giving it less time delay. 

x(k)=x(k)+LMk)-m) (117) 


50 

except for L p =TL, making the observer unstable at relatively 
low sample rates. Yet in cases where L is a nonlinear function, 
this may be the only way of discrete implementation. For the 
sake of further discussion, the past method is referred to as the 
5 Euler approximation. 

Applying ZOH 


This is referred to as a current discrete estimator G. F. Fran- 
klin, J. D. Powell, and M. Workman, Digital Control of 
Dynamic Systems, 3rd ed. , Menlo Park, Calif. : Addison Wes- 10 
ley Longman, Inc., 1998, pp. 328-337. When the sampling 
rate is low, this could play a significant role in enhancing the 
stability of a closed loop system. A block diagram is illus- 
trated in 4000 of FIG. 40. The estimator is then rewritten to 
output the new state 15 


<J> = e AT 


Z A k T k 


r = 



drB=> 


Z 


A k T k+ 1 d 


H = C,J = 0 


(122) 


y d (k)=[I-LJt]m+[~LJ,LMk) (118) 

where u/k)=[u(k), y(k )] T is the combined input and y d is the 2 o 
output. The only difference for the predictive estimator is that 
y^k)=x(k). 

Discrete Parameterization of the ESO 

For the sake of simplicity, the current estimator gain vector 
L c is determined by placing the poles of the discrete charac- 2 5 
teristic equation in one location. 

\(z)=\zI-(<S>-<t>LJf) l=(z-|3) 3 (119) 

The relation between the discrete estimator poles and the 
continuous observer poles is given as 30 

p=e“ UJor . (120) 

For example, applying Euler to (109) and solving (119) for 
L c yields 


to (109) produces a more accurate estimation than Euler. 



T 2 


t 2 


1-/F 


1 T “T 


b ^r 


3 

<J> = 

2 

0 1 T 

, r = 

2 

bT 

» L c = 

(1-/?) 2 (1+/?) — 


.0 0 1 


0 




H = [ 1 0 0], J = [0] 


Simulation and Analysis of the Discrete ESO 
Various discretization methods are analyzed through simu- 
lation of various plants. The ESO is first applied in open loop 
to a simple motion system plant model 



'I T 0 


0 ' 

<D = 

0 1 T 

, r = 

bT 


o 

o 


_ 0 _ 


l -/? 3 

(2-3 yS + y? 3 )! . 

a-P)f 2 


h = [ l o 0 ], y = [ 0 ] 


35 y=50y+500«+100w (124) 

(121) where w is a 2.5 Hz square wave starting at 0.3 sec. and u is a 
trapezoidal profile that lasts 0.125 sec. The estimator param- 
eters are ct^^OO and T=0.005. A tracking error plot is shown 
4Q in 4100 of FIG. 41 that compares the predictive and current 
discrete methods using both Euler and ZOH. The transient 
and steady state parts of each trajectory are evaluated using 
integral absolute error and then summarized in Table V. 


where T is the discrete sample time. However, note that T is 45 
denoted as the matrix transpose. In the past, the ESO was 
implemented by integrating each state equation in ( 1 1 0) using 
Euler (J. Han, “Nonlinear Design Methods for Control Sys- 
tems”, Proc. \4th IFAC World Congress, 1 999; Z. Gao, “Scal- 
ing and Bandwidth-Parameterization Based Controller Tun- 50 
ing,” American Control Conference , pp. 4989-4996, June 
2003; Z. Gao and S. Hu, “A Novel Motion Control Design 
Approach Based on Active Disturbance Rejection,” Proc. of 
the 40th IEEE Conference on Decision and Control , p. 4974, 
December 2001 ; Y. Hou, Z. Gao, F. Jiang, and B. T. Boulter, 55 
“Active Disturbance Rejection Control for Web Tension 
Regulation,” IEEE Conference on Decision and Control, 
2001; B. Sun, “Dsp-based Advanced Control Algorithms for 
a DC-DC Power Converter,” Master's Thesis, Cleveland 
State University , June 2003; R. Kotina, Z. Gao, and A. J. van 60 
den Bogert, “Modeling and Control of Human Postural 
Sway,” XXth Congress of the International Society of Biome- 
chanics, Cleveland, Ohio, July 31 -Aug. 5, 2005; R. Mikloso- 
vic and Z. Gao, “A Dynamic Decoupling Method for Con- 
trolling High Performance Turbofan Engines,” Proc. of the 65 
1 6th IFAC World Congress , Jul. 4-8, 2005. The problem with 
this method is that it produces the same matrices as (121) 


TABLE V 


OPEN LOOP TRACKING ERRORS 

Discretization 

Transient Integral Absolute Error 

Method 

y 

y’ 

F 

ZOH Current 

1.49E-6 

12E-3 

5.90 

Euler Current 

1.51E-6 

51E-3 

7.16 

ZOH Predictive 

136E-6 

41E-3 

8.41 

Euler Predictive 

134E-6 

87E-3 

9.66 

Discretization 

Steady State Integral Absolute Error 

Method 

y 

y' 

F 

ZOH Current 

0.10E-5 

9E-3 

2.83 

Euler Current 

0.10E-5 

20E-3 

2.83 

ZOH Predictive 

9.55E-5 

30E-3 

4.32 

Euler Predictive 

9.55E-5 

40E-3 

4.32 


When the step size T=0.005, the Euler approximation 
becomes unstable and therefore was not shown. However, the 
four methods shown use discrete pole placement and do not 
become unstable until T=0.066. From the table, the second 
most important option appears to be the current discrete 
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method for tracking accuracy. The table also shows that ZOH 
is better than Euler and, interestingly, dominant in estimating 
transient velocity. 

Next, the ESO is applied in closed loop to (124) and to a 
more complex simulation of an actual servo-motor. 

*4=80(75 U-0.075/J, I V m I <160, It/ 1 <8 
4=-2500( V m -0AI a - 1 .2 y) 

y=ll.l(100w+1.54) (125) 10 

With 0.4=30 and 0.4=3 00, the sample period is increased to the 
point of instability and then tabulated in Table VI. 
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In (109), previous ESO design, and in disturbances are 
considered to be piece-wise constant with h=l or a series of 
steps. Now an ESO with h=l, 2, 3 can respectively track a 
square, triangular, or parabolic disturbance. A sinusoid is a 
5 different matter because it is infinitely differentiable. How- 
ever, increasing h increases the degree of the polynomial and 
improves tracking of a sinusoid or any time varying distur- 
bance. An ESO with h extended states for a relative n th order 
plant is denoted as an ESO w h . 

The new form is represented in continuous state space 

x=Ax+Bu+Ef h) 


TABLE VI 


MAXIMUM CLOSED LOOP STEPSIZE 

Discretization 

Simple Plant 

Servo-motor 

Method 

(18) 

(19) 

Euler 

26E-4 

30E-4 

Approximation 
Euler Predictive 

37E-4 

57E-4 

Euler Current 

47E-4 

68E-4 

ZOH Predictive 

85E-4 

140E-4 

ZOH Current 

150E-4 

300E-4 


The results show that the most important option for low 
sampling time requirements is ZOH, followed by the current 
discrete method. In this regard, the current discrete ESO with 
ZOH appears to be six to ten times better than the Euler 
approximation used in previous literature. The servo system 
in (145) was also simulated, resulting in an improvement of 
5.3 times. In summary, the current discrete ESO with ZOH 
should be used for improved tracking accuracy as well as 
closed loop stability. 

Generalization of the Extended State Observer (GESO) 
Although a second order example was used in the previous 
section, (110) through (120) and (122) are applicable to a 
plant of arbitrary order with any number of extended states. 

For example, a class of general n th order plants similar to 
( 1 07) i s represented as 

y^—gty, ■ ■ • ,y (r, ~ l \t)+w+bu (126) 

where y (w) denotes the n th derivative of the output and 
g(y, . . . , y ( ” -1) , w, t) represents the internal dynamics. Two 
critical parameters are relative order n and high frequency 45 
gain b. Combining the unknowns into one generalized distur- 
bance f(y, . . . , y (w-1) , w, t) results in 

< y(”)=j ! ( y , . . . ,y ( ' n ~ i \w,t)+bu. (127) 

Note that when represented with an equivalent input distur- 50 
bance d=j7b, the design model becomes 

P d (s)=b/s n . (128) 

As a signal, the type of disturbance can be characterized 
similar that of system type in a classical control. This speci- 55 
fication is outlined in G. F. Franklin, J. D. Powell, and A. 
Emami-Naeni, Feedback Control of Dynamic Systems, 4th 
ed., Upper Saddle River, N.J.: Prentice-Hall, Inc., 2002, pp. 
239-242, 601 -604 as the degree of a polynomial that approxi- 
mates a signal, which directly relates to the number of times 60 
it is differentiated before reaching zero. Sometimes distur- 
bances are represented by a set of cascaded integrators 1/s* 
with unknown input. Under this assumption, the plant is 
represented in 4200 of FIG. 42 by two sets of cascaded 
integrators; one for the design model and another for the 65 
disturbance model. It will also be shown that this assumption 
leads to an estimated disturbance equivalent to that of a DOB. 


y=Cx+Du (129) 

where the state includes the disturbance f and its derivatives 
to be estimated. 


X = [x { , ... ,x n ,x n+l , ... , x n+h ] T (130) 

_ 5 / (0) , , f h ~ l) ] T 


Since the new form consists of cascaded integrators, the A 
matrix simplifies to an n+h square matrix with ones on the 
super diagonal. Each element of A is defined as 


1, i = j-l 

0, otherwise. 


(131) 


Since the input is added after n-integrators, the first state is 
defined as the output, and the derivative of the last state is j {h \ 
the other matrices become 

B=F„-1 b 0 A ] r , C=[ 1 0 n+h _J, E=p n+h _ i If (132) 

where 0 h represents a lxh zero vector and D=0. 

For the sake of simplicity, the observer gain vector is deter- 
mined by placing all of the poles of the characteristic equation 
in one location. 

X(5)=|5/-G-£C)I=(5+co £) )” +;i (133) 

As a result, each element in L becomes 

l i = C n+Ki e>J, *=1,2, . . . ,n+h (134) 

where the binomial coefficients are 


The ESO can also be represented in filter form 

f'^sVQyMl-QJPM i= 0, . . . ,n-l 

l 0) =d[bQ/P d - l y-u)]J=Q , . . . ,h - 1 (135) 

where binomial filters 


Q,(s) = 


fin+h,n+h-i-l (-S') 
Pn+h,n+h( s ) 


, QfW = 


fin+h,h-j- 1 6 >) 
f3n+h,n+h 0 s ) 


(136) 


consist of numerator and denominator polynomials that are 
functions of a single tuning parameter co 0 =1/t. 
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TABLE VH-continued 


J (137) 

Pij{s) = 1 + 2_J C ir(TS) r 
r= 1 

This form shows that additional extended states raise 

observer order, n+h, and increase the slope of the cutoff 
frequency. It also shows that the estimated disturbance is 
equivalent to a DOB, i.e. a filtered version of the actual f. 

S=bQj(P d - l y-u) (138) 

Discrete Implementation of the GESO 

Applying ZOH to (129) using (122) produces an n+h 
square <D matrix where each element is defined as 

yj-i, i*j (139) 

0, otherwise 


CDESO ESTIMATOR GAINS FOR ZOH AND FOH 


n + h 

K 


4 

[1 -P>, (1 -p)Hu+p(14+ 11/S))T 

10 


a-jS) 3 (i+ffi2 (i-0 4 T] T 


5 

[1 - y? 5 , (1 - p)\\+ p)(5 + pa + 5p)) A , 

15 


(l-fl 3 (7 + «10 + 7/J))-A I , 



5 1 F 

(1-Ml+«2T3,(1-«V| 


for Y*. = T*Zk! . The T matrix reduces to 

r=[bv„ ■ • • by! 0 h ] T . (140) 

If FOH is preferred, the only change is in the T and J 
matrices, which become 


A simulation of an industrial motion control test bed is used 
to demonstrate the control design procedure and its simplic- 
ity, resulting performance, and overall effectiveness in the 
absence of a simulation model. The servo amplifier, motor, 
and drive train are modeled with a resonant load as 

V m =4( L-2-054) , \ V C \ <4. 5, 1 V m I ■ <1 0 


r = 


( 2” +1 - 2 ) 
n + 1 



' bjn by i 

n + 1 2 



(141) 4=2500(^-44-0.24,), \I a Id 

r =0.54-4-4 

30 m a 1 

4=0.0005(4-44+0.0001 (x ra -44 


If Euler is preferred, the O matrix where each element is 
defined as 


( 1, i = j (142) 

<f>U = | T > * + 1 =j 

{ 0, otherwise 

and the T matrix reduces to 

r=[0 M _i bT0 h ] r . (143) 

Discrete Parameterization of the GESO 

For the sake of simplicity, the current estimator gain vector 
Lc is determined by placing the poles of the discrete charac- 
teristic equation in one location. 

X(z)=lz/-(<t>-<X>4/7) I =(z-(3)" +a (144) 

As a result, the current estimator gain vector is listed in Table 
VII as a function of n+h. A current discrete ESO with h 
extended states for a relative n th order plant is denoted as 
CDESO^. 


x„,=2500T w 

4=1754 

where V c , x /5 and T d are the control input voltage, output load 
position, and torque disturbance, respectively. Backlash of a 
±0.31 fxm/sec. dead -bandwidth onx is also applied. The con- 
trol design method using the ESO is fairly straight forward 
with only a few physical intuitions. In the most basic sense, a 
servo motor can be considered as a double integrator. 


45 */(*) ^ h _™_ (I 46 ) 

V c (s) * 52 ■ 

It is put into the new canonical form where f(t) represents any 
of the discrepancies or dynamics not modeled in (38). 

50 

x(t)=f(t)+b m V c (t) (147) 

First, a CDES0 2 A is used to estimate x z (t), x z (t), and f(t) in 
discrete time. Then the estimated disturbance is fed back to 
cancel itself 

55 


TABLE VII 

CDESO ESTIMATOR GAINS FOR ZOH AND FOH 


V c (k) = 


«<)(*) - /(*) 
b m 


(148) 


n + h 


1 

2 


3 


[i - pf 

[1 -f, (i -p)\ 1 (i -pf^ 


which reduces the system to a double integrator, x/t)«u 0 (t). 
Finally, a parameterized control law is used to control the 
augmented system where r(k) is a reference motion profile. 

u 0 (k)^2(r(k)-xM)-2MMk) (149) 

The observer and control laws in (148) and (149) are selected 
with a sample rate of 10 kHz to control the motion system’s 
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model in (145). The gain b m =25 is crudely estimated as the 
initial acceleration from a step response. Disturbance rejec- 
tion was tested by applying various torque disturbances at 
time t=l second and 0.1% white noise is injected into the 
output. Keeping the control signal within ±4.5V and its noise 
level within ±100 mV, ad and co 0 were increased to 50 and 
150, respectively. The results for a type 1 square, type 2 
triangular, and type oo sinusoidal torque disturbance are 
shown in 4300 of FIG. 43, 4400 of FIG. 44, and 4500 of FIG. 
45, respectively. Robustness was tested by increasing the load 
by a factor of nearly 8. There was no noticeable difference. 
The results show that two extended states reduce the error 
compared to one extended state. In FIG. 44, two extended 
states drive to zero the error created by type 2 disturbances. 
Although a sinusoidal disturbance is infinitely differentiable, 
three extended states significantly reduce the steady state 
error in FIG. 45. 

Various discrete implementations of the extended state 
observer are studied and compared. It is shown that the cur- 
rent discrete formulation is superior to the predictive one in 
reducing the delay associated with the sampling process. It is 
also demonstrated that the ZOH implementation improves 
estimation accuracy and stability without additional com- 
plexity to the user. To facilitate the ESO implementation for 
practitioners, the algorithm is derived symbolically with a 
single tuning parameter, i.e. the bandwidth of the observer. 
Another significant development is the generalization of the 
ESO for various types of systems and disturbances. Finally, a 
filter version shows that the estimated disturbance is equiva- 
lent to the DOB structure. Unlike the DOB, however, the ESO 
estimates suitable derivatives of the output, allowing for a 
straightforward controller design. The motion control prob- 
lem is complex with many uncertainties, yet preliminary 
results show that this observer can achieve high performance 
over a wide range of system dynamics while remaining easy 
to use. 

Tracking Control Applied to ADRC 

The various preferred embodiment controllers and observ- 
ers described herein can be used in conjunction with tracking 
components to further improve their function and perfor- 
mance. 

The immediate application of the tracking enhancement 
can be applied to ADRC controllers. Due to the current and 
eventual wide-spread application of ADRC, the powerful 
tracking enhancement also has great immediate and future 
potential. Future uses also include the specific application of 
the tracking control method to new controllers. 

The already proven ADRC control structure which works 
well for steady state set point control can now be extended to 
handle transient tracking control by use of the preferred 
embodiment tracking strategy. This enhancement allows 
these advanced control methods to be a practical solution for 
industry to transparently implement a high performance 
tracking controller into their systems. The problem it solves is 
that the usability of the controller will no longer suffer dra- 
matically as a result adding complexity to achieve higher 
tracking performance. This means a significant reduction in 
time to design, implement, tune, and maintain each drive in 
every plant of every company. 

The preferred embodiment tracking strategies have been 
applied to ADRC in the form of prefilters and/or feed forward 
tends to make the desired closed loop transfer function of 
ADRC approximately equal to one or, more generally, have a 
relative order equal to zero. Test results in simulation and 
hardware have shown error reduction up to eighty fold. 

The tracking enhancement was applied to a controller 
using ADRC and a Current Discrete Extended State Estima- 


tor (CDESO). The controller was tested in a realistic simula- 
tion and in hardware in motion system servo drive. 

Tracking ADRC Applied to a Second Order Plant 

For the sake of clarity, let us first consider a general second 
5 order plant where u and y are the input and output, respec- 
tively, and b is a constant. 

y=g(y,y,t)+w+bu (150) 

Combining the internal dynamics g(y,y,t) with an external 
io disturbance w to form a generalized disturbance f(y,y,w,t), 
the system is rewritten as 

y=f(y,y,w,t)+bu (151) 

An augmented state space model is constructed 

15 

x = Ax + Bu + Ef ( 152 ) 


y = Cx 


20 

0 

1 

O' 


'O' 


'O' 

A = 

0 

0 

1 

, B = 

b 

,E = 

0 


.0 

0 

0. 


. 0 . 


_1 . 


c =[ 1 0 0] 

25 

where x=[x l ,x 2 ,x 3 \ T =[y,y,f] T includes the disturbance. 

An ESO is then created from (152) to estimate the states 

x =Ax+Bu+L (y- Cx) (153) 

30 where x^x^^Xgj^yj}^?] 2 ^ For the sake of simplicity, the 
observer gain vector L is determined by placing the poles of 
the characteristic equation in one location. 

K(s)=\sI-(A-LC)\=(s+(a 0 ) 3 (154) 

35 9 . 7 - 

L =[3 ( 0 o ,3 ( x ) 0 2 , 0 )/] 1 

A disturbance rejection control law is applied to the plant in 
(151) to dynamically cancel f(y,y,w,t) using its estimate x 3 . 

4Q u=(u 0 -£ 3 )/b (155) 

This reduces the plant to a double integrator at low frequen- 
cies. 

y=u 0 (156) 

45 A simple control law is then applied 

Uo =k p(r-Xl)- k d*2 ( 1 57 ) 

to form the following closed loop transfer function. 


50 


Gry{s ) 5S 


k P 

s 2 + kjs + k p 


(158) 


For the sake of simplicity, it is set equal to a desired closed 
55 loop transfer function that provides a smooth step response. 


60 


= 


(S + (x ) c ) 2 


(159) 


The resulting controller gains become 

k p =(o c 2 ,k^=2to c . (160) 

65 The problem with (157) is in the phase lag it produces in 
(158). Therefore in situations requiring precise command 
following, it is proposed that the inverse of the closed loop 
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transfer function, shown in square brackets in (1 61), is added As a result, each element in L becomes 
to the reference input of the control law as a prefilter 

(« + !)! . . ( 170 ) 


to compensate for the predicted phase lag by making the new 
closed loop transfer function G^l, thus producing a much 
smaller error e=r-y than the original controller. This configu- 
ration is shown in 4600 of FIG. 46. 

A simpler way of implementation is in reducing the new 
control law in (161) to 

Uo=k p {r-x l )+k d {r-x 2 )+r (162) 

where velocity and acceleration feed forward are utilized. 
Here, the first two terms drive the error and its derivative to 
zero while the last term provides a desired control input u* 0 
such that y follows r. This equivalent tracking control con- 
figuration is shown in 4700 of FIG. 47. Even though this 
example applies tracking to a parameterized controller, note 
that the concepts in (161) and (162) will work for any linear 
time invariant controller regardless of its parameter values. 
Therefore the application of tracking to ADRC is independent 
of the application of parameterization to ADRC, whether it is 
by means of a prefilter as in (161) or by a single control law 
with feed forward terms as in (162). 

A compromise in performance between the point-to-point 
controller in (157) and the tracking controller in (162) is 
reached when using 

Uo=k p {r-^)+k d {r-^ (163) 

to produce a closed loop transfer function with a relative 
degree of one. 


A disturbance rejection control law is applied to (16) to 
dynamically cancel f(y,y,w,t) using its estimate in x w+1 . 

U=(u 0 -x n+l )/b (171) 

and reduce the plant to cascaded integrators at low frequen- 
cies. 

y ( "W 0 (172) 

A point-to-point control law is then applied 

u 0 =k 0 (r-x i )-k l x 2 - . . . (173) 

to form the following closed loop transfer function. 


For the sake of simplicity, it is set equal to a desired closed 
loop transfer function that provides a smooth step response. 


and the controller gains are determined as 


w 0 = kp 


(\s 2 +k d s + k p ] 

k r “* 1 

U K p J 


V-xi \-k d x 2 


(161) 5 ilin+l-W 0 * 0 ' 1 ’ 


10 


15 


20 


30 
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k d s + k p (164) 

s*+k d s + k p 

40 

Tracking ADRC Applied to an nth Order Plant 

When (151) is extended to arbitrary order, the system is 
represented by 

45 

y in) =f(y,y, ■ ■ ■ ,y (n ~ l) ,w,t)+bu. (165) 

Where y (M) denotes the n th with derivative of y. 

An ESO is constructed 

■ 50 

x=Ax+Bu+L(y-Cx). (166) 

The elements of the A matrix have ones on the super diagonal 


For precise tracking, the inverse of (174) is added to the 
reference input of (173) as a prefilter, making G As a 
result, the tracking control law becomes 

u 0 =k 0 (r-^ l )+ . . . +^- 1 (r (M_1) -4)+/" ) . (177) 

A single control law is formed by combining (177) with (171) 
u=k(x*-x) (178) 

where the new gains k=[k 0 , . . . , kj/b and the feed forward 
terms x*=[r, r , r (n> ]. 

Discrete Implementation of Tracking ADRC 
In hardware, a discrete ESO is created. 

x(k+ 1 ) =^>x(A:) +r u (k)+L p (y(k) -Hx(k ) ) (179) 

where x(k)=x(k)+L c (y(k)-yHx(k)) is the current time update. 
The matrices are determined by applying zero order hold. 


a i,j 


1 i=j - 1 
0, otherwise. 


and the other matrices become 


(167) 55 


Z A k T k y A k T k+l 

k=0 k=0 


H - C 


(ISO) 


B=l 0„_i b Of, C =[ 1 OJ (168) 

where 0 n represents a 1 xn zero vector. 

For the sake of simplicity, the observer gain L=[l l5 \ 2 , , 
l„ +1 fis determined by placing the poles of the characteristic 
equation in one location. 


60 The discrete estimator gain vector L p is determined by plac- 
ing the poles of the discrete characteristic equation in one 
location 

X(z)= \zl- (<J>-<J>Z,,J7) I =(z- 13)” + 1 (181) 

65 Here, the relation between the discrete estimator poles and the 
continuous observer poles is given as 


n + 1 


|3=e-“° r . 


k(s)=\sl- ( A-LC)\=(s+(x)J 


(169) 


(182) 
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For example, the matrices for the second order system 
become 


<I> = 


T 2 


T 2 


1 T — 


b— 


2 


2 


0 1 T 

. r = 

bT 

> L P - 

.001. 


0 



3-3/? 

(l-0 2 (5+0^; 


H = [l 0 0]. 


(183) 5 


10 


A discrete control law is formed 

u(k)=k(?*(k)-m) (i84) 15 

where the feed forward terms x*(k)=[r(k),r(k), . . . , r^(k)] r . 

The feed forward terms contain the reference input r and its 
derivatives. Since the reference input is generated by an algo- 
rithm and not a measured signal, it is often noise free. As a 
result, discrete differences are used for derivatives. A typical 20 
example of velocity and acceleration feed forward calcula- 
tions are given. 

i{k)={r{k)-r{k-\))IT 


: m=(m-hk-WT (185) 

Since the states of the ESO are subtracted from these sig- 
nals in the control law in (184), any difference in derivation 
between the ESO states and the feed forward terms causes 
dynamic errors. As a result, an estimator similar to the dis- 30 
crete ESO is used to estimate the feed forward terms, thus 
reducing the error and increasing performance. A model rep- 
resenting n+1 cascaded integrators and h- 1 extended states is 
used to form a discrete estimator where only the model output 
signal r is available. 35 

x*(k+l )=m* (k)+L p (y(k)-m*(k)) (186) 

where x*(k)=x*(k)+L c (y(k)-yElx*(k)) is the current time 
update. 

Another problem arises in discrete implementation of the 40 
ESO and feed forward estimators when signals within them 
grow large and create numerical errors and even instability 
due to increasing input signals. As a result, the two estimators 
are combined to form a single estimator that uses error as an 
input. Since the control law merely subtracts the states of each 45 
estimator, the matrices T>,L p , L c , H are only functions of n+h, 
and T only resides in the state feedback observer, (179) is 
subtracted from (186) to form 


Zik + 1 ) = <m) - ru(k) + Lp(e(k) - Hz(k )) 



' Z\{k) 


m-m 

z{k) = 

Zn{k) 

= 

r (n - l \k)-y (n ~ l \k) 


. Zn+l (k) _ 


r M -hk) 


50 

(187) 
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where e(k)=r(k)-y(k). This eliminates numerical errors and 
instability, cuts computation in half, and keeps all signals 60 
within the new estimator and control law small. 

The current discrete estimator form of (187) becomes 


As proof of concept, a simulation example is given. The 
setup in (1 45) is used to track a motion profile reference signal 
with a final time t =1 second. The system is simulated with 
and without a tracking controller and the results are compared 
in 4800 of FIG. 48 . Note that the tracking and reference traces 
are so close that they are on top of each other. With tracking 
control, the maximum error is reduced by a factor of 80 times 
during the initial transient shown in the first 1 .5 seconds with 
only a slight increase in maximum control effort. There is also 
no adverse affect on the disturbance rejection property of the 
controller, shown when a square wave disturbance is added 
after 1 .5 seconds. 

This discrete tracking ADRC algorithm does not require an 
explicit mathematical model to achieve high performance. It 
has one or two tuning parameters that can be adjusted quickly, 
meaning that the level of expertise, time, and resources typi- 
cally needed to construct a model, design a controller, and 
maintain performance is no longer required. This reduces 
manufacturing costs. 

Multivariable ADRC 

A general control method is given that can be applied to any 
MIMO system with the number of inputs greater than or equal 
to the number of outputs, not just jet engines. For proof of 
concept, it is then applied to dynamically decouple and con- 
trol turbofan engines, i.e. jet engines. Since the jet engine is 
one of the most complex systems in existence, it is to be 
appreciated that this will provide a good example. If one is 
able to control it without knowing its mathematical model, 
which may be several thousand lines of code, then there is 
good chance this method will work on almost any plant. For 
example, this control system can be utilized with chemical 
processes, flight control of airplanes and missiles, CNC 
machine control, robotics, magnetic bearing, satellite attitude 
control, and process control. 

Consider a system formed by a set of coupled n th order 
input-output equations 


y ( i ] =fi+biU 


(189) 


yf = f q +b q u 

where y/' !) denotes the n‘ h derivative of y,. The input U= 
K, • • u p ] r , the output y=[y 1; . . .,y q f, andb, |b ll5 . . . , b ijP ] 
for i=l, 2, . . . , q and q=p. Each equation consists of two 
terms, the instantaneous b z U and the dynamic £(Y, Y, . . . , 
yC” -1 ). All interactions between equations, internal dynam- 
ics, and external disturbances are considered part of f z (Y, 
Y, . . . , Y (w_1) , t). The system is rewritten 


Y^ u) =F+B 0 U (190) 

where Y w =[y 1 < ' 0 , . . . , y ? ( ” ) ] r , and F=[J : 1 , . . . , .f,,U and 
B 0 =[bi r , . . . , b ? r ] r - Assuming that n is known and that B is an 
qxp approximation of B 0 where both are fall row rank, a 
generalized disturbance is defined as H=F+(B 0 -B)U. The 
system reduces to 


Y (n) =H+BU. 


(191) 


f(k+l)=[®-L/k]z(k)-T u(k)+L p e(k) 
z(k)=[I-LJI]z(k) + LXk) 

where z(k) is the current time update and L^=OL c . 


The idea is to estimate H and cancel it in real time, reducing 
65 the plant to a set of cascaded integrators. In order to represent 
the plant with a set of state equations, let X=[X 1 r , X 2 r , . . . , 
X w+ /]=[Y r , Y t , . . . , Y (w " 1)r , H r , . . . , U (h ~ 1)T ] T such that 
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x l = x 2 
*„-l = x„ 

X n = X n+l +BU . 

X n + 1 — %n +2 + H 

x n+h = H {h) 


Solving for each gain as a function of m oi results in 

(192) 

, (n + h)\ j (199) 

5 * JKn + h-J)^’ 

With B + defined as the right inverse of B, a disturbance 
rejection control law is applied to (191), effectively cancel - 
ling H at low frequencies. 

U=B + (U 0 -X n+ 1) (200) 


In state space form, the plant is represented by 

X^AX+BU+Erf^ 

Y=CX (193) 

where X=[X 1 r , X 2 r , . . . , X w+ ^ r ] r , 0^ and \ q are qxq zero and 
identity matrices, and A is an q(n+h) dimensional square 
matrix. 



A 

u 

■ 

■ V 






o. 

% 


% 


^qyxj 



A = 

0 <? 

% 


■■ U 

, B = 

B 

,E = 



. % 

Qq 

0 , ■ 

_ 


Qqhxp 


. ^<7 . 


c=[l q % % ■■■ oj 


This allows a kind of feedback linearization and decoupling 
to occur which reduces the plant to a set of parallel n-integra- 
tor systems at low frequencies. 

u 0 ( 201 ) 

At this point, any number of control methods may be used. A 
simple control law with no integrators is proposed 

20 U 0 =K 0 ( Y*-Xy )-KyX 2 - . . . -K n _yX n (202) 

where Y* is the desired trajectory for Y and the controller 
gains K 0 , K l5 . . . , K w-1 are qxq matrices in general. However, 
for the sake of tuning simplicity, the controller gains are 
defined to form q parallel control loops for j=0, 1, . . . , n-1 . 

Kj=dia%(kjy, k ja , . , k M ) (203) 

Each control loop has its n poles placed in one location, for the 
sake of further simplification. 

30 


An observer is then designed from the state space model 
where L |I,,, . . . , L„ +h \ T 


A c (s) = \sl - A + = Y] (s + u c ,i) n 


(204) 


X=AX+BU+L(Y-t) 


35 


Solving for each gain as a function of oo c i results in 


(195) 


In (1 95), the state equations of the multivariable ESO become 

40 


ka = 




(205) 


k l =X 2 +L l {Y-Y) (196) 

*»-l = *„+Z*-l (>"-?) 

• X n = X n+l +L n (Y-Y) + BU- 
*»+l =X n+2 + L n+l (Y-Y) 

X nth = L n +h[Y-Y) 

The observer gains L 1; L 2 , . . . , L n+h are qxq matrices, in 
general. However, for the sake of tuning simplicity, the gains 
are defined to form q parallel observer loops for j=l, 2 , . . . , 
n+h. 

Lj=dmg(l JA ,l ja , . . . ,l M ) (197) 

Each loop then has its n+h poles placed in one location, for the 
sake of further simplification. 

* f (198) 

A 0 (s) = \sl - A+ LC\ J J (s + oj o/i ) n+h 


Typically, a nonsingular B -1 can be approximated by a 
diagonal matrix of reciprocal elements, since inaccuracies in 
B can be accounted for in H. 

45 The observer is simplified to remove B by substituting 
(200) into (196). 


X l = X 1 +L l (Y-Y) 

50 v ' 

U-i =*„+A,-i (Y-f) 

■ k n = u 0 + l„(y-y) 

55 

U+l = x„ + 2 +l„ +i (y-y) 

, Xn+h = Ln+h{Y - 0 

60 

The commonly used SISO form of ADRC is, in fact, the 
q=l case. 

Multivariable Tracking ADRC 
A tracking controller can be used in place of (202) to 
65 improve the tracking error. 


u 0 =k 0 ( Y*-Xy )+ . . . +i£ M _ 1 (y* ( ”- 1) -Y„)+r* ( ” ) 


(207) 
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Multivariable Discrete ESO 

To output the states, the ESO in (1 94) and (1 95) is rewritten 
as 

±=[A-LC]X+ [BX] [U, Yf 

Y C =X (208) 

where Y c is the state output. It can then be discretized, form- 
ing a multivariable CDESO 

Jt(k + 1 )=[W-LjI]MMTL I J[U(k)Y(k)] T 

m=[ ( 209 ) 

where L C =[L C;1 , . . . , K,„ +h ] T , U/=diag(lc,- 15 \c ja , lc^), 
and similar notation is true for L p . The simplest way to deter- 
mine the matrices is the case when (209) is equivalent to q 
parallel SISO loops. To show an example of how the matrices 
are directly extended from their SISO counterparts, the matri- 
ces for a CDES0 2 1 as a result of ZOH become 





T 2 






>7 


BT 

5> = 

0, 

T 

',T 

,r = 

B 


0 q 

o<? 

^ _ 


0 qxq 


IClj 


i-ys? 

i- c 2,i 

= 

(i-A-) 2 (i+A-)3 
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Ic^j 


(1-A) 3 i 

ft 


H — \ 0<7 0 ,] 


Turbofan Model and Design Specifications 

FIG. 30 shows an engine schematic from a turbofan engine 
3000. In this example, a Modular Aero -Propulsion System 
Simulation (MAPSS) package, developed by Parker and Guo, 
(2003) at the NASA Glenn Research Center, is employed. 
The package is used because it is comprehensive enough to 
simulate any two spool jet engine. A component-level model 
(CLM) within MAPSS consists of a two-spool, high pressure 
ratio, low bypass turbofan with mixed-flow afterburning. 

FIG. 31 illustrates a top-level control diagram 3100 of the 
turbofan engine 3000. The model consists of hundreds of 
coupled equations and look-up tables that ensure mass, 
momentum, and energy balances throughout while modeling 
gas properties effectively. Mathematical details can be found 
in (Mattingly, J. D. (1996). “Elements of Gas Turbine Pro- 
pulsion”, McGraw-Hill, Inc.; Boyce, M. P. (2002). “Gas Tur- 
bine Engineering Handbook,” Second Edition, Butterworth- 
Heinemann; Cumpsty, N. (2002). “Jet Propulsion: A Simple 
Guide,” Cambridge University Press. 

In general, the CLM is defined by two nonlinear vector 
equations 

Xcuvrt ( x clm> u clM’P ^altyX/w) 

y clm = S ( x clm> u clM’P > attorn) ( 21 1 ) 

that are functions of a 3x1 state vector (x CiM ), a 7x1 input 
vector (u Ci ^), a 10x1 health parameter vector (p), altitude 
(alt), and Mach number (xm) . A 22x 1 vector of sensor outputs 
(y clm) is combined to calculate thrust (fn), fan stall (sm2) and 
over-speed (pcn2r) margins, engine temperature ratio (etr), 
and pressure ratios of the engine (eprs), liner (lepr), a core 
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(cepr). These performance parameters form the controlled 
output. 

Y= [fn, eprs, lepr, etr, sm2,pcn2r, cepr ] T (212) 

5 Each of the seven inputs (u CLM ) is controlled by a separate 
SISO actuator consisting of a torque motor and servomecha- 
nism with saturation limits for position, velocity, and current. 
The first three actuators drive the fuel flow (wf36), variable 
nozzle exit area (a8), and rear bypass door variable area (al 6), 
io respectively. These actuator inputs form the control signal. 

U= [wf3 6 act ,a& act ,a.l 6 act ] T (213) 

The remaining four actuators drive stator and guide vane 
angles using steady state schedules within the primary control 
15 loop, ensuring safe operating limits. 

The goal of the control system is to achieve a fast thrust 
response with minimal overshoot and zero steady state error, 
while maintaining safe rotor speeds, pressure and tempera- 
ture limits, and stall margins. In MAPSS, the supplied multi- 
2 q mode controller consists of four multivariable PI regulators, 
each controlling only three outputs at one time. 

Y 1 =[fn, eprs, lepr] r 

Y 2 =[fn, etr, lepr] r 
25 

Y 3 =[fn,smn2,lepr ] T 

Y 4 =[pcn2r, cepr, lepr] r (214) 

The first regulator controls eprs at low speeds, while the 
30 second regulator controls etr at high speeds. The third and 
fourth regulators actively control limits associated with the 
fan components, namely the fan stall and over-speed margins 
when their limits are approached. Limits associated with the 
engine core are met by acceleration and deceleration sched- 
35 ules on fuel flow (Kreiner, A. and K. Lietzau (2003). “The Use 
of Onboard Real-Time Models for Jet Engine Control.” MTU 
Aero Engines , Germany). These schedules along with actua- 
tor limits are then placed to constrain the outgoing control 
signal. 

40 Multivariable ADRC Design Procedure with a Jet Engine 
Example 

A generic design procedure is given for multivariable 
ADRC, using the application to the MAPPS jet engine 3000 
in FIG. 30 as an example. Test conditions are then discussed 
45 for the jet engine, followed by simulation results that compare 
the new algorithm to the current one, showing that similar 
performance can be achieved with much less design effort. 

The design procedure for applying any of the new distur- 
bance rejection techniques is uniquely characterized by the 
50 plant representation. 

Y /i) =H+BU (215) 

In (215), the size of the input vector U and output vector Y 
should be known. The design procedure involves the deter- 
55 mination of n and B as well as a method of tuning the con- 
troller. 

A generic procedure is given, followed by a detailed expla- 
nation of each step with examples specific to the jet engine 
application. 

60 1 . Determine the number of inputs p and outputs q of the 

system. Use multimode control if q>p. 

2. Determine the high frequency gain B of the system. 

3 . Determine the relative order n of the system. If unknown, 
begin by assuming n= 1 . 

65 4. Determine the number of extended states h. A value of 

h=l is usually sufficient. 

5. Apply the new algorithm to the system. 
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6. Run a closed loop simulation or hardware test in order to 
tune the controller and observer bandwidths. 

7. Run a closed loop simulation or hardware test in order to 

tune B + , the right inverse of B. Repeat step six as 
required. 5 

Step One: The first step is to determine the number of 
control inputs p and controlled outputs q of the plant. If q>p, 
then multi-mode control should be used to make q^p,- for 
each of the j sub-controllers. It is preferable to make q^p, 
which produces a square B matrix that allows the diagonal 
elements to become tuning parameters. Note that when B is a 
square matrix, B + =B _1 . A diagonal B matrix also permits the 
new technique to be reduced to multiple SISO techniques. 

The jet engine in MAPPS, for example, has three actuator 15 
inputs that control seven performance parameter outputs. As 
a result, the jet engine controller consists of four separate 
regulators, each controlling only three outputs at a time. In 
this research, a simple form of multivariable ADRC using 
Euler integration is applied to the three-input three-output 20 
low speed regulator section and tested in simulation. This 
approach will isolate the affects that blending of multiple 
modes may have on the results. 

Step Two: The second step is to determine the high fre- 
quency gain B of the system. This matrix will drastically 25 
change for different values of relative order n. The trick is that 
n must be known in order to determine B, and B must be 
known to determine n, a circular argument. As a result, steps 
two and three are interchangeable and an iterative process 
may be used in finding n and B. Nevertheless, there is only a 30 
problem if both are unknown. If this is the case, use the 
identity matrix for B to first determine n and then iterate. The 
B matrix can also be tuned in step seven. In practice, B + needs 
to be initially within fifty percent of its true value and such a 35 
broad range is frequently known. However, if it is unknown or 
the system is too complex, then various system identification 
techniques can be used. 

In MAPSS, the control signals are scaled to produce the 
proper units for each actuator input allowing each control 40 
signal to be within the same relative range. Thus, a logical 
starting point for B in the low speed regulator is the identity 
matrix. 

Step Three: The third step is to determine the relative order 
n of the plant. The overall structure of the observer and 45 
controller depends on n which may or may not be the actual 
order of the system, depending on which dynamics are domi- 
nant. The idea is to find n=l , 2, or 3 that produces the smallest 
H, the smallest control signal, or the best closed loop results. 
Sometimes n is known or can be derived from a model of 50 
physical relationships. If none of the above techniques work, 
the last resort is trial and error. In this situation, start by 
assuming the system is first order and complete the remaining 
steps. Then assume the system is second order and repeat the 
proces s to see if the results are better. Then try third order, etc . 55 
Lower order is usually better. 

Another consideration is that n can be determined for each 
plant input-output pair. For a particular output, the input 
yielding the lowest order with the highest gain is the most 
direct form of control and therefore should be used. 60 

Since not all of the engine’s states are measurable in 
MAPPS, the model for low speed regulation is represented as 
a nonlinear input-output vector function. Without explicit 
knowledge of system order, the simplest and lowest order 
case is first attempted. 65 
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When a 3x3 matrix B is used to approximate the actual high 
frequency gain B 0 , the signal H is defined as 

HsF(Y,U,t)-B 0 U (217) 

The system then reduces to a form that has distinct terms to 
represent any internal or external dynamics and an instanta- 
neous input. 

Y~H+BU (218) 

After running simulation tests at higher orders, first order was 
found to be sufficient for MAPPS. This also makes sense 
since the CLM is represented as a first order state space 
equation and the actuator dynamics are fast enough to be 
neglected. 

Step Four: The fourth step is to determine the number of 
extended states h. This affects the overall structure of the 
ESO. For ADRC, select h=l , 2, or 3 based on the system type 
of the disturbance H or an external disturbance. The assump- 
tion h=l will suffice in most cases and therefore is used in the 
remaining examples for the sake of clarity. Similar is true for 
determining the number of extended states m in the control 
law when using generalized PID. 

Step Five: The fifth step is to apply the new algorithm to the 
system. The overall configuration is shown in FIG. 32. The 
structure of the observer and controller depends on the integer 
values selected for n and h. The most frequently used cases 
are now explicitly given for multivariable ADRC with Euler 
integration. In general, U is a pxl vector, Y is a qxl vector, B 
is a qxp matrix, and Lj and K . are qxq matrices. 

When n=l and h=l, the ESO equations become 


x 1 =Y 2 +z 1 (y-Y 1 )+ J g£/ 


X 2 =Z 2 (7-Y 1 ) 

(219) 

and the controller is represented by 


U=B + (K p (Y*-Jt l )-Jt 2 ) 

(220) 


By applying (219) and (220), the control configuration is then 
shown in FIG. 33. Notice that since the input to B + is essen- 
tially BU, it is used as an input to the ESO instead of multi- 
plying U by B. In doing so, there is only one matrix to adjust 
containing elements of B and it acts to scale the plant and 
allow the rest of the algorithm to be designed for a unity gain 
plant. 


When n=2 and h=l, the ESO equations become 


X^+L^Y-Jt,) 


X 2 =1 3 +Z 2 (7-X 1 )+5C/ 


X 3 =Z 3 (7-X 1 ) 

(221) 

and the controller is represented by 


U=B + (K p (Y*-X l )-KjC 2 -X 3 ) 

(222) 

By applying (221) and (222), the control configuration is then 
shown in FIG. 34. 

The first order plant in (2 1 8) used for low speed regulation 
is represented by state equations where the extended state X 2 
is assigned to track the general disturbance H 

x 1 =x 2 +bu 


x 2 =h 

(223) 


and the 3x1 state vectors are defined as X=[X 1 ,X 2 ] r =[Y r , 
H T ] T . By also defining the estimated states as X=[X 1 r ,X 2 r ] r , 
an ESO is designed from (223). 

±=% 2 +L i (Y-X i )+BU 
X 2 =Z 2 (7-Y 1 ) 


Y„=F(Y lr U,t) 


(216) 


(224) 
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A disturbance rejection control law is then implemented 

U=B + (U 0 -X 2 ) (225) 

to decouple the plant, reducing it to three parallel integrators 
at low frequencies 5 

Y~U 0 (226) 

whereby a simple proportional control law is applied. 

Uq = K 1 {Y*-X 1 ) (227) 10 

The entire algorithm consists of (219) and (220) and is shown 
implemented in FIG. 33. 

X^+BU+L^Y-X,) ^ 

X 2 =L 2 (Y-X{) 

U=B + (Z 1 (y*-^ 1 )-X 2 ) (228) 

The original jet engine regulators incorporate PID control- 20 
lers that are subject to integrator windup because the integra- 
tor input is a function of the controller error R-Y where 
convergence to zero is affected by plant saturation. An inter- 
esting benefit of replacing these regulators with ADRC is that 
the integrators within ADRC do not wind up because their 25 
inputs are a function of the observer error Y-X l where con- 
vergence to zero is not affected by plant saturation. As a 
result, additional anti -windup mechanisms are no longer 
required. 

Step Six: The sixth step is to run the closed loop system in 30 
order to tune the controller and observer bandwidths. In gen- 
eral, L . and K ; . are qxq matrices. However, when L f and K j are 
selected as diagonal matrices, the ADRC algorithm reduces to 
a set of SISO controllers, one to control each output. An 
example is given in FIG. 35 for a three output system. A 35 
tuning procedure for this type of configuration is proposed. 

When n=l and h=l, the resulting observer and controller 
gain matrices become 

Z, 1 =diag(2(0 <>jl ,2o) 05 2, . . . , 2oo 0>? ) 40 

Z 2 =diag(( 0 o4 2 , <d 0?2 2 , . . . , (v/) (229) 

^=diag((o c4 ,co c?2 , . . . co^) (230) ^ 

When n=2 and h=l, the resulting observer and controller 
gain matrices become 
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TABLE VIII 

NORMALIZED SETTLING TIMES FOR VARIOUS ORDERS 
n+h 12345 

t„ 3.9124 5.8342 7.5169 9.0842 10.5807 

Since a step is the fastest possible profile, the minimum 
settling time of the system for a given bandwidth then 
becomes 

t=tjto c (233) 

When faced with time specifications, this can serve as a 
starting point for tuning oo c or it can determine if a solution is 
even feasible. 

When a profile with settling time X is used, the total settling 
time of the system is then approximated by 

t,~ynt s (234) 

A procedure to tune co c . and co 0 < for each output i is now 
given, dropping the i subscript for the sake of clarity. The idea 
is to set the controller bandwidth as high as possible. A profile 
is then typically used to achieve a slower settling time or to 
meet control signal constraints, but a step is used if the fastest 
possible response is desired without regard to exact trajec- 
tory. 

1 . Set oo c using (233) according to initial specifications. 

2. Set cl) o =2~10(jl) c as a rule of thumb. The exact relation 
will be based on the proximity of the desired closed loop 
bandwidth to the dominant poles of the system, resonant 
frequencies, and noise. 

3. Run the closed loop system and increase oo c and oo 0 
together to a point just before oscillation appears in the 
control signal. 

4. Adjust the relation between oo c and oo 0 to meet the design 
specifications for noise level and disturbance rejection. 

In MAPPS, all three observer bandwidths are set equal, 
i=oo 0 2 =a) o, 3 =a) o> f° r the sake of simplicity and proof of 
concept. This makes the observer gain matrices a function of 
one parameter. 

L 1 =2o) 0 I 3 ,L 2 =2co £> 2 I 3 (235) 

All three controller bandwidths are also set equal, oo c 1 =co c 
2 =(d C; 3 =( jo c , for the same reason, making the controller gain 
matrix a function of one parameter as well. 


L 1 =diag(3(0 Ojl ,3(0 o?2 , . . . , 3(0 o J 
L 2 =diag(3 ( 0 o : 2 ,3 ( 0 o2 2 , . . . ,3co 0 ^ 2 ) 

L 3 =diag(( 0 o4 3 ,( 0 o>2 3 , . . . ,a 0 J) (231) 

K p =ding(to c /,a c /, . . . ,co^ 2 ) 

^=diag( 2(0^, 2(0^, . . . ,2co c ^) (232) 

Preliminary system identification is often unnecessary 
because the only design parameters are oo c and b and they have 
a direct impact on the bandwidth and overshoot of the output, 
meaning they can easily be adjusted by the user. In practice, 
1/b needs to be initially within fifty percent of its true value, 
the total inertia in a second order system, and such a broad 
range is frequently known. 

When a step input is applied to the system, normalized 
settling times for co c =l are shown in Table VIII. 


K^I, (236) 

Step Seven: The seventh step is to run the closed loop 
system in order to tune B + . Although B is determined using 
system identification techniques in step two, it acts as a con- 
trol signal gain that directly affects overshoot of the closed 
loop system. As a result, B can also be tuned by adjusting its 
elements to the point just before overshoot appears in the 
output. It is preferable to have as many inputs as outputs to 
facilitate tuning. This is demonstrated by first expanding the 
vector products U=[u l5 . . . , u^andb^b^, . . . , b^] 7 * in the 
following n th order input-output equations. 


yf^n+biU (237) 


yW=f q+ b q U 
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to produce: 

y[ n) =fi+ b lt i mi +- ... + ( 238 ) 

=fl+ b 2 ,l Mi + ... 4- b^M^ 

10 

In the i** state equation, the input u z is used to control the 
output y i9 and the rest of the inputs are combined into a new 
disturbance h z . 

hi = fi + b[2 u 2 + ■ ■ ■ + bl,q u q (239) 15 

^1 = h + ^ 2,1 w 2 + ^ 2 , 3 W 2 + ■ ■ ■ + b 2 , q U q 

hq = fq + b q ,\U\ + . . . 4- bq-l f q-\Uq-l 20 

The plant is rewritten where only the diagonal elements of B 
are considered. 

25 

y { ? = hi+b U iui (240) 

y { 2 = h 2 + b 22 u 2 

: 30 

y q ^ — hq + b qq li q 

The inverse of each diagonal element then becomes a tuning 
parameter of an individual SISO control loop. 

5" 1 =diag(6 1 " 1 ,62 _1 , . . . ,b 9 _1 ) (241) 

An example is given in FIG. 36 for a three-output system. 

Since each jet engine regulator has three inputs and three 
outputs and the system parameters are unknown, the plant 40 
used for low speed regulation was represented by (240) and 
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engine operating within safe limits for several thousand flight 
cycles. With repeated use, the engine components wear and 
performance is degraded. For example, turbine blades erode 
and clearances open up. In order to achieve the same level of 
thrust as a new engine, a deteriorated engine must run hotter 
and/or faster. This shift from nominal operation increases 
with use, and eventually reaches the point where performance 
can not be maintained without compromising the safety of the 
engine or the life of its components. The degradation in per- 
formance is simulated in MAPSS by adjusting ten health 
parameters. 

In most turbofan engines, thrust is calculated as a function 
of regulated and non-regulated variables, since it cannot be 
directly measured. Although regulated variables are main- 
tained at their set points regardless of engine degradation, 
non-regulated parameters shift from their nominal values 
with deterioration. As a result, the closed-loop performance 
of the current model-based controller suffers as the engine 
wears. One of the objectives here is to control the transient 
thrust response of a deteriorated engine, making it behave as 
close to a new engine as possible. 

Gas path analysis is a diagnostic technique that is used to 
estimate and trend health parameters by examining shifts in 
component health based upon gas path sensor measurements, 
i.e. pressures, temperatures, rotor speeds, and the known 
aero-thermodynamic relationships that exist between them. 
The health parameters follow an average degradation profile 
over the life span of the engine 

Pi=Vi( 1 -exp {-bit^+cjeff (242) 

where a z , b z , and c. are constants for each health parameter and 
^represents the physical age of the engine in effective flight 
cycles rather than it chronological age. The initial exponential 
rise is intended to simulate rub-in and new engine deteriora- 
tion mechanisms. As the engine ages, the health parameter 
degradation tends to become more linear. 

The component degradation values in percent resulting 
from health parameter changes are shown in Table IX. They 
reflect moderate to beyond severe degradation such as what 
might occur when the engine is due for an overhaul or when 
the engine is placed in a harsh desert environment. 


TABLE IX 


DEGRADATION VALUES DUE TO HEALTH PARAMETER CHANGES 


Flight Cycles 

Fan LPC HPC HPT LPT 


teff 

Flow % 

% 

Flow % 

% 

Flow % 

% 

Flow % 

% 

Flow % 

% 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3000 

-2.04 

-1.50 

-2.08 

-1.46 

-3.91 

-2.94 

+1.76 

-2.63 

+0.26 

-0.538 

3750 

-2.443 

-1.838 

-2.56 

-1.748 

-6.448 

-4.555 

-1.96 

-2.925 

+0.3 

-0.673 

4500 

-2.845 

-2.175 

-3.04 

-2.035 

-8.985 

-6.17 

+2.16 

-3.22 

+0.34 

-0.808 

5250 

-3.248 

-2.512 

-3.52 

-2.323 

-11.52 

-7.785 

+2.37 

-3.515 

+0.38 

-0.943 

6000 

-3.65 

-2.85 

-4 

-2.61 

-14.06 

-9.4 

+2.57 

-3.81 

+0.42 

-1.078 


the diagonal elements of B _1 were tuned. The identity matrix 
was selected as a starting point. The relative signs of each 
diagonal element were next determined, followed by magni- 
tude adjustment. 

Jet Engine Control Simulation Results 

Turbofan engine performance varies from engine to engine 
due to manufacturing tolerances and deterioration from 
extended use. Even though degradation may eventually 
require an engine to be overhauled as limits are reached, the 
engine control system should be robust enough to keep the 


The level of degradation is characterized by effective 
cycles X eff where zero cycles is a new engine with no degra- 
60 dation, 3000 cycles is moderate degradation, 4500 cycles is 
heavy degradation, and 5250 cycles is severe degradation. 

Test operating points were selected to cover a large portion 
of the entire MAPSS flight envelope and most of the subsonic 
range. They are shown in Table X. Test point #1 represents 
65 ground idle conditions where the pla is stepped from 2 1 to 35 
for takeoff. The rest of the test points represent the majority of 
subsonic power conditions. 



US 8,060,340 B2 


72 


71 


TABLE X 


TEST OPERATING POINTS WITHIN THE MAPSS ENVELOPE 

Pt. 





Op 






1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

alt 

OK 

20K 

20K 

36,089 

36,089 

36,089 

20K 

40K 

40K 

40K 

xm 

0 

0.5 

0.8 

0.5 

0.8 

1.0 

0.3 

0.3 

0.5 

0.8 

pla 

21-35 

30-35 

30-35 

30-35 

30-35 

32-37 

30-35 

30-35 

30-35 

30-35 


The redesigned low speed regulator in (228) was digitized 
using Euler integration. MAPSS is a multi-rate simulation 
package where the engine sample time is fixed at 0.0004 
seconds and the controller sample time is fixed at 0.02 sec- 15 
onds. For proof of concept, the new ADRC controller and the 
supplied nominal controller were simulated at each of the first 
three operating points in Table X. The results from the nomi- 
nal controller are used as a reference to compare the perfor- 
mance of the ADRC controller with. The goal here was not to 20 
show that one controller is better in performance over the 
other but merely that they are comparable in performance and 
ADRC is very simple to design, especially since the exact 
method of tuning the nominal controller is not known. The 
new ADRC controller was then simulated at the next three 25 
operating points in Table X to show how it is able operate over 
a substantial range of the low speed regulator. 

All simulations were conducted for each of the six levels of 
degradation in Table IX. As shown in Table XI, they are 
labeled as Run 1 through Run 6 in each simulation. 


TABLE XI 


DEGRADATION TEST RUNS 

Run# 

1 

2 

3 

4 

5 

6 

V 

0 

3000 

3750 

4500 

5250 

6000 


Although the high speed regulator and the other two fan 4Q 
safety regulators were not tested, similar performance is 
expected. The results are shown in 3700 of FIG. 37 and in 
3800 of FIG. 38 for test point #1. Note that the trajectories at 
different levels of degradation are virtually indistinguishable 
from one another. Other test point yielded similar results with 45 
no change in the new controller parameters. 

Although the ADRC controller responded a bit faster and 
with less overshoot to the change in demand levels than the 
nominal controller, the real significance is in the simplicity of 
design of the new controller and how it was able to control 5Q 
engine thrust without being affected by degradation over a 
wide range of operation. The design procedure of the nominal 
controller basically involves running the CLM at several 
operating points to calculate a set of gains from Bode and 
Ny qui st arrays at each operating point. The eighteen gains are 5 5 
each scheduled by six parameters, amounting to a total of 1 08 
possible adjustments that can be made when configuring a 
single regulator on an actual engine. During the simulations, 
these gains change by as much as 200 percent. 

In contrast, the five ADRC gains remained constant 6Q 
throughout all simulations. 

(O c =8, (i) 0 =16, 5 _1 =diag(0.2, -0.5, -0.5) (243) 

There was no scheduling. Each gain was quickly tuned on the 
CLM just as if it would be on an actual engine. The engine was 65 
then simulated at multiple operating points to verily the per- 
formance of the new controller. 


Preliminary results of these simulation tests on a rather 
complex turbofan model show the power of the dynamic 
decoupling method proposed here. Mathematical models are 
often inaccurate when representing nonlinear multivariable 
systems. Gain scheduling helps in this area, but makes tuning 
even worse than it was before. Where modern multivariable 
control schemes are limited, this approach appears well suited 
for complex nonlinear systems with incomplete model infor- 
mation. The ultimate goal is to offer a degree of tunability to 
account for variations between engines without sacrificing 
performance, while being robust enough to withstand slow 
degradations from aging or damage. 

Health Monitoring and Fault Detection using the Extended 
State Observer 

This research combines the unique concept of design 
model disturbance estimation with health monitoring and 
fault diagnosis. The tools developed in the previous descrip- 
tion can be directly applied to health monitoring with mini- 
mal model information. The unique application uses the ESO 
as a disturbance estimator with minimal plant information to 
estimate system dynamics and disturbances. The estimated 
disturbance is then used for health and fault diagnosis. Most 
dynamic health and fault monitoring estimators require sig- 
nificant model information to work effectively. Since the ESO 
use a simple design model that works on a wide variety of 
plants, estimator design can be reduced to single parameter 
tuning. In 4900 of FIG. 49, the general concept of fault 
monitoring is illustrated. 

Decomposing Disturbances for Fault Diagnosis 

The following gives a more detailed explanation of the 
concept of using the ESO for health and fault monitoring. A 
wide range of input, u, output, y, systems can be described by 
the differential equation in Han’s canonical model form 

■ ■ ■ ,y (ri ~ i ^,w)+bu. (244) 

Here y (w) represents the n th derivative of y, where f is a 
lumped nonlinear time varying function of the plant dynamics 
as well as the external disturbance w. Based on the input 
output data, detail 5000 of FIG. 50 illustrates the concept of 
generating the unknown dynamics of f from the input output 
characteristics. Once estimated, f can be analyzed for health 
diagnosis, fault detection and performance analysis. 

The unknown portion, f, contains modeling inaccuracies, 
f m , effects of faults jy, and external disturbances f d . Detail 
5100 of FIG. 51 illustrates how f can encapsulate a number of 
effects related to health monitoring and fault detection where 
f is lumped unknown dynamics, f m is unmodeled dynamics, 
f d is unmodeled external disturbances, f s is static nominal 
plant inaccuracies, f t is time varying plant degradations, jyis 
faults due to large model structural changes, f p is time varying 
model parameters and f h is health degradation. 

For the most part, literature has approached the unknown 
plant effects that compose j separately. In each case, assump- 
tions are made that the other effects are negligible. Likewise, 
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an overarching framework to integrate the problem of control, 
fault diagnosis and health monitoring and accommodation is 
an open problem for research. 

Disturbance Estimation for Control Health Diagnosis 

Recently large saving has been provided by six sigma 5 
techniques to actively monitor the health of a control loop 
performance by (C. McAnamey and G. Buckbee, “Taking it 
to the boardroom: Use performance supervision information 
for higher-level management decisions,” In Tech. ISA , Jul. 31, 
2006). Using the disturbance estimation concept for health to 
diagnostics in combination with closed loop control is an 
effective means to provide health diagnosis without extensive 
model information. 

Once the model is formulated in the input, u, output, y, and 
disturbance, f, formulation, the control problem can be refor- 15 
mulated. 

Canceling the unknown disturbance and dynamics begins 
with estimating f. The main idea is to use input output data 
and minimal dynamic information to estimate j and cancel it. 

20 

Kf (245) 

Once f is estimated, the disturbance is rejected to behave 
like the forced design model plant, with a new input, u 0 . 

find u s.t. (246) ^5 

At this point, the unknown disturbances and plant dynam- 
ics have been removed and a conventional controller based on 
the design model can be designed so the output, y, follows the 
reference, r. 30 

find u 0 s.t. y-»r. (247) 

This overview of concepts suggests that there are three 
independent mathematical expressions that solve the control 
problem: 1) the estimation law (245), 2) the rejection law 35 
(246), and 3) the nominal control law (247). This division is 
illustrated in 5200 of FIG. 52 . 

Most control paradigms include the estimation and rejec- 
tion laws lumped together in the control law. Since f is the key 
to this control paradigm, this research investigates the active 40 
estimation of f for health and fault monitoring. 

Health Monitoring by Disturbance Estimation 

The general methodology to apply the Extended State 
Observer to the health monitoring problem follows: 

1. Determine the appropriate coupled input and outputs. 45 
For single input and single output systems this step is not 
required. However, for effective estimation of multi- 
input multi -output systems, each input needs to be pared 
with an output that is dynamically linked in some man- 
ner. Cross coupling between separate input-output pairs 50 
is included in the disturbance estimate of each pair. This 
way each input-output pair can be considered indepen- 
dently. 

2. Determine the order of each input output coupling. The 
order can be determined from intuition about the physi- 55 
cal process or by trial and error. 

3. Build a matching Extended State Observer to estimate 
states and disturbances 

4. Select tuning parameters for stable output tracking. 

5. Determine nominal conditions for the estimated distur- 60 
bance, f. 

6. Monitor the estimated disturbance for variations from 
the nominal conditions. 

7. If model information is known, specific fault informa- 
tion can be extracted through the estimation of f. Since 65 
the dynamic information is estimated in f, this can usu- 
ally be composed of an algebraic equation. 
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The systems, methods, and objects described herein may 
be stored, for example, on a computer readable media. Media 
can include, but are not limited to, an ASIC, a CD, a DVD, a 
RAM, a ROM, a PROM, a disk, a carrier wave, a memory 
stick, and the like. Thus, an example computer readable 
medium can store computer executable instructions for one or 
more of the claimed methods. 

What has been described above includes several examples. 
It is, of course, not possible to describe every conceivable 
combination of components or methodologies for purposes of 
describing the systems, methods, computer readable media 
and so on employed in scaling and parameterizing control- 
lers. However, one of ordinary skill in the art may recognize 
that further combinations and permutations are possible. 
Accordingly, this application is intended to embrace alter- 
ations, modifications, and variations that fall within the scope 
of the appended claims. Furthermore, the preceding descrip- 
tion is not meant to limit the scope of the invention. Rather, 
the scope of the invention is to be determined only by the 
appended claims and their equivalents. 

All documents cited herein are, in relevant part, incorpo- 
rated herein by reference; the citation of any document is not 
to be construed as an admission that it is prior art with respect 
to the present invention. 

While the systems, methods and so on herein have been 
illustrated by describing examples, and while the examples 
have been described in considerable detail, it is not the inten- 
tion of the applicants to restrict or in any way limit the scope 
of the appended claims to such detail. Additional advantages 
and modifications will be readily apparent to those skilled in 
the art. Therefore, the invention, in its broader aspects, is not 
limited to the specific details, the representative apparatus, 
and illustrative examples shown and described. Accordingly, 
departures may be made from such details without departing 
from the spirit or scope of the applicant’s general inventive 
concept. 

The invention claimed is: 

1. A computer-implemented method for controlling a 
velocity within a dynamic system, comprising: 

specifying a velocity value that is defined as v(t)=f(t)+bu 
(t), where f(t) represents the combined effects of internal 
dynamics and external disturbance of the plant, u(t) is a 
control signal, and b is a constant to an approximate 
value; 

converting the velocity value into a first order state space 
model; 

estimating the value of f(t) by a linear extended state 
observer, which is a function of a single performance 
parameter; 

canceling the effect of ;f(t) on the velocity by the estimate 
from the linear extended state observer. 

2. The method according to claim 1, further comprising: 

assigning eigenvalues of the observer as a function of a 

single tuning parameter; 

reducing the velocity value to an approximate integral 
plant, 

v(t)=(f(t)-z 2 (t))+u 0 (t)=u 0 (t),and 

controlling the approximate integral plant via u 0 (t)=kp(r 

(t)-Zi(t)). 

3. The method according to claim 2, further comprising: 

creating an approximate closed-loop transfer function with 

no finite zeros; 

solving each of the controller gains to correspond to the 
single tuning parameter by a model transfer function; 

setting each of the controller gains to correspond to the 
single tuning parameter. 
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4. The method according to claim 1, wherein the dynamic 
system is a web processing system. 

5. The method according to claim 4, wherein the velocity 
value is at least one of a carriage roller velocity, an exit-side 
roller velocity, and a process-side roller velocity. 

6. A method of monitoring the health of a system compris- 
ing: 

receiving sensor data from a sensor operably connected to 
the system to measure an output (y); 

storing an input (u) applied to the system that is coupled to 
the output (y); 

processing the input (u) and the output (y) using an 
extended state observer, the extended state observer 
designed to estimate at least one state and a disturbance 
(f) associated with a model that couples the input (u) 
with the output (y); and 

comparing the disturbance (j 1 ) against a nominal distur- 
bance to create a variance. 

7. A method of claim 6, further comprising a disturbance 
rejection controller that is adapted to reject the disturbance 

CD- 

8. A method of claim 7, wherein the disturbance rejection 
controller rejects the disturbance such that the disturbance 
and plant dynamics associated with the system are removed 
from the system. 

9. A method of claim 8, wherein a controller is applied to 
the system such that the output (y) follows a reference (r). 

10. A method of claim 6, further comprising: 

modeling a fault in the system associated with a character- 
istic variance; and 

determining the fault by comparing the variance with the 
disturbance. 

11. A method of 10, further comprising: 

rejecting the disturbance from the system by feeding the 
disturbance into a disturbance rejection controller to 
create a new command signal adapted to reject the dis- 
turbance. 

12. A method of claim 6, wherein the disturbance (j 1 ) com- 
prises, unmodeled dynamics ( f m ), unmodeled external dis- 
turbances (f d ), status nominal plant inaccuracies (f s ), time 
varying plant degradations (f t ), faults due to large modal 
structural changes (jy), time varying model parameters (f ), 
and health degradation ( f h ). 

13. A health monitoring system for a plant comprising: 

a computer adapted to receive a signal from a sensor oper- 
ably connected to the plant; 

an output (y) of the plant that is determined from the signal; 

an input (u) is determined from a control signal applied to 
the plant; 

an extended state observer adapted to estimate a state and a 
disturbance (J) associated with the input (u) and the 
output (y); 

a monitor adapted to compare the disturbance (j 1 ) with a 
nominal disturbance; and 

a model of specific fault information that is matched to the 
monitor to determine a specific fault. 

14. A system of claim 13, further comprising: 

a disturbance rejection controller that is adapted to gener- 
ate a rejection signal that removes unknown distur- 
bances and plant dynamics. 

15. A system of claim 14, further comprising: 

a controller wherein the controller outputs a new input (uo) 
such that the output (y) follows a reference (r); and, 

a control signal applied to the plant via an output from the 
computer comprising the sum of the new input(uo) and 
the rejection signal. 
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16. A computer- implemented method to generalize an 
extended state observer, comprising: 

representing a plant with a continuous-time differential 
equation of an nth order plant, y ( ” )= f(y, y, . . . , y( (w_1) ), 
5 w,t)+bu, where y(n) denotes the nth derivative of y, u is 

a control signal, and b is an estimate of a value; 
constructing an n+h order state space model of the nth 
order plant using h-cascaded integrators to represent a 
disturbance f and its h derivatives; 

1 0 discretizing the state space model by applying one of Euler, 

zero order hold, or first order hold methods; 
creating a predictive discrete estimator from the discretized 
state space model; 

15 creating a current discrete estimator from the discretized 
state space model; 

implementing the current discrete estimator in a computer 
component that receives a signal from a sensor that is 
operably connected with the plant; and 
20 processing the signal using the discrete estimator to esti- 
mate the value of the disturbance f and its h derivatives. 

17. The method according to claim 16, further comprising: 
discretely parameterizing estimator gains, such that imple- 
mentation of the current discrete estimator is based upon 

25 a single tuning variable. 

18. The method of claim 16, wherein the plant has more 
than one input and more than one output. 

19. A method to enhance performance of an active distur- 
bance rejection controller by providing transient tracking 

30 control, comprising: 

representing a plant with a continuous-time differential 
equation of an nth order plant, y (w) =f(y, y, . . . , y (w-1) , w, 
t)+bu y(II)=f(y,:/, . . . ,y(II-I), w, f)+bu, where yen) 
denotes the nth derivative of y, u is a control signal, and 
35 b is a known value; 

constructing an extended state observer to estimate gener- 
alized disturbance, output y, and n-1 derivatives of the 
output; 

applying a disturbance rejection control law to cancel the 
40 generalized disturbance using an estimated disturbance 
value from the extended state observer; 
reducing the plant to n cascaded integrators; 
applying a point-to-point control law to the reduced plant 
to form a desired closed-loop transfer function; 

45 adding an inverse of the closed-loop transfer function to a 

reference input of the controller to form a new closed 
loop transfer function equal to one, or its relative order 
equal to zero; 

implementing the extended state observer, the disturbance 
50 rejection control law, and the point-to-point control law 
in a computer component that receives a signal from a 
sensor that is operably connected with the nth order 
plant; 

processing the signal using the extended state observer, the 
55 disturbance rejection control law and, the point-to-point 
control law to estimate a control signal u; and, 
applying a control command representing the control sig- 
nal u to the nth order plant. 

20. The method according to claim 19, further comprising: 
60 applying a tracking control law comprised of feed forward 

terms directly to the reduced plant to form a desired 
closed loop transfer function equal to one, or its relative 
order equal to zero. 

21. The method of claim 19, wherein the plant has more 
65 than one input value and more than one output value. 

22. A computer-implemented method to discretely imple- 
ment an extended state observer, comprising: 
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representing a plant with a continuous-time differential 
equation of an nth order plant, y ( " )= (y,y,w,t)+bu where f 
is a function of an internal system dynamics, an external 
disturbance w, and b is a constant; 
constructing a n+ 1 order state space model of the nth order 
plant; 

discretizing the state space model by applying one of Euler, 
zero order hold, or first order hold methods; and 
creating a predictive discrete estimator from the discretized 
state space model; 

creating a current discrete estimator from the discretized 
state space model; 

implementing the current discrete estimator in a computer 
component that receives a signal from a sensor that is 
operably connected with the plant; and 
processing the signal using the current discrete estimator to 
estimate the value of the function f . 

23 . The method according to claim 22, further comprising: 
discretely parameterizing estimator gains, such that imple- 
mentation of the current discrete estimator is based upon 
a single tuning variable. 

24. A method for designing a system to control a multiple- 
input, multiple-output system, comprising: 

discretizing a system model to describe one or more dis- 
tinct states, where each input has a distinct output and 
disturbance; 

constructing an extended state estimator from the dis- 
cretized system model; 

implementing the extended state estimator in a computer 
component that receives a signal from a sensor that is 
operably connected to the system; 
determining one or more correction terms as a function of 
a single tuning parameter; and 
utilizing the correction terms with the extended state esti- 
mator and the signal to estimate system states and 
extended states of one or more orders. 

25. The system according to claim 24, wherein the system 
is at least one of a chemical process, a mechanical process, or 
an electrical process. 

26. A method to provide health monitoring to a system, 
comprising: 

determining appropriate coupled control inputs and con- 
trol outputs; 

determining order of each input/output coupling; 
building a matching extended state observer to estimate 
states and a disturbance; 
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adjusting a value of at least one tuning parameter that 
provides stable output tracking; 
determining at least one nominal condition for an esti- 
mated disturbance; 

5 implementing the matching extended state observer in a 
computer component that receives a signal from a sensor 
that is operably connected to the system; 
processing the signal using the matching extended state 
observer to estimate the disturbance; 
to monitoring a variation between the estimated disturbance 
and the at least one nominal condition; and 
extracting fault information from the variation. 

27. A method for controlling a turbofan, comprising: 
creating a model of a portion of a turbofan system as a 

15 nonlinear input-output vector function; 

approximating a general disturbance of the modeled sys- 
tem; 

reducing the system to a second model that distinguishes 
between an instantaneous input and one or more 
20 dynamic variables to be estimated in real time; 

representing the system by one or more state vectors, 
wherein an extended state is assigned to track the general 
disturbance; 

determining a disturbance rejection control law; 

25 implementing the disturbance rejection control law on a 
computer component that receives a signal from a sensor 
that is operably comiected to the turbofan; 
utilizing the disturbance rejection control law to decouple 
the system and reduce it to one or more parallel integra- 
30 tors; 

controlling the simplified parallel integrator system; and 
generating a control signal that is applied to the turbofan. 

28. A method of adding disturbance information into the 
linear extended state observer is described, comprising: 

35 representing a plant with a continuous-time differential 

equation of an nth order plant, where the nth derivative 
of the output y(t) equals a generalized disturbance f(t) 
plus an input bu(t) where b is a constant; 
constructing a state space model of the plant; 

40 creating an extended state ob server based on the state space 

model, having correction terms that are a function of a 
single parameter; and 

adding a term to the extended state observer that is com- 
posed of the derivative of f(t) if it is known or partially 
45 known. 



